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On  Certain  of  the  More  Complex  Stress  'Distributions  in 
Engineering .  Materials.— Report  of  Committee  (Prof.  E.  G. 
GOKER,  Chairman;  Profs.  L.  N.  G.  FILON  and  A.  ROBERTSON, 
Secretaries;  Prof.  A.  BARB,  Dr.  C.  CHREE,  Dr.  GILBERT  COOK, 
Prof.  W.  E.  DALBY,  Sir  J.  A.  EWING,  Mr.  A.  E.  FULTON,  Mr.  J.  J. 
GUEST,  Dr.  B.  P.  HAIGH,  Profs.  Sir  J.  B.  HENDERSON,  C.  E. 
INGLIS,  F.  C.  LEA,  A  E.  H.  LOVE,  and  W.  MASON,  Sir  J.  E. 
PETAVEL,  Dr.  F.  ROGERS,  Dr.  W.  A.  SCOBLE,  Mr.  R.  V.  SOUTH- 
WELL, Dr.  T.  E.  STANTON,  Mr.  0.  E.  STROMEYER,  Mr.  J.  S. 
WILSON. 

Introduction. 

THE  Committee  submit  as  their  Report  the  following  contributions  embodying  : 
(1)  A  review  of  recent  advances  in  the  special  domains  of  (a)  Stress  Concentra- 
tions, (6)  the  Theory  of  Plane  Stress  in  perforated  plates;  (2)  special  researches 
carried  out  by  members  of  the  Committee  during  the  past  year;  (3)  an  account 
by  Prof.  B.  P.  Haigh  of  his  theory  of  failure  of  elasticity. 

I.  Stress   Concentrations   due  to   Notches  and  like    Discontinuities.     By  Prof. 

E.    G.  Coker,  F.R.S.,  and  Dr.  Paul  Heymans. 

II.  The  Distribution  of  Stress  in  a  Flanged  Pipe.     By  Prof.  Gilbert  Cook,  D.Sc. 

III.  On    Stresses    in    Multiply-connected    Plates.     By   Prof.    L.    N.    G.    Filon, 
F.R.S. 

IV.  Stress  Concentrations  in  Theory  and  Practice.     By  A.  R.  Griffith,  M.Eng. 

V.  The  Strain  Energv  Function  and  the  Elastic  Limit.     By  Prof.  B.  P.  Haigh, 

D.Sc.,   M.B.E.  " 

VI.  Alternating    Combined    Stress    Experiments.     By    W.    Mason,    D.Sc.,    and 

W.  J.  Delaney,  B.Eng. 

VII.  On    some    Problems    relating    to    the    Design    of    High-speed    Discs.     By 
R.  V.  Southwell,  M.A. 

VIII  On  the  Stability  of  a  Rotating  Shaft,  subjected  simultaneously  to  End 
Thrust  and  Twist.  By  R.  V.  Southwell,  M.A.,  and  Barbara  S.  Gough. 

IX.  The  Stresses  in  Cylinders  and  Pipes  with  Eccentric  Bore.  •  By  G.  B. 
Jeffery,  M.A.,  D.Sc. 

The  Committee  ask  for  re-appointment,  with  a  grant  of  100Z. 

I. 

Stress  Concentrations  due  to  Notches  and  Like  Discontinuities. 

By  Prof,  E.  G.  COKEB,  F.R.S. ,  and  Dr.  PAUL  HEYMANS, 
University  College,  London. 

IT  is  well  known  that  discontinuities  in  loaded  members  cause,  in  general, 
severe  stresses  in  their  neighbourhood,  and  therefore  do  not  allow  full  economy 
of  material  to  be  realised  in  cases  where  such  discontinuities  are  inevitable  in 
practical  construction.  In  the  majority  of  such  cases  it  is  difficult  and  some- 
times impossible  to  calculate  the  maximum  stresses  produced,  and  experience 
in  the  behaviour  of  similar  members  is  then  usually  relied  on,  in  the  absence 
of  an  accurate  knowledge  of  the  stress  distribution  which  exists.  The  stress 
distribution  caused  by  discontinuities  has  also  become  of  importance  in  recent 
years  in  the  testing  of  materials,  and  a  considerable  field  of  inquiry  has 
developed  on  the  behaviour  of  materials  subjected  to  impulsive  loads,  which 
produce  stresses  at  definitely  shaped  discontinuities  sufficient  to  cause  rupture 
at  the  chosen  section,  and  in  this  way  a  classification  of  materials  is  obtained 
which  has  proved  of  great  value  in  engineering  practice,  although  with  our 
present  knowledge  such  testing  operations  can  only  be  regarded  as  of  an 
empirical  nature. 
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.In,  ^he.  present  report  *a  brief  account  will  be  attempted  of  some  experimental 
in^stigp'bi^f  Avhi£B  Jiavje  been'made  on  discontinuities  by  photo-elastic  methods 
descYlb'ed*m«*pfevtolis.se,porU>(J),  whereby  stress  distributions  have  been  deter- 
mined sufficiently  completely  to  allow  of  a  fairly  accurate  value  being  assigned 
to  the  maximum  stresses  experienced  under  given  loads.  In  all  these  cases  the 
member,  a  flat  bar,  has  been  subjected  to  plane  stress,  and  the  results  are, 
therefore,  only  applicable  strictly  to  this  kind  of  stress  except  in  some  cases 
which  will  be  noted. 

An  interesting  case  of  a  discontinuity  of  a  simple  kind  arises  where  it  is 
necessary  to  alter  the  cross  section  of  a  member  in  a  definite  manner,  as  for 
example  in  a  flat  tension  member  or  a  cylindrical  rod.  This  case  occurs  in  a 
variety  of  instances  in  practice,  and  the  salient  fact  in  connection  with  this 
type  of  discontinuity  is  the  increase  of  stress  which  occurs  near  the  join  of  the 
smaller  to  the  larger  section  when  the  contours  are  formed  by  circular  arcs, 
which  are  tangential  to  the  boundaries  of  the  smaller  section  and  intersect  the 
larger  at  an  angle. 

The  maximum  stresses  reached  are  shown  experimentally  not  to  occur  exactly 
at  the  joins  (2),  but  slightly  beyond,  at  places  where  a  tangent  to  the  contour 
makes  a  small  angle  to  the  axis  of  the  specimen.  Thus,  for  example,  in  a  case 


(6976.A) 


where  the  smaller  breadth  of  a  flat  bar  is  .4488  in.  and  the  larger  0.855  in. 
connected  by  arcs  of  0.3  in.  radius,  the  maximum  stress  is  found  at  a  point 
of  the  curved  contour  where  a  tangent  line  makes  an  angle  of  about  8°  to  the 
line  of  pull,  and  the  stress  intensity  is  nearly  20  per  cent,  greater  than  the 
uniform  stress  in  the  smaller  section. 

The  distribution  in  this  member  at  the  contour  is  shown  in  the  accompany- 
ing fig.  1,  from  which  it  will  be  observed  that  a  maximum  stress  of  1,470  Ib. 
per  sq.  in.  is  reached  when  the  uniform  stress  in  the  smaller  section  is 
1,230  Ib.  per  sq.  in.  As  might  be  expected,  this  maximum  stress  increases  with 
decrease  of  radius  at  the  re-entrant  angle,  and  in  a  very  nearly  similar  contour, 
but  with  a  re-entrant  angle  of  -^  in.  radius,  an  increase  of  stress  of  57.5  per  cent. 
was  observed. 

The  type  of  contour  shown  in  fig.  1  is  probably  more  frequently  met  with 
in  British  practice  than  any  other.  It  is,  for  example,  adopted  by  the 
Engineering  Standards  Committee  as  a  standard  form  of  end  (with  a  radius  of 
one  inch  for  the  connecting  arc)  for  tension  test  members  of  boiler-plates  and 
like  material. 

As  will  be  observed  from  the  contour  stresses,  the  distribution  in  the  region 
of  the  discontinuity  is  extremely  variable  and  in  the  interior  of  the  plate  is  of 
a  complex  nature. 
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Along  the  central  line,  for  example,  there  is  a  gradual  change  .frpm  ihigh  tq 
low  tensional  stress  p  111  the  direction  of  pull,  accompanied  by.  a '  v>*?4We. ..C^QS?, 
stress  q,  which  latter  has  a  maximum  tensional  vf\iuev"ve'ry  near 'to  Ihe  section- 
where  the  width  of  the  member  begins  to  increase,  but  it  soon  changes  sign 
and  becomes  a  small  cross  compression  stress,  and  ultimately  vanishes  when' the 
stress  in  the  larger  section  becomes  a  uniform  tension.  A  map  of  the  stress-  in 
the  direction  of  the  line  of  pull  for  this  case  is  drawn  to  a  distorted  scale  both 
horizontally  and  vertically  (tig.  2),  in  order  to  show  the  distributions  at  sections 
one-tenth  of  an  inch  apart,  and  from  this  it  will  be  observed  that  the  greatest 
variations  of  stress  occur  at  the  contours,  and  in  general  within  the  region  of 
complex  stress  the  maximum  stresses  occur  at  the  sides  in  the  smaller  section 
and  along  the  axial  line  in  the  larger  section.  The  type  of  stress  distribution 
shown  here  does  not,  however,  appear  to  persist  beyond  the  elastic  limit  in 
ductile  material,  and  the  variations  of  stress  shown  in  the  figure  probably  tend 
to  equalise,  and  hence  fracture  does  not  necessarily  begin  to  take  place  at  the 
point  on  the  contour  where  the  maximum  stress  is  indicated.  In  brittle 
materials,  however,  where  little  change  occurs  in  the  type'  of  stress  distribution, 
it  is  found  that  fracture  occurs  very  frequently  at  this  cross  section.  This 
latter  result,  which  is  often  ascribed,  and  generally  erroneously  so,  to  imperfect 
centering  of  the  specimen  in  the  testing  machine,  is  therefore  more  probably 
explained  by  the  facts  of  the  distribution  observed. 

Another  fact  which  emerges  from  the  experimental  observations  is  the 
penetration  of  complex  stress  into  the  parallel  part  of  the  reduced  section,  and 
it  is  easy  to  show  for  this  case,  where  the  model  is  one  of  three-tenths  scale  of 
an  Engineering  Standards  Committee  test  bar,  that  complex  stress  occurs  for  a 
distance  of  .185  in.  within  the  parallel  part,  so  that  in  a  specimen  of  this  size 
and  contour,  and  a  parallel  part  0.37  in.  long,  there  is  no  pure  tension  stress 
at  any  cross  section  between  the  enlarged  ends  except  possibly  the  central  one. 
In  a  geometrically  similar  test  bar  with  connecting  arcs  of  one  inch  radius 
it  would  appear,  therefore,  that  a  total  length  of  1.23  in.  of  the  parallel  part 
is  in  a  state  of  complex  stress,  and  since  it  is  found  experimentally  that  complex 
stress  distribution  at  the  ends  does  not  vary  to  an  appreciable  extent  as  the 
length  of  the  parallel  part  increases  beyond  this,  it  may  be  inferred  that,  within 
the  elastic  limit,  a  sufficiently  short  tension  test  member  under  load  may  be  in 
a  state  of  complex  stress,  which  has  little  or  no  resemblance  to  pure  tensional 
stress  at  any  part  of  it,  and  it  is  therefore  likely  to  give  fallacious  results  for 
practical  applications. 

In  the  example  cited  above  the  limiting  length  of  the  parallel  part  for  which 
there  is  no  pure  tension  stress  across  the  sections  between  the  shoulders  is 
1.23  in.,  and  for  such  a  case  complex  stress  occurs  at  all  sections  except 
possibly  the  central  one. 

The  danger  of  accepting  tension  tests  on  short  specimens  as  representing  the 
true  behaviour  of  the  material  in  pure  tension  is  clear,  especially  as  it  is  not 
•difficult  to  show  (2)  that  in  turned  specimens  of  the  same  contour  as  a  flat  bar 
specimen  the  complex  stress  in  the  former  extends  further  into  the  gauge  length 
than  in  the  latter. 

In  another  standard  form  widely  used  in  Continental  practice  the  enlarged 
ends  are  gradually  tapered  towards  the  gauge  length  and  have  straight  line 
contours,  which  join  the  parallel  part  at  a  very  obtuse  angle.  At  this  junction 
photo-elastic  observations  show  that  there  is  a  small  amount  of  stress  concen- 
tration probably  not  more  than  10  per  cent,  greater  than  the  stress  intensity  in 
the  gauge  length  and,  moreover,  very  localised  in  extent. 

In  still  another  form,  used  by  Professor  Dalby,  mainly  to  avoid  contact 
stresses  and  local  indentations  of  extensometer  screws,  the  extremities  of  the 
gauge  length  are  defined  by  two  thin  collars  turned  on  the  specimen  with  con- 
necting arcs  to  the  main  body  of  0.04  in.  radius.  In  this  case  the  stress 
concentration  is  very  local,  and  is  rather  more  than  30  per  cent.  (3)  of  the  mean 
average  stress,  and  is  symmetrically  disposed  with  reference  to  the  collar.  It 
is,  moreover,  practically  independent  of  the  radius  of  the  fillet  within  wide 
limits,  and  in  that  respect  differs  from  the  stress  at  a  re-entrant  angle,  where 
the  radius  of  the  fillet  'is  the  main  factor  which  determines  the  stress  concen- 
tration. 

Observations  of  stress  distribution  by  photo-elastic  means  indicate,  in 
general,  that  the  stress  concentrations  described  above  tend  to  equalise  more  or 
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less  in  the  plastic  stage.  Thus  in  a  nitro-cellulose  model,  of  the  form  used  by 
the  Engineering  Standards  Committee,  it  is  observed  that  the  varied  colour 
effects,  which  mark  the  penetration  of  complex  stress  within  the  gauge  length, 
tend  to  recede  towards  the  enlarged  ends  as  the  material  passes  into  the  plastic 
state,  and  in  general  at  fracture  the  line  of  greatest  colour  variation  is  concen- 
trated near  to  the  section  where  the  curved  contours  meet  the  parallel  part  of 
the  gauge  length,  so  that  it  seems  possible  in  ductile  material  that  nearly  the 
wh9le  of  the  parallel  part  of  the  gauge  length  ultimately  comes  into  a  state  of 
uniform  stress. 

Discontinuities  in  Materials  Subjected  to  Impact  Tests. 

The  practical  necessity  of  ascertaining  the  ability  of  materials  to  withstand 
suddenly  applied  loads  has  led  to  the  introduction  of  shock  tests  of  various 
kinds,  and  these  are  generally  arranged  in  such  a  manner  that  fracture  of  the 
specimen  shall  occur,  generally,  with  a  single  application  of  a  load  applied  by 
a  falling  weight  or  swinging  pendulum.  Such  tests  are  essentially  different 

Fig.3(l.STRESS  DISTRIBUTION  AT  THE  CONTOUR 
OF  A  VEE  SHAPED  NOTCH. 

1100 
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from  static  tests  of  materials,  and  in  general  classify  them  in  a  different  manner, 
for  it  is  quite  possible  to  have  materials  which  give  similar  results  under  static 
loading  and  yet  behave  very  differently  under  impact  tests. 

Although  there  cannot  be  any  doubt  of  the  importance  of  impact  tests  in 
affording  an  element  of  value  in  determming  the  properties  of  a  material  to 
resist  stress,  yet  there  is  a  considerable  difference  of  opinion  as  to  the  exact 
nature  of  the  information  afforded,  and  this  is  hardly  surprising  when  th« 
number  of  variables  is  considered  which  enter  into  the  problem. 

The  various  machines  designed  for  impact  testing  differ  greatly  in  the  manner 
and  swiftness  with  which  the  load  is  applied,  in  the  rigidity  of  the  frame  and 
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clamping  devices  for  the  test  specimen,  while  the  form  of  the  piece  also  varies 
very  much  especially  as  regards  the  notch  cut  in  it  to  fix  the  place  of  failure. 
These  and  other  circumstances  tend  to  explain  why  it  is  so  difficult,  at  present, 
to  correlate  the  results  obtained  in  a  complete  and  satisfactory  manner. 

An  important  element  for  consideration  is  the  stress  distribution  produced 
under  load  in  various  standard  types  of  discontinuity,  and  it  is  proposed  to 
examine  some  simple  cases  here,  and  discuss  their  bearing  on  the  problem. 

It  is  evident  from  other  cases  of  discontinuity  already  described  that,  within 
the  elastic  limit,  the  contour  of  the  notch  is  of  chief  importance  as  regards 
stress  concentration,  and  that  whether  the  notch  be  subjected  to  bending  or 
tensional  stress,  the  salient  facts  will  be  much  the  same  at  the  contour.  There 
is,  however,  a  considerable  difference  in  the  experimental  difficulties,  and  very 
little  to  be  gained  in  a  preliminary  inquiry  by  examining  a  notch  under  the 
former  type  of  stress.  In  all  these  experiments,  therefore,  notches  are  cut  of 
considerable  size  in  a  very  wide  plate  of  transparent  material,  and  this  plate 
is  subjected  to  uniform  tensional  stress  in  a  testing  machine  in  order  to  find 
the  stress  at  the  notch  contour  and  across  the  principal  section  through  the  line 
of  symmetry. 

To  obtain  symmetrical  conditions  about  the  line  of  pull  a  second  notch  is  cut 
in  the  plate  at  the  opposite  edge. 

In  some  standard  forms  used  in  impact  testing  the  notch  is  of  V  form  with 
sides  inclined  at  45°,  and  with  a  definite  radius  of  curvature  at  the  apex. 


A 

FigA.MAXIMUM  STRESSES  AT  THE  APEX 
OF  A  45°  VEE  SHAPED  NOTCH  OF 
JT  FIG.3a.FORVARYING  RADII.  . 
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976.0.)             Radius  at  Apex  of  Notch  (Millimeter*) 

Such  a  notch  is  shown  on  the  accompanying  fig.  SA,  one  centimetre  deep 
with  a  radius  of  g  cm.  at  the  apex  and  cut  in  a  plate  13  cm.  wide.  Under  load 
there  is  no  great  difficulty  in  measuring  the  distribution  of  stress  at  any  point 
of  a  plate  of  this  size,  and  when  this  is  carried  out  along  the  minimum  cross 
section  perpendicular  to  the  line  of  pull,  a  measure  of  the  general  accuracy  of 
the  observations  is  obtained,  since  the  value  fp.dx  along  this  line  should  be 
equal  to  the  load  registered  by  the  testing  machine. 

The  values  of  p  and  the  cross  stress  q  in  this  case  are  shown  in  the  accom- 
panying diagram  (fig.  SB),  from  which  it  appears  that  the  stress  rises  to  a 
maximum  of  1,100  Ib.  per  sq.  in.  accompanied  by  a  cross  stress  q,  which 
is  very  small  except  near  the  ends,  where  a  maximum  value  of  about 
200  pounds  per  square  inch  is  reached  at  a  distance  of  between  0.3  to  0.4  of  a 
centimetre  from  the  edge  of  the  notch.  The  mean  stress  as  given  by  the  test- 
ing machine  is  230  Ib.  per  sq.  in.,  while  the  integrated  value  of  the  curve 
of  distribution  corresponds  with  240  Ib.  per  sq.  in.,  or  about  4.8  per  cent. 
in  excess.  Assuming  the  correctness  of  the  testing  mahcine  reading,  the  stress 
concentration  at  the  notch  is  4.78  times  the  average  stress  at  the  cross  section, 
and  for  an  angular  distance  of  about  30°  on  each  side  of  this,  the  stress  at  the 
contour  does  not  fall  below  4.40  times  of  this  value,  but  beyond  this  the  stress 
falls  rapidly  to  zero  along  the  notch  contour.  If,  however,  the  radius  of  the 
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notch  is  decreased  to  -fa  cm.,  the  maximum  stress  rises  to  nearly  six  times  the 
mean  stress,  as  fig.  4  shows,  while  if  the  radius  is  very  large,  say  10  mm.,  the 
maximum  stress  is  approximately  2.2  times  the  mean  stress. 

It  seems  evident  from  these  experiments  and  others  described  below  that 
an  impact  test  on  a  notched  specimen  affords  a  valuable  means  of  discriminating 
between  ductile  and  brittle  material,  for  whatever  be  the  radius  of  curvature 
at  the  apex  of  the  notch  in  a  ductile  material,  it  acts  primarily  as  an  indicator 
of  the  place  where  fracture  is  to  commence,  but  once  plastic  stress  begins,  its 
influence  in  maintaining  the  high  stress  concentration  observed  here  must  recede 
in  importance,  although  it  is  quite  .possible  its  form  may  be  a  factor  in  the 
final  result.  In  hard  and  brittle  materials,  however,  the  radius  at  the  bottom 
of  the  notch  is  apparently  a  dominant  factor  since  the  stress  distributions 

Fig.5a.  STRESS  DISTRIBUTION  AT  THE  CONTOUR 
OF  A  NOTCH  OF  THE.  CHARPY  TYPE. 
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observed  here  can  suffer  little  change  in  type  throughout,  and  when  due  allow- 
ance is  taken  of  static  breaking  stress  in  tension,  such  materials  must  have 
a  comparatively  low  impact  value  on  account  of  the  maintenance  of  the 
relatively  high  stress  concentration. 

As  a  discriminating  test  of  quality  an  impact  test  may  possibly  be  too 
severe,  and  good  material  may  be  rejected  by  it,  since  the  radius  of  curvature 
of  the  notch  is  a  factor  of  variable  importance  according  as  the  material  is  or 
is  not  ductile.  In  very  hard  materials  there  seems  good  evidence  for  bringing 
this  factor  into  account  in  grading  materials  by  impact  tests. 

This  is  more  particularly  emphasised  when  tests  under  different  conditions 
of  notch  form  are  considered,  as,  for  example,  in  one  of  the  Charpy  type,  in 
which  all  the  conditions  of  experiment  are  exactly  the  same  as  before,  with  a 
notch  1  cm.  deep  rounded  off  at  the  apex  by  an  arc  of  ^  cm.  radius, 
(fig.  SA).  Measuring  the  general  accuracy  of  the  experimental  values  by  a 
comparison  of  the  average  stress  intensity  of  282  Ib.  per  sq.  in.  recorded  by 
the  testing  machine  with  the  integrated  stress  intensity  curve  AB  of  fig.  SB. 
1921  Y 
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we  find  in  this  case  a  value  of  289  Ib.  per  sq.  in.,  or  2.48  per  cent,  in  excess, 
with  a  maximum  stress  of  1,000  Ib.  per  sq.  in.  at  the  notch  contour,  and  a 
concentration  stress  factor  of  3.54,  very  much  less  in  fact  than  is  afforded 
by  the  standard  notch  of  the  form  used  in  Izod  machines  and  described  above. 

This  difference  in  concentration  of  stress  is  even  more  noticeable  when  the 
notch  is  of  rectangular  cross  section  (fig.  GA)  with  angles  rounded  off  to  a 
radius  of  0.2  cm.  for  a  notch  1  cm.  deep.  In  this  case  the  mean  stress  of 
230  Ib.  per  sq.  in.  recorded  by  the  testing  machine  agrees  remarkably  well 

Fig.  6a.  STRESS  DISTRIBUTION  ATTHE  CONTOUR 
OF  A  STANDARD  RECTANGULAR  NOTCH  WITH 
ROUNDED  CORNERS.    ^—-  70O , 
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with  the  mean  value  of  229  Ib.  per  sq.  in.  afforded  by  the  curve  AB  of  stress 
distribution  across  the  principal  cross  section  (fig.  GB),  and  we  now  find  the 
maximum  stress  at  points  in  the  contour  where  a  tangent  line  makes  a  slight 
inclination  to  the  line  of  pull.  Its  value  is  640  Ib.  per  sq.  in.,  and  gives  the 
low  stress  concentration  factor  of  2.78. 

These  results  seem  to  show  that  the  form  of  the  notch  in  an  impact  test 
not  only  exercises  a  very  considerable  influence  on  the  results,  but  any  one 
type  of  notch  may  cause  that  influence  to  be  exerted  in  a  variable  manner 
according  to  the  material  tested. 

As  a  test  it  is  undoubtedly  complex  and  even  more  difficult  to  analyse  in  the 
beam  form,  and  its  simplification,  if  that  were  possible,  would  be  a  great 
advantage. 

It  might,  for  example,  be  worth  while  to  examine  the  possibilities  of  an 
impact  tension  test  specimen  formed  by  drilling  out  a  very  large  hole  centrally 
placed  in  a  flat  tension  member,  since  this  form  of  discontinuity  has  been 
shown  (4)  to  cause  an  approximately  linear  stress  distribution  across  the 
minimum  section  with  a  maximum  value  at  the  inner  contour,  and  a  minimum  at 
the  outer  contour,  which  is  nearly  but  not  quite  zero.  In  hard  materials, 
therefore,  where  the  elastic  stress  distribution  probably  remains  practically 
unchanged  up  to  fracture,  the  maximum  stress  is  very  approximately  twice  the 
mean  value,  while  in  ductile  materials  photo-elastic  indications  appear  to  show 
that  the  stress  distribution  becomes  fairly  uniform  at  fracture. 


COMPLEX  STRESS  DISTRIBUTIONS  IN  ENGINEERING  MATERIALS. 


299 


REFERENCES. 

(1)  'Stress   Distributions   in   Engineering  Materials.'     B.A.    Report,   1914. 

(2)  '  Photo-Elastic  Measurements  of  the  Stress  Distribution  in  Tension  Members 

used  in  the  Testing  of  Materials.'  By  Prof.  E.  G.  Coker,  F.R.S.  Min. 
Proc.  Inst.C.E.  Vol.  ccviii.  (1918-19).  Part  ii. 

(3)  *  Contact  Pressures  and  Stresses.'     Prof.  E.  G.  Coker,  K.  C.  Chakko,  and 

M.  S.  Ahmed,  Proc.  Inst.  Mech.  Engrs.     March  1921. 

(4)  *  Photo-Elastic  and  Strain  Measurements  of  the  Effects  of  Circular  Holes  on 

the  Distribution  of  Stress  in  Tension  Members.'  By  Prof.  E.  G.  Coker, 
K.  C.  Chakko,  and  Y.  Satake.  Trans.  The  Inst.  of  Engrs.  and  Ship- 
builders in  Scotland.  Vol.  Ixiii.  Part  i. 


II. 

The  Distribution  of  Stress  in  a  Flanged  Pipe. 

By  GILBERT  COOK,  D.Sc. 

The  problem  of  the  elastic  deformation  of  a  thin  cylindrical  shell  subjected 
to  internal  pressure  when  certain  boundary  conditions  restricting  the  displace- 
ment at  the  ends  are  to  be  satisfied  has  been  investigated  theoretically  by  Love,1 
and  as  applied  to  the  case  of  a  short  boiler,  in  more  detail  by  Nicolson.2 
The  particular  types  of  end  constraint  considered  were  those  in  which  the  ends 
of  the  tube  remain  circular  and  suffer  no  radial  displacement.  In  a  discussion 
of  certain  experiments  on  short  boilers3  it  has  been  shown  that  the  stress 
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FIG.  7. 


distribution  in  the  vicinity  of  the  end  constraints  is  of  a  somewhat  peculiar 
nature,  inasmuch  as  over  a  certain  region  the  radial  deformation  and  the 
corresponding  circumferential  stress  are  greater  than  that  which  would  obtain 
in  an  infinitely  long  tube,  but  little  experimental  evidence  of  this  fact  appears 
to  have  been  forthcoming.  A  third  type  of  constraint  which  possesses  consider- 
able practical  interest  is  that  produced  by  a  flange  or  collar  at  any  point  in 
the  length  of  the  tube.  At  this  point  the  inclination  of  the  generators  of  the 
cylinder  to  the  axis  is,  by  symmetry,  zero.  The  radial  displacement  is  not  zero, 
but  is  determinate,  and  the  constraint  is  of  a  type  which  may  be  more  accurately 
and  conveniently  reproduced  experimentally  than  those  referred  to  above. 

The  general  equation  of  equilibrium  for  any  type  of  end  constraint  may  be 
very  simply  obtained  as  follows  : 

Let  t  be  the  thickness  of  the  tube  (fig.  7),  assumed  small  in  comparison  with 
II,  the  mean  radius.  Let  z  be  the  radial  displacement,  q  the  circumferential 


1  Mathematical  Theory  of  Elasticity.     Third  Ed.     Art.  340. 

2  Nicolson,  '  The  Strength  of  Short  Flat-ended  Cylindrical  Boilers.' 
X.E.  Coast  Inst.  Engineers  and  Shipbuilders.     Vol.  vii.,  p.  205. 

3  Engineering.     Vol.  51   (1891),  p.  468. 
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stress,  and  /  the  longitudinal  stress  (tensile)  in  the  shell  at  any  point  distant 
x  measured  along  the  axis  from  some  fixed  point  0.  Let  P  be  the  internal 
pressure. 

Then  the  circumferential  strain  ^=£  —  ^=,    where  E  is  Young's  Modulus,  and 

i\     EJ    msii 

Poisson's  Ratio. 
m 

So  that 


Consider  the  equilibrium  of  a  longitudinal  strip  of  the  shell  subtending  an 
angle  8<J»  at  the  axis.  This  may  be  considered  to  act  as  a  beam  supported  in 
some  arbitrary  manner  at  the  points  of  constraint,  and  to  sustain  a  load  in  the 
plane  of  bending  which  will  be  the  resultant  in  that  plane  of  the  internal  pressure 
and  the  circumferential  stress.  The  internal  pressure  P  will  produce  a  uniform 
load  intensity  equal  to  PR  $<f>  per  unit  length,  while  the  component,  in  the  plane 
of  bending,  of  the  load  due  to  the  circumferential  stress  q  will  be  equal  to  qttify. 
The  component  of  the  longitudinal  stress  will  be  negligible.  The  resultant 
intensity  of  load  per  unit  length  in  the  plane  of  bending  is  therefore 


The  equation  connecting  load  and  displacement  in  a  beam  is 

-S- 

where  ,6  is  the  flexural  stiffness,  which  for  a  thin  plate  is  equal  to 


So  that  we  have         B 


Inserting  the  value  of  q  from  equation  (1),  we  obtain  the  equation  of  equilibrium  • 


-_-  (9\ 

m2-!  12  ^    R2  mR 

If  it  be  assumed  that  the  longitudinal  stress  is  that  due  to  the  internal  pressure 
when  the  ends  of  the  tube  are  closed, 

.PR 

f—2i 

and  the  equation  becomes     E    T_i  '  TO 


where 

It  reduces  to  the  form 

where 


m 


This  statement  involves  the  assumption  that  the  circumferential  stress  in 
the  tube  has  no  effect  on  the  flexural  stiffness.  A  more  rigorous,  but  less  simple, 
analysis  leading  to  equation  (2)  is  given  by  Nicolson,  loc.  cit.,  pp.  205-214. 
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and 

The  solution  of  this  equation  is 

•p'r>  2 
z=cos  nx  (A  cosh  rw;4-B  sinh  •/&:*;)  + sin  nx  (C  cosh  nx  +  D  sinh  nx)-\ — ^r-         .         (4) 

Consider  a  length  AB    (=1)   of  the  pipe  between  successive  flanges,  it    being 
assumed  for  the  present  purpose  that  the  flanges  are  spaced  equidistantly  along 
the  pipe,  and  let  the  origin  for  a;  be  at  A  (fig.  7). 
It  is  evident,  from  symmetry,  that 

—  =o  for  x=o  and  x—l 
ax 

and  also  z—z0  for  x=o  and  x=l 

where  ZQ  is  the  radial  displacement  at  the  flange,  as  yet  undetermined.  These 
four  conditions  enable  the  four  arbitrary  constants  A,  B,  C,  and  D  to  be 
expressed  in  terms  of  known  quantities  and  ZQ,  and  equation  (4)  becomes 

T>1T?  2        /"pl"R2  \       / 

2=-—  —I  ---  —  ZQ  j   (  cosh  nx  cos  nx— H  sinh  nx  cos  nx+H.  cosh  nx  sin  nx 

—  L  sinh  nx  sin  nx) (5) 

TJ     cosh  nl—  cos  nl 

where  il=-r-r — r^ — = ; 

sinh  wZ+sin  nl 

L_sinh  nl— sin  nl 
sinh  wZ+sin  nl 

When  the  distance  I  between  successive  flanges  is  large,  H  and  L  both  approxi- 
mate to  a  value  1  and  the  equation  giving  the  radial  displacement  becomes 


Z=Z\  —  (ZI—ZQ)  (cosh  nx—  smhnx)  (cos 
or  z=zi—  (zi—  ZQ)  e-"*(sin  n#-fcos  nx)        .....         (6) 

p/r>-2 

where  zi=—  =  —  ,  the  radial  displacement  in  a  uniform  tube  without  flanges. 
jfit 

It  remains  to  determine  zot  the  radial  displacement  at  the  flange.  For  this 
purpose  it  will  be  assumed  that  the  flange  and  the  part  of  the  tube  beneath 
it  form  a  homogeneous  ring  of  longitudinal  thickness  c  and  external  radius  Rj. 
The  forces  acting  on  the  ring  tending  to  increase  its  diameter  will  be  those 
due  to  (1)  the  internal  pressure  acting  on  its  inner  surface,  and  (2)  the  reactions 
between  the  ring  and  the  adjacent  tube.  The  latter  will  consist  of  a  shearing 
force  acting  radially  outwards  on  the  flange  over  the  whole  circumference  of 
the  tube.  Denoting  this  shearing  force  per  unit  length  measured  along  the 
circumference  of  the  pipe  by  F,  the  resultant  outward  force  per  unit  length  of 
the  circumference  will  be 

Pc+2F 

and  the  equivalent  pressure  per  unit  area  of  the  internal  surface  of  the  flange, 
tending  to  increase  the  diameter,  will  be 

PCJ-2F 
c 


It  will  be  assumed  that  the  longitudinal  stress  in  the  flange  is  conhned  to,  and 
uniformly  distributed  over,  a  ring  of  thickness  equal  to  the  thickness  of  the 
tube.  The  flange  may  then  be  regarded  as  a  compound  cylinder  under  internal 

pressure   PH  --    the  inner  portion  whose  thickness  is  t  having  a  longitudinal 

C 

stress  equal  to  that  in  the  tube  due  to  the  pressure  P.  Observing  that  at  the 
common  surface  the  radial  stress  and  radial  displacement  have  equal  values  in 
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the  two  portions,  the  radial  displacement  at  radius  R  can  be  shown,   by  the 
ordinary  theory  of  the  stress  in  a  thick  cylinder,  to  be  approximately  (i.e.  neg- 

lecting squares  of  ^) 

JCb 


p  2m-l  ,  2F  1  R(R-<)     Ri2(m+ 
'  * 


2m         c  J         Em     R^R-f  *)—  R2(R-<) 
2m-l  ,  2F 


Now  the  shearing  force  at  any  point  of  a  beam  is  equal  to 

fg* 

^s 

where  ft  is,    as  before,  the  flexural  stiffness.     At  any  point   of  a  longitudinal 
strip  (of  unit  width)  of  the  tube  wall  under  consideration,  it  will  be  given  by 

™    m2        *3     d»z 


Performing  the  required  differentiation  on  the  expression  for  z  given  by  equa- 
tion (6),  and  putting  jc  =  0,  the  expression  (8)  becomes 


o  vnr  —  1 

which,  from  a  consideration  of  the  signs,  will  be  the  shearing  force  at  the  end 
of  the  strip  acting  inwards  on  the  flange.  The  value  of  F,  acting  outwards, 
will  therefore  be 


The  substitution  of  this  value  of  F  in  equation  (7)  enables  zn  to  be  determined. 

1+2  _m*_     EM3,V 
2m 


P    2m—  1,2     m2       E 
'  ^' 


Thus 


20  = 


12  _nP_     E 
^ 


and  this  value  of  zn  substituted  in  equation  (6)  will  enable  the  radial  displace- 
ment at  any  point  to  be  determined,  and  from  this  the  circumferential  stress 
by  equation  (1). 

It  will  be  seen  from  equation  (6)  that  the  expression  for  the  radial  displace- 
ment contains  a  periodic  function  of  the  distance  from  the  flange,  and  that  the 
longitudinal  profile  consists  of  a  series  of  waves  whose  amplitude  rapidly 
diminishes  with  the  distance  from  the  flange. 

Thus,  if  the  equation  be  written  in  the  form 

z=zi—  (zi—  zoM*) 
<(>(x)  reaches  a  negative  maximum  when  x=*9    —  ,    ~.  .  .  . 

?&        7t         ¥1* 

For  steel,  taking  -  =0'295,  these  values  become 
m 


7-32  x/<R,  12'20-v/JR  .... 
and  <j>(a;)=--0432,  --0019,  -'00008.  .  .  . 

It  would  appear  therefore  that  the   constraint  imposed  by  a  flange  produces 
an  appreciable  increase  in  radial  displacement  over  a  region  extending  from 


i.e.  x  =1-83 
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In  the  limiting  case  in  which  displacement  at  the  flange  is  entirely  prevented 
the  increase  in  stress  may  thus  amount  to  4.3  per  cent,  of  the  normal  stress. 
The  elastic  deformation  at  the  flange  will  reduce  this  figure,  but  it  is  sufficiently 
great  to  be  measurable  experimentally. 

A  tube  of  mild  steel  was  accordingly  machined  accurately  to  an  internal 
diameter  of  1  in.  and  an  external  diameter  of  1.05  in.,  a  collar  being  left 
on  the  tube  having  an  external  diameter  of  1.5  in.  and  a  thickness  of 
.038  in.  Connections  were  provided  at  the  ends  for  attachment  to  a  pump 
and  gauge,  and  the  tube  was  free  to  take  up  the  longitudinal  stress  produced  by 
the  internal  pressure  applied.  The  increase  in  diameter  was  measured  by  an 
instrument  specially  designed  for  the  purpose,  and  similar  in  principle  to  that 
previously  used  in  another  connection5  by  the  author  and  Professor  Robertson. 
An  outline  sketch  of  the  instrument  is  shown  in  fig.  8.  A  light  frame  A  is 
held  on  the  cylinder  by  means  of  two  flat  springs  B  so  that  the  steel  point  P1  is 
set  at  one  end  of  the  diameter  whose  extension  is  to  be  measured.  In  contact 
with  the  cylinder  at  the  other  end  of  the  diameter  is  a  second  steel  point  P2 
carried  by  a  level  L,  which  is  free  to  turn  about  the  knife  edge  D,  and  is  held  in 
position  by  the  light  spring  Si  (consisting  of  a  rubber  thread).  Between 
V's  in  the  outer  end  of  the  lever  L  and  a  vertical  bar  E  supported  as  shown 
in  the  sketch  from  the  frame  A,  is  held  a  small  diamond-shaped  rocking  piece 
R  possessing  two  knife  edges,  one  edge  resting  in  the  V's,  and  the  other  against 


FIG.  8. 

the  flat  surface  of  the  vertical  bar  E.  The  rocking  piece  carries  on  the  con- 
tinuation of  its  axis  a  mirror  Mj  and  a  counterpoise.  The  angular  rotation  of 
this  mirror  is  proportional  to  the  change  in  diameter  between  the  steel  points 
PiP2,  and  is  measured  by  means  of  a  telescope  and  scale.  A  mirror  M2  (not 
shown  in  the  right-hand  view)  fixed  to  the  vertical  bar  E  serves  to  detect  any 
displacement  of  the  instrument  as  a  whole.  In  these  experiments  the  scale 
was  placed  at  a  distance  of  15  feet  from  the  mirrors,  and  by  using  a  telescope 
having  a  linear  magnification  of  about  25,  and  a  fine  cross  wire,  the  scale  could 
easily  be  read  to  0.1  mm.  Successive  readings .  of  the  scale  for  the  same 
pressure  range  differed  by  not  more  than  0.1  to  0.2  mm.  on  a  total  reading 
(for  parts  of  the  tube  remote  from  the  flange)  of  about  4*0  mm.  The  lever 
ratio  being  2.85,  and  the  rocking  piece  being  0.25  in.  bro^id,  diametrical  exten- 
sions could  be  measured  to  an  accuracy  of  the  order  of  10"6  in.  To  secure 
this  degree  of  accuracy  it  was  found  necessary  to  mount  the  apparatus  on  a 
massive  pedestal  supported  by  felt  pads  to  damp  out  small  tremors,  and  to 
maintain  a  constant  temperature.  It  was  not  possible  to  take  accurate 
measurements  within  half  an  hour  of  handling  the  instrument. 

The  pressure  gauge  used  was  a  standard  steel-tube  test  gauge  maVle  by 
the  Budenberg  Gauge  Co. 

In  the  test  the  increase  in  diameter  was  measured  at  various  distances  from 
the  flange  when  the  pressure  was  raised  from  100  to  900  Ib.  per  sq.  in., 
and  the  results  are  tabulated  below,  together  with  the  values  calculated  by  the 
theory  described  above.  For  the  latter  purpose  an  accurate  determination 

5  Engineering,  December  15,  1911. 
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of  the  elastic  constants  for  the  material  was  made.  Young's  Modulus  (for 
longitudinal  stress)  was  found  from  the  longitudinal  extension  in  a  cylindrical 
bar,  and  Poisson's  Ratio  from  a  measurement  of  the  lateral  contraction  of  the 
same  bar  by  the  instrument  used  on  the  tubes,  axial  loading  shackles  being 
employed  to  secure  uniformity  of  stress.  The  values  thus  found  were 

E=29'88X  106  Ib.  per  sq.  in. 

1=0-295. 
m 

The  results  are  also  plotted  in  fig.  9,  and  it  will  be  seen  that  the  experimental 
results  follow  closely  the  theoretical  curve.  The  small  discrepancy  which 
appears  at  points  remote  from  the  flange  is  entirely  accounted  for  by  the  fact 
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that  the  actual  extensions  were  measured  on  the  outside  of  a  tube  of  appreciable 
thickness,  whereas  the  theory  assumes  uniform  stress  across  the  walls.  Thus, 
the  radial  extension  at  the  outside  of  a  plain  tube  is  given  by 

E 

where  R'  and  R0  are  the  external  and  internal  radii  respectively,  giving,  with  the 
present  data,  a  value  23.95xlO~5,  almost  identical  with  that  observed  in  the  test. 
The  discrepancy  near  the  flange  is  probably  due  to  errors  in  the  assumptions 
involved  in  the  calculation  of  the  displacement  of  the  flange.  The  point  of 
maximum  displacement  is  situated  at  a  distance  of  about  .28  in.  (as  nearly  as 
can  be  measured)  from  the  flange,  the  calculated  distance  being  .276  in. 
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The  actual  increase  at  this  point  is  4.5  per  cent,  of  the  difference  between 
the  displacement  remote  from  the  flange  and  at  the  flange,  which  approaches 
very  closely  the  calculated  value  of  4.3  per  cent. 

Tests  have  also  been  carried  out  on  tubes  of  greater  thickness  and  with 
broader  flanges,  the  results  obtained  being  very  similar  to  those  described 
above,  but,  as  would  be  expected,  the  agreement  of  theory  and  experiment  is 
not  so  close. 

The  results  of  this  investigation  have  an  important  bearing  on  the  strength 
of  steel  pipes  reinforced  by  steel  bands,  sudh  as  are  occasionally  used  in  pipe 
lines  for  hydro-electric  installations.  Experimental  work  in  this  connection  is 
at  present  *in  progress. 


Radial  Displacement, 

Distance  from 

for  Pressure  800  Ib.  per  sq.  in. 

Flange  (ins.) 

Calculated 

Observed 

(ins.  x  10-5) 

(ins.  x  10-5) 

0 

10-60 



•025 

11-53 

12-05 

•05 

13-60 

14-20 

•075 

15-94 

16-30 

•100 

18-16 

18-80 

•125 

20-31 

20-80 

•150 

21-92 

21-95 

•175 

22-95 

22-80 

•200 

23-78 

23-30 

•225 

24-22 

23-80 

•250 

24-46 

23-8o 

•30 

24-52 

23-90 

•35 

24-32 

23-70 

•40 

24-10 

23-50 

•45 

24-06 

23-35 

•50 

24-02 

23-35 

in. 


On  Stresses  in  Multiply-connected  Plates. 

By  Professor  L.  N.  G.  FILON,  M.A.,  D.Sc.,  F.R.S.,  University  College,  London. 

1.  It  is  a  well-known  theorem,  due  to  Airy,  that  in  a  problem  of  plane  strain,  the 
axes  of  x  and  y  being  in  the  plane  of  the  strain,  the  stresses  can  be  expressed  in  terms 
of  a  single  function  E  by  the  equations 


—     82E     — 
##=-— 

Kti" 


*-»     S2E 


(1) 


us  denoting  the  stress,  parallel  to  *,  across  a  face  perpendicular  to  n. 

This  result  can  be  extended  to  the  case  where  the  stress  zz  normal  to  the  faces  of 
a  thick  plate  vanishes  throughout.  In  this  case,  P,  Q,  S  denoting  the  mean  stresses 
taken  acro^  the  faces  of  the  plate 


p=?!5    Q= 

8?/2 


S=- 


Sa% 


(2) 


In  either  case  the  function  E,  which  is  called  Airy's  stress  function,  satisfies  the 
equation 

V4E-0 (3) 
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If  A,  n-  are  the  elastic  constants  of  Lam6,  we  have,  in  the  first  case,  the  stresses 
»'x,  xy,  yy  given  in  terras  of  the  displacements  u,  v  by  the  equations 


Sy 


, 

" 


(4) 


In  the  second  case  P,  Q,  S  are  given  in  terms  of  the  mean  displacements  U,  V  by 
precisely  similar  equations,  A  being  replaced  by  A',  where  A'  =  2A/x/(A  +  2u).  The 
two  problems  are  thus  analytically  the  same,  and  in  what  follows  we  shall  confine 
ourselves  to  the  second  one,  which  is  known  as  that  of  generalised  plane  stress. 

2.  The  conditions  at  a  boundary,  where  n  is  the  direction  of  the  outwards  normal 
whose  direction  cosines  are  Z,  m,  are  given  by 


X  and  Y  being  the  mean  stresses  across  the  boundary  parallel  to  the  axes. 
Equations  (5)  can  be  written  in  the  form 

S/SEv 

**(to)    .....      (G) 

ds  being  an  element  of  the  boundary. 

Since  the  differential  equation  (3)  satisfied  by  E  and  the  boundary  conditions  (6) 
for  E  do  not  -involve  the  elastic  constants,  it  would  at  first  sight  appear  that  the 
solution  for  E,  for  any  given  set  of  boundary  stresses,  cannot  involve  the  elastic 
constants  and  that  the  stress-distribution  in  a  plate,  under  given  applied  boundary 
stresses  in  the  plane  of  the  plate,  is  independent  of  tha  material  of  the  plate, 
provided  only  it  be  elastic. 

This  proposition,  if  true,  is  of  great  practical  importance,  because  it  shows  that  if 
we  can  find  the  stress  distribution  in  a  plate  by  any  method  whatever,  for  example 
by  the  exploration  of  a  plate  of  xylonite  by  means  of  polarised  light,  the  results,  so 
far  as  the  stresses  are  concerned,  can  be  transferred  immediately  to  a  plate  of  any 
other  material,  such  as  steel,  provided  the  applied  stress-system  is  the  same. 

It  was  first  shown  by  J.  H.  Michell  (1)  that  the  theorem  in  question  does  not 
hold  when  the  area  of  the  plate  is  multiply  connected,  that  is,  when  it  contains  one 
or  more  Tioles,  so  that  closed  circuits  can  be  drawn  inside  the  plate,  which  cannot  be 
made  to  shrink  into  a  point  by  continuous  deformation  without  passing  outside  the 
material  of  the  plate. 

The  problem  was  discussed  by  A.  Timpe  (2)  in  connection  with  a  ring  boundary 
bounded  by  concentric  circles,  and  Timpe  was  the  first  to  point  out  that  the  failure 
of  the  theorem  was  due  to  the  existence  of  types  of  strain  in  such  a  multiply- 
connected  plate,  created  by  removing  a  thin  strip  or  wedge  of  material,  thus  cutting 
a  thin  channel  or  gap,  with  straight  edges,  from  the  hole  to  the  outer  boundary  and 
then  closing  up  the  gap  and  cementing  its  faces  together.  We  then  obtain  a  stressed 
ring,  although  no  forces  are  applied  to  the  circular  boundaries. 

More  recently  Prof.  Vito  Volterra  (3)  has  given  a  general  theory  of  such  types  of 
strain,  to  which  he  has  given  the  name  of  distorsions,  and  for  which  Prof.  A.  E.  H. 
Love  (4)  has  proposed  the  English  term  dislocations.  Prof.  Volterra's  account  is  not 
restricted  to  two  dimensions,  and  he  has  incorporated  in  it  various  scattered  results 
obtained  by  Weingarten,  Cesaro,  and  Tedone.  He  does  not,  however,  discuss  the 
question  whether  the  solution  is  independent  of  the  elastic  constants,  as  this  question 
does  not  really  arise  in  the  three-dimensional  case.  So  far  as  I  have  been  able  to 
find  out,  MichelPs  has  been  the  latest  statement  on  this  particular  subject. 

3.  This  result  of  Michell's  necessarily  throws  doubt  upon  the  general  applicability 
of  investigations  with  xylonite  models  to  ordinary  engineering  structures,  since  the 
ratio  of  the  elastic  constants  is  different  in  xylonite  and  in  the  usual  engineering 
materials.  M.  Mesnager  (5),  in  answering  a  similar  objection  to  his  use  of  glass  models, 
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appears,  however,  unaware  that  the  unique  determination  of  the  stress  function, 
independently  of  the  elastic  constants,  is  subject  to  exceptions. 

The  following  sections  give  a  short  account  of  the  theory  of  such  dislocations,  in 
a  plane,  on  lines  somewhat  more  direct  than  those  followed  by  Michel  1  and  Volterra. 
One  important  result  obtained  by  this  method  of  attack  is  to  show  that  the  correc- 
tions to  the  stresses  under  a  given  system  of  loads,  when  we  pass  from  one  material 
to  another,  can  be  obtained  directly  from  exploration  of  the  stresses  in  one  kind  of 
material  only,  by  a  subsidiary  experiment  with  a  dislocated  plate,  provided  Poisson's 
ratio  is  known  for  both  materials;  this  correction,  for  any  given  shape  of  plate  and 
applied  force  system,  can  thus  be  determined  experimentally  by  optical  means,  even 
when  a  mathematical  solution  is  not  available. 

It  seemed  also  of  interest  to  determine,  for  a  shipe  of  plate  which  allows  of 
mathematical  solution,  an  upper  limit  to  the  divergences  to  be  expected  in  the 
stresses,  when  we  pass  from  one  material  to  another.  This  has  been  done  here  for  a 
ring-shaped  plate ;  the  final  results  are  simple  in  form  and  give  indications  which 
should  be  valuable  to  the  engineer.  In  particular  they  show  that  if  the  width  of  the 
ring  be  moderate  in  comparison  with  its  radius,  the  corrections,  in  practice,  will  be 
very  small,  at  any  rate  at  all  points  where  the  stresses  are  of  any  magnitude. 

4.  In  what  follows  we  shall  have  to  consider  functions  of  which  the  derived 
functions  of  a  certain  order  are  one-valued,  whilst  the  functions  themselves  and 
their  derived  functions  of  lower  order  are  many-valued. 

Let  <J>  be  a  function  of  x,  y,  and  let  its  differential  coefficients  be  denned  in  a  non- 
ambiguous  manner  all  over  the  area  contained  by  a  closed  circuit  APBQA  (fig.  10), 


FIG.  10. 

* 

which  can  be  made  to  shrink  to  a  point  without  passing  outside  the  material.     Such 
a  circuit  is  said  to  be  reducible. 

Then  if  ($>),  be  the  value  of  <t>  selected  at  A  at  the  beginning  of  the  circuit 
and  (».,  the  value  reached  at  A  after  describing  the  circuit 


—  (</>),  =       —  -- 
JA  8* 


-  ds  taken  round  the  circuit. 


But  since   _,    —  .  have  been  defined  in  a  non-ambiguous  manner,  then,  in  the  figure 


(bearing  in  mind  the  sense  of  describing  the  contour) 


Q  5//  V5.A 

where  dxt  dy  are  now  positive  throughout. 
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Thus 


taken  over  the  whole  area  of  the  circuit. 
Similarly 


fA  5<f>   ,  f  5    /5d>\ 

-i:  dy  =  +    —  (  J:  ) 

JA   fy  J8#  Vfy/ 


Thus 


=  o, 

if  inside  the  circuit  —-J-  exists  and  is  denned  in  an  unambiguous  manner. 


Now  it  is  clear  that  if 


I 


IAPB     os  JBQA  °* 

then,  reversing  the  sign  of  the  second  integration 


JAPB   8*          JAQB  5* 


OB  —  <J»A)  along  APB  =  (4>B  -  <£A)  along  AQB, 

that  is  we  obtain  the  same  value  of  <£  at  B  if  we  build  it  up  by  integration  along 
any  two  reconcilable  paths,  starting  from  a  definite  selected  value  at  a  point  A. 

The  value  of  (/>  thus  obtained  at  any  point  which  can  be  reached  from  A  by  a 
continuous  path  is  thus  denned  in  an  unambiguous  manner,  provided  the  connectivity 
of  the  path,  i.e.,  the  manner  in  which  it  winds  amongst  the  holes  in  the  material, 
is  specified. 

Now,  if  A  is  a  point  of  a  reducible  circuit,  all  paths  AR  from  A  to  a  point  inside 
or  on  the  circuit,  are  clearly  reconcilable.  We  thus  get  a  function  <?>  unambiguously 

defined  in  the  same  manner  that  ~  and  ~  were  unambiguously  defined. 

ox  oy 

It  follows  Uiat  we  can  in  this  way  construct  the  function    <f>   from  its   first 

differential  coefficients,  provided  —  ,    —  and  —  2-  exist  and  can  be  unambiguously 

So?      8 


defined  in  the  region  considered. 

If  we  continue  our  function  <j>,  obtained  in  this  way,  round  an  irreducible  circuit 
(enclosing  say  a  hole  in  the  plate),  then  on  reaching  again  the  same  point  A  we 
get,  in  general,  a  value  4>2  different  from  </>,.  We  get,  however,  the  same  value  <f>z 
for  any  other  irreducible  circuit  reconcilable  with  the  first  one,  so  that  $>2  is  a  definite 
function  of  the  co-ordinates  x,  y  of  A  determined  by  the  connectivity  of  the 
irreducible  circuit. 

We  shall  call  <f>2  —  <j>,  the  cyclic  function  of  </>  for  this  particular  connectivity  and 
denote  it  by  Cy(<f>)  ;  it  will  usually  be  convenient  to  omit  the  indication  of  the 
circuit  typical  of  the  connectivity  :  this  may  be  specified  by  a  suffix  denoting  the 
circuit  in  question,  where  required. 

It  will  follow  directly  from  the  definition  that 


<*(£)  •" 

so  that  Cy  and  —  (  or  —  )  are  commutative  symbols  of  operation. 
5a?  \      oyj 

5.  In  what  follows  the  stresses  and  strains  will  be  restricted  to  be  essentially 
one-valued,  that  is,  no  restriction  has  to  be  placed  on  the  functions  which  represent 
them  in  order  to  make  them  one-valued  over  a  given  region  and  this  quite 
irrespectively  of  multiple  connection  of  the  space  considered. 
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It  is  easily  shown  that  the  displacements  are  obtained  in  terms  of  the  stress 
function  as  follows 


(7) 

where 

V2*=0 .     (8) 

J5|L=V2E (9) 

and  ff  =  \\'(\'  +/*),  so  that  (1  —  <r)  =  (1  +  *?)-1,  where  TJ  is  Poisson's  Ratio. 

We  note  first  of  all  that  — -,    — ,  all  exist  and  are    one- valued,  and  that 

5a?2      5y2      5a?5?/ 

the    same    holds  good  of    the    third  differential    coefficients.     From  this,   by  the 
theorem  ot  §  4,  it  follows  that   _,    —  are  acyclic  for  a  reducible  circuit  and  have 

a  definite  cyclic  function  for  an  irreducible  circuit.     By  repeated  application  of  the 
theorem  the  same  result  holds  good  for  E. 

*"ni          XT?  • 

The  cyclic  functions  for   -    and  —  are  very  readily  found  from  equations  (6). 
Sx  oy 

Taking  a  circuit  enclosing  one  hole  and  one  only,  we  have,  integrating  equations 
(6)  along  this  circuit 

(t"W\  A 

—  )=-    Yrf*=  —  Total    force    resultant  on  the  circuit    parallel    to   y*=  Force 
oa?/         J 

resultant  on  boundary  of  enclosed  hole  parallel  to  y  =  Y0,  say. 

(XTT\ 
—  1  =  -~X0,.  where  X0  =  Force  resultant  parallel  to  x  acting  on  the 

boundary  of  enclosed  hole.     This  result  is  due  to  Michell. 
In  a  similar  manner  it  can  be  shown  that 

.         .         .         .         .     (10) 


where  M0  is  the  couple  resultant  at  the  origin  of  the  forces  acting  on  the  inner 
boundary  of  the  enclosed  hole.  The  value  of  Cy(E)  at  any  point  is  thus  the  moment 
about  that  point  of  the  forces  applied  to  the  boundary  of  the  hole  about  which  the 
typical  circuit  is  taken. 

We  have  next  to  consider  the  cyclic  functions  of  _i  and  2?.     First  of  all  we  note 

Sy         Sat 

that  since  V  2E  is  essentially  one-valued,  so  is   —  I-  and  therefore  also    5L  and 


.  And  since  V2ij/  =  0  .*.—  r  +  —  ^-  =  0,  so  that,  J2-  being  essentially  one-  valued 
3 


so  is  —t  and,  similarly,  so  is   -T. 
ox  oy? 

Thus  the  third  differential  coefficients  of  ij>  are  essentially  one-valued,  and 
accordingly  so  will  the  fourth  differential  coefficients  be.  Hence  by  the  theorem  of 
§4  the  second  differential  coefficients  have  cyclic  functions  (which  involves  being 
acyclic  for  reducible  circuits  and  being  unambiguously  defined  for  reconcilable 
paths).  The  first  differential  coefficients  and  i|>  itself  will  therefore  have  cyclic 
functions  for  any  typical  irreducible  circuit. 

Applying  Cy  to  equations  (8)  &nd  (9)  and  remembering  that  V2E  is  acyclic,  we 
have 

0        .  ....     (11) 

.         ...         -         .     (12) 
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These  two  equations  require  that 


where  A,  B,  C,  D  are  constants.     Thus 

A.<?  +  B 


"(2)- 


6.  Referring  now  to  equations  (7)  and  applying  the  operation  Cy,  we  have 


-0         --  (     } 

2,uCy  (  V)  =     X0  +  (  I  -  <r)[Atf  +  B]  j 

Thus  Cy(U),  Cy(V)  are  equivalent  to  a  rigid  body  displacement  consisting  of  a 
(small)  translation  of  components  (  -  Y0  +  1  —  <rC)/2M  parallel  to  x  and  (X,,  +  (  1  -  <r)  B)/2,u 
parallel  to  y  and  a  small  rotation  (1  -o-)A/2/*  about  the  origin. 

This  is  Weingarten's  theorem  (6).  So  soon  as  it  is  given  that  U  and  V  have 
definite  cyclic  functions  the  result  is  obvious.  For  since  U,,  V,  and  U2V,  are  both 
solutions  of  the  equations  of  elasticity,  U2  —  U,,  V2  —  V,,  that  is  Cy(U)  and  Cy(V)  are 
also  solutions.  But  the  stresses  obtained  from  these  must  be  zero,  since  the  stresses 
are  essentially  acyclic.  Thus  Cy(U)  and  Cy(V),  by  a  well-known  result,  are 
necessarily  a  rigid  body  displacement. 

We  can  now  introduce  Timpe  and  Volterra's  interpretation  of  Cy(U)  and  Cy(V). 
Reduce  the  multiplicity  of  connection  by  drawing  a  barrier  CDE  (fig.  11)  from  the 


FIG.  11. 

boundary  of  any  inner  hole  H  to  the  outermost  boundary  of  the  plate. 

This  cut  need  not  be  straight. 

Now  give  the  cut  a  very  small  rigid-body  displacement  so  that  it  takes  up  the 
position  C'D'E'.  Cut  the  plate  so  that  CDE  is  one  boundary  of  the  gap  and  C'D'E' 
the  other.  If  necessary,  wedges  of  material  may  have  to  be  cemented  on  to  effect 
this.  Now  cement  the  two  boundaries  together,  so  that  the  plate  once  more  forms  a 
whole.  Then,  in  consequence  of  this  dislocation,  strain  is  introduced,  corresponding 
to  displacements  having  for  cyclic  function  with  respect  to  the  hole  H  the  rigid  body 
displacement  which  transformed  CDE  into  C'D'E'.  With  regard  to  any  other  hole 
the  displacements  are  still  acyclic. 

In  general  there  will  be  a  cyclic  function  for  each  hole,  each  such  function 
involving  three  constants,  so  that  if  there  are  n  holes,  so  that  the  multiplicity  of 
connection  of  the  plate  is  n  +  1,  these  are  Sn  cyclic  constants.  By  cutting  n 
channels  of  appropriate  shape  and  recementing  we  obtain  a  plate  under  internal 
strain  and  stress,  in  which  the  displacements  are  affected  by  the  correct  cyclic 
functions.  If  now  external  forces  are  applied  to  the  boundary  of  the  cemented 
plate,  these  introduce  a  strain  which  involves  only  acyclic  displacements,  and 
so  leaves  the  cyclic  functions  unaltered. 

We  thus  obtain  a  representation  of  the  most  general  type  of  strain  of  such 
a  plate  as  a  combination  of  an  external  force  system  with  a  set  of  dislocations. 
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7.  Let  us  now  consider  in  what  way  the  elastic  constants  are  introduced  into 
the  solution. 

Imagine  a  plate  of  the  same  size  and  shape  as  the  one  considered,  but  made 
of  some  ideal  material  whose  elastic  constants  have  fixed  numerical  values.  Then, 
if  the  given  forces  are  applied  to  this  ideal  plate,  there  will  necessarily  on  physical 
grounds  be  a  solution  which  will  lead  to  acyclic  displacements  U0,  V0.  Let  E0  be 
the  value  of  E  and  ^0  the  value  of  ^  corresponding  to  this  solution. 
Then,  as  before, 


u0,  ff0  referring  to  the  'idenl  '  material,  so  that  these  are  fixed  numbers. 

Now,  if  the  plate  is  simply  connected,  every  circuit  drawn  on  it  is  reducible, 
hence  E0  and  ^0  are  necessarily  acyclic. 

If  we  now  consider  the  same  plate,  made  of  any  actual  material,  and  write  down 
the  displacements 


S// 

(Iti) 


These  displacements  will  clearly  lead  to  the  same  stresses  as  before  —  since  the 
stresses  depend  on  E  only,  and  E  is  here  E0.  Also  since  E0  and  i^0  are  both  acyclic 
so  are  U  and  V.  We  have  therefore  the  solution  required. 

Thus,  for  any  material,  E  =  E0,  and  so  is  independent  of  the  elastic  constants. 

We  have  now  to  inquire  how  this  is  modified  when  the  plate  is  multiply 
connected,  there  being  no  dislocations.  This  requires  that  U  and  V  shall  be  acyclic. 

First  consider  the  case  where  the  forces  applied  to  each  boundary  separately 

7\  R*        XT?1 

reduce  to  a  couple.  This  requires  (§5)  that  __,  _  are  acyclic.  Then,  consider- 
ing the  plate  of  ideal  material 


(1-0  Cy  =  2/ioCy  (V0)  +  Cy         >     =  0. 


Hence  -J^»  and  -^°   are  again  acyclic,  equations  (16),  therefore,  lead  to  acyclic 

values  of  U  and  V  for  the  actual  plate  whatever  <r  may  be,  and  we  can  take  E  =  E0 
and  thus  again  the  stresses  are  independent  of  the  elastic  constants. 

If,  however,  the  forces  applied  to  each  boundary  do  not  reduce  separately  to 

SF  SR! 

a  couple,  then  —  _!£  and  —  5  are  cyclic,  and  in  order  to  make  the  displacements 

U0,  V0  acyclic,  it  is  necessary  that 


Hence,  for  the  actual  material,  when 
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Thus  we  cannot,  in  this  case,  take  E  =  E0  and  ty  =  tj/0  and  the  function  E  (and  therefore 
also  the  stresses)  must  involve  the  elastic  constants. 

8.  Let  us  now  return  to  our  plate  of  ideal  material  and  introduce  in  it  a  simple 
dislocation  with  respect  to  any  given  hole  corresponding  to  a  unit  relative  translation 
of  the  faces  of  the  cut  parallel  to  x.  (The  unit,  of  course,  must  in  this  case  be 
a  small  length,  so  that  the  squares  of  the  strains  introduced  are  negligible.)  No 
forces  are  to  be  applied  to  any  boundary. 

Corresponding  to  this  dislocation  of  the  ideal  plate,  there  will  be  an  internal 
strain  and  stress  denned  by  the  stress  function  E,,  and  let  <J/j  be  the  corresponding 
value  of  ^.  Owing  to  the  absence  of  external  forces,  E,  and  its  differential 
coefficients  are  acyclic  throughout.  The  corresponding  displacements  U,,  V,  are 
cyclic  for  a  circuit  enclosing  the  hole  considered.  Also 

Cy(U,)  =  l,        Cy(V,)  =  0. 
Hence,  applying  Cy  to  (15). 

--  .        .        .'(,8) 


Similarly  if  E,,  t|/2  belong  to  a  simple  dislocation  with  respect  to  the  same  hole 
corresponding  to  a  unit  relative  translation  of  the  faces  of  the  cut  parallel  to  y 


In  like  manner,  unit  translational  dislocations  with  reference  to  thfe  other  holes 
will  lead  to  stress  functions  ES,E4,  .  .  .  E^.^E-j,,  and  corresponding  ^  functions 
Wv  •  •  •  *«—i.*aii. 

Build  up  a  solution 

E  =  E0  +  a,E,  +  c^Eg  +  e^E,  +  o4E4  +    .  .  . 


where  a,  ...  aan  are  constants  to  be  determined. 

If  this  solution  is  to  lead  to  acyclic  values  for  the  displacements  U,  V,  we  must 
have,  for  every  irreducible  circuit 

' 


If  we  take  Cy  to  refer  to  an  irreducible  circuit  enclosing  the  hole  to  which  solu- 
tions 1  and  2  refer,  but  none  of  the  others,  then  Cv(E,)  =  Cy(K2)  =  Cy(Es)  =  Cy('E4) 
=  .  .  .  .  =Cy(B2»  =  0,  and  also,  since  the  displacements  U,,  V8,  >ay,  are  acyclic  except 
for  a  circuit  enclosing  that  hole  to  which  solutions  3  and  4.  refer, 

•••  =*• 

Equations  (20)  then  become,  on  substitution, 


whence,  using  .(18)  and  (19), 

!_  Cvf^°\-  Cv(^°\-  2/tt°ai  } 

1  ~         y\     it*.     I     ~          ^\     X..     I          1    _yr-        I 

.         .         .     (21) 


Since    the  cyclic  functions  on   the  left-hand   side  are  known  to  reduce  to  mere 
constants,  equations  (21)  determine  a,  and  o2. 

Similarly,  taking  irreducible  circuits  about  the  other  holes,  we  obtain  ets,  o4,  etc., 
and  we  have  finally  an  acyclic  solution  which  holds  good  for  the  actual  plate,  but 


COMPLEX  STRESS  DISTRIBUTIONS  IN  ENGINEERING  MATERIALS.      313 

this  solution  is  at  once  seen  to  involve  the  ratio  of   the  elastic  constants,  since  a 
enters  into  the  values  of  a,,  c^  determined  by  equations  (21). 

These  values  can  be  simply  expressed  in  terms  of  the  force  resultants  X0,  Y0 
applied  to  the  corresponding  hole.     For  using  (17)  we  have 


9.  These  results  are  going  to  enable  us,  not  only  to  obtain  an  upper  limit  to  the 
divergences  in  the  stresses  in  two  circular  rings  of  different  materials  under  the  same 
applied  forces,  owing  to  the  divergence  in  the  ratio  of  the  elastic  constants,  but, 
what  is  even  more  important  to  correct  definitely  results  obtained  on  a  multiply- 
connected  xylonite  model  and  make  them  applicable  to  a  material  like  steel,  even 
when  the  mathematical  solution  cannot  be  attained. 

Let  us  take  this  latter  application  first. 

Suppose  that  our  '  ideal  '  material  is  now  taken  to  be  xylonite,  and  that  in  a  plate 
of  xylonite  the  stresses  P,  Q,  S  have  been  measured  by  the  optical  and  transverse 
contraction  methods  described  by  Professor  Coker  [(7)  (8)  and  (9)],  so  that  we  may 
consider  P0,  Q0,  S0  as  known  (the  suffixes  having  the  same  meaning  as  before). 
Suppose  further  that  by  a  separate  experiment  with  a  xylonite  model  suitably  cut 
and  re-cemented,  the  stresses  P,,  Qlf  S,  ;  P2,  Q2,  S2  ;  etc.,  corresponding  to  specified 
unit  dislocations,  have  been  similarly  explored.  As  the  only  dislocations  considered 
are  translational,  the  experiment  is  very  easily  carried  out  in  the  majority  of  cases 
by  cutting  a  straight  channel  perpendicular  to  the  direction  of  translation  required 
and  bringing  the  edges  of  the  cut  together. 


FIG.  12. 

Now  let  the 'actual'  material  be  any  ordinary  engineering  material,  e.g.  steel, 
then  the  stresses  in  this  material,  due  to  the  same  force  system  already  applied  to 
the  celluloid,  are  given  by 


S  =  S0  +  0,8,4-0282+     .  .  . 
where  a,,  a2,  as,  .  .  .  are  given  by  equations  (22)  and  similar  equations. 

Now  the  total  force  resultants  (per  unit  thickness)  applied  to  the  various 
boundaries  are,  in  most  problems  to  which  this  method  can  be  applied,  directly 
known,  and  the  values  of  TJ,  T/O,  /tt0  for  the  materials  concerned  can  be  found  without 
difficulty  by  direct  experiment  when  they  are  not  alreadj7  known.  Thus  the  o's  are 
calculated  without  difficulty  and  the  stresses  in  the  steel  plate  at  all  points  can  be 
deduced.  Accordingly  the  exploration  of  a  xylonite  model  can,  even  in  this  case,  be 
made  to  give  complete  information  about  the  stresses  in  a  plate  of  any  material, 
even  though  the  mathematical  solution  is  beyond  our  present  powers  of  analysis. 

10.  We  will  now  apply  this  method  to  find  the  magnitude  of  the  corrections 
which  have  to  be  applied  to  the  stress-system  observed  in  a  xylonite  ring,  the  radii  of 
the  inner  and  outer  boundaries  being  a  and  b.  In  order  to  do  so  it  is  necessary  first 
of  all  to  work  out  the  stresses  for  the  two  elementary  translational  dislocations 
represented  in  fig.  12. 

1921  z 
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We  can  further  simplify  by  taking  the  axis  of  y  parallel  to  the  force  resultant 
on  the  inner  boundary.  We  have  then  X0  =  0,  and  therefore  a2  =  0,  so  that  only  the 
dislocation  due  to  translation  parallel  to  x  need  be  considered  (fig.  12,  right-hand). 

The  stress  function  for  this  dislocation  has  been  given  by  Timpe  (loo.  cit.,  p.  31). 
It  can  also  be  deduced  from  Volterra's  general  results.  In  this  case 


E,  -  K  sin  e  - 


log  r] 


leading  to  the  following  stresses,  r,  9  being  polar  co-ordinates 
^       K  sin  6 


We  find  also 

f,  =  2K{ia>3-^2+(«2  +  &2)  (*  log  r-y«-»)}, 
whence 


and  therefore 


so  that  equation  (18)  gives 

2/t0/(l  -  <r0)  =  -  4irK  (a2  +  J2), 
or 


and,  using  (22) 

«K=-  (*-"o)Yo 
4ir(a2  +  &2)' 

The  corrective  terms  to  the  stresses,  which  we  may  call  for  shortness  Arr,  Ar0,  A00> 
are  then 


(24) 


Note  carefully  that  these  corrective  terms  depend  only  upon  the  total  force 
resultants  applied  to  the  hole  and  not  in  any  may  upon  the  manner  in  which  these 
force  resultants  are  distributed;  and  this  remark  holds  good,  not  merely  for  the 
circular  ring,  but  for  the  general  case  of  a  plate  of  any  shape. 

11.  We  have  now  to  consider  where  these  corrective  stresses  are  greatest. 
It  is  clear  from  equations  (24)  that  for  each  stress  the  maximum  occurs  for  the 
same  value  of  r  on  each  radius  vector.     These  maxima  will  themselves  reach  their 

greatest  values:  (1)  in  the  case' of  Arr  and  A00  when  0=  ±7r/2,  that  is  along  the 

diameter  parallel  to  the  resultant  applied  force  ;  (2)  in  the  case  of  Ar0  when  0  =  0  or 
IT,  that  is  along  the  diameter  perpendicular  to  the  resultant  applied  force. 

Consider  now  the  maximum  corrective  stress  A00.  It  will  occur  when  <j>  =  «2J2/r3  —  3r 
+  (a*  +  J2)/r  has  its  greatest  numerical  value.  Now  d<f>jdr  -  -  SaWfr*  -  3  -  (a2  +  62)/r2 
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and  is  essentially  negative.  Thus  the  greatest  algebraic  value  of  (ft  occurs  at  the  inner 
radius  r  —  a  when  <f>  =  2(7>2  — «2)/«  and  the  least  algebraic  value  at  the  outer  radius 
r**b  where  <f>  =  (a2  —  &2)/J.  We  thus  see  that  </>  changes  sign  between  the  two  and 
that  the  greatest  numerical  value  is  at  r  =  a,  leading  to  : 


from  which  it  appears  at  once  that  this  corrective  stress  becomes  small  if  b[a  is  near 
unity,  that  is,  if  the  thickness  of  the  ring  is  small  in  comparison  with  its  radius. 
If  we  take  rj0  =  0.2,  T?  =  0.3  and  b\a  as  large  as  1.25,  we  get  A00max.  =  (0.022)  YJZira,  i.e. 
about  2%  of  a  stress  equal  to  the  total  load  distributed  uniformly  over  the  internal 
circumference.  This  last  stress  will  in  general  itself  be  small  compared  with  the 
bigger  existing  stresses,  so  that,  in  this  case,  the  correction,  at  any  rate  near  critical 
or  dangerous  points,  will  be  entirely  negligible. 

Coming  now  to  Arr  and  ArQ,  we  have  to  make  <j>  =  (ft2  —  r*)(r2  —  #2)/r3  a  maximum. 
Here  we  have  clearly  no  change  of  sign  inside  the  ring  and  a  positive  maximum 
occurs  between  r  =  a  and  r  =  b. 

The  value  of  r  corresponding  to  this  maximum  is  readily  found  from  : 

i/(02~+T 


This,  however,  leads  to  a  somewhat  awkward  algebraic  expression  for  the 
maximum  stress  corrections,  involving  lengthy  radicals,  and  since  we  are  really 
chiefly  concerned  with  finding  an  upper  limit  for  the  stress  corrections,  we  notice 
that  we  necessarily  increase  <p,  if  in  the  denominator  we  replace  r3  by  its  least 
possible  value  a3.  The  greatest  value  of  (b2  —  r2)  (r2  —  a2)  is  then  well  known  to 
occur  when  r2  =  (a2  +  J2)/2,  so  that  the  maximum  value  of  <J>  is  certainly  less  than 

i(y-g*)2t  leading  to 


Since  this  contains  .the  square  of  W—a2  it  leads  in  general  to  an  even  smaller 
correction  than  for  the  stress  66.  Taking  the  same  values  of  TJ,  T?O  and  bja  as  before, 
\ve  find 


<  (0'00154)Y0/2*YJ 
which  will  usually  be  entirely  negligible. 

The  investigation  of  the  case  of  the  circular  ring  therefore  indicates  that  the 
correction  due  to  variation  in  the  ratio  of  the  elastic  constants  is  usually  very  small. 
The  more  general  theorems  previously  obtained  show  how,  even  when  this  correction 
cannot  be  neglected,  or  computed  mathematically,  it  can  nevertheless  be  allowed 
for  if  a  suitable  experimental  exploration  is  undertaken.  This  justifies  the  use  of 
xylonite  models  for  the  exploration  of  stress,  even  when  such  models  are  multiply 
connected. 
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IV. 

Stress  Concentrations  in  Theory  and  Practice. 

By 'A.  A.   GRIFFITH,  D. Eng.,  of  the  Royal  Aircraft  Establishment. 

1.  Introduction. 

It  has  long  been  realised,  as  a  result  both  of  experience  and  theoretical  con- 
siderations, that  sharp  re-entrant  corners,  rapid  changes  of  section  and  the  like, 
are,  in  general,  undesirable  in  members  subjected  to  considerable  stresses,  by 
reason  of  the  local  weakening  which  they  cause.  At  the  same  time,  very  few 
attempts  have  been  made  to  take  account  of  these  stress  concentrations  in 
design,  partly  because  of  the  inherent  difficulty  of  the  subject,  but  also  very 
largely  because  calculations  based  on  the  mathematical  theory  of  elasticity  have 
in  some  cases  appeared  to  overrate  greatly  the  weakening  effect  which  is 
produced. 

During  recent  years,  however,  the  greater  stringency  of  the  working  con- 
ditions imposed  on  machines  and  structures  has  so  largely  increased  the  number 
of  failures  arising  from  stress  concentration  as  to  render  imperative  a  better 
understanding  of  the  mechanism  of  such  failures,  and  particularly  of  the  causes 
which  limit  the  validity  of  the  mathematical  methods  of  attack. 

In  the  present  paper  it  is  proposed  to  review  briefly  the  existing  state  of 
knowledge  of  this  subject,  and  to  discuss  some  of  the  problems  which  still 
await  solution. 

2.  Methods  of  Estimating  Stresses. 

The  following  methods  are  at  present  available  for  estimating  the  magnitude 
of  stress  concentrations  : — 

I.   Direct  Mathematical  Calculation.. 
II.  The  Photo  elastic  Method 

III.  The  Soap-film  Method. 

IV.  The  Thermal  Method. 
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So  far  as  stress  concentrations  are  concerned,  the  most  important  numerical 
results  which  have  been  obtained  by  method  I.  are  those  relating  to  the  stresses 
in  a  plate  containing  an  elliptic  hole  of  any  eccentricity  (1)  (and,  approximately, 
those  in  a  plate  from  whose  edge  springs  a  small  semi-elliptic  crack  or  groove) ; 
and  in  some  cases  the  stresses  in  bent  and  twisted  cylindrical  shafts,  the 
boundaries  of  whose  cross-sections  include  re-entrant  portions.  These  results 
may  further  be  extended  to  find  approximately  the  stress  concentrations  due  to 
small  semi-elliptic  holes  or  cracks,  and  also  scratches  or  grooves,  formed  in  the 
surface  of  members  which  are  not  cylindrical,  but  in  which  the  general  stress 
distribution  has  been  found  by  other  methods. 

In  a  further  application,  the  general  two-dimensional  solution  in  elliptic  co- 
ordinates (1)  may  be  employed  to  find  the  stresses  in  a  flat  plate  having  the 
form  of  a  hyperbolic  cylinder  of  any  eccentricity,  when  it  is.  subjected  to  a 
tensile  load  directed  along  the  imaginary  axis  of  the  hyperbola. 

Method  II.  (2)  may  be  used  to  find  directly  the  stresses  in  a  flat  plate  of  any 
shape,  to  which  given  edge  tractions  are  applied.  As  with  the  previous  method, 
concentrations  due  to  small  cylindrical  holes  or  cracks,  cut  normally  into  the 
surface  of  a  member,  may  be  found,  but  in  this  case  there  is  110  restriction 
on  the  shape  of  the  hole.  Further,  the  effect  of  cutting  any  scratch  or  groove, 
in  a  surface  along  which  the  principal  stresses  are  parallel  and  perpendicular 
to  the  direction  of  the  scratch  or  groove,  may  be  determined. 

There  is  another  possible  application  of  this  method  which,  so  far  as  the 
author  is  aware,  has  not  yet  been  used.  If  a  thin  plate,  initially  flat,  be  bent 
hy  couples  applied  at  its  edge,  the  equations  satisfied  by  the  principal  curvatures 
of  its  surface  are  mathematically  identical  with  those  satisfied  by  the  principal 
stresses  in  a  plate,  of  the  same  shape,  subjected  to  appropriale  edge  tractions. 
Hence,  if  a  photo-elastic  experiment  be  performed  in  which  the  applied  edge- 
tractions  represent,  on  some  convenient  scale,  the  prescribed  conditions  of 
curvature  at  the  boundary  of  the  bent  plate,  then  the  measured  stresses  at  any 
point  in  the  stretched  plate  will  represent,  on  the  same  scale,  the  curvatures  at 
the  corresponding  point  of  the  bent  .plate.  Hence  the  stresses  in  the  latter  may 
be  found. 

Turning  now  to  III.  (3),  it  is  possible  by  this  method  to  find  the  shearing 
stresses  in  a  bent  or  twisted  cylindrical  beam  or  shaft,  having  a  cross-section  of 
any  given  shape.  The  weakening  effect  of  any  small  groove,  formed  in  the 
surface  of  a  member  which  is  not  cylindrical,  may  be  found  in  those  cases  where 
the  stress  in  the  neighbourhood  is  a  shearing  stress  acting  on  planes  respectively 
parallel  and  perpendicular  to  the  groove. 

The  soap-film  solutions  for  small  surface  grooves  are  complementary  to  those 
obtained  by  the  photo-elastic  method,  so  that  a  combined  method  may  be  used 
to  find  the  effect,  of  such  a  groove  in  any  case  whatever. 

The  thermal  method  depends  on  the  fact  that  there  is  a  sudden  generation 
of  heat  in  soft  metals  such  as  mild  steel  when  the  stress  reaches  the  yield  point. 
The  resulting  rise  of  temperature  may  be  detected  by  means  of  a  thermo-couple 
applied  to  the  surface  of  the  material.  It  follows  that  advantage  may  be  taken 
of  this  phenomenon  to  determine  the  magnitude  of  stress  concentrations,  and 
that,  unlike  the  other  available  methods,  there  is  no  restriction  to  two-dimen- 
sional cases.  So  far,  however,  very  few  attempts  have  been  made  to  develop 
this  method. 

In  this  branch  of  the  subject  progress  is  most  urgently  needed  in  the  esti- 
mation of  stresses  in  solids  of  revolution.  On  the  purely  mathematical  side 
the  most  promising  line  of  attack  seems  to  be  the  investigation  of  solutions 
applicable  to  ellipsoids  and  hyperboloids  of  revolution,  and  it  may  be  regarded 
as  probable  that  an  important  advance  in  this  direction  will  shortly  be  made. 
As  regards^  solids  of  revolution  in  general,  what  is  required  is  some  method 
of  the  type'of  the  soap-film  and  photo-elastic  methods. 

3.  Examples  of  Stress  Concentrations. 

3.  Examples  of  Stress  Concentrations. — It  is  not  proposed  to  attempt  here 
anything  like  a  comprehensive  summary  of  the  detailed  work  which  has  been 
performed  by  the  foregoing  methods.  Nevertheless,  it  is  desirable  for  the 
purposes  of  the  present  paper  to  illustrate  by  means  of  examples  the  general 
nature  of  the  results  obtained. 
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In  what  follows  it  will  be  convenient  to  use  the  term  '  concentration  factor  * 
to  denote  the  ratio  of  the  true  calculated  stress  to  the  stress  which  would  be 
calculated  if  the  concentration  of  stress  were  to  be  neglected. 

A  very  simple  example  is  that  of  a  small  round  hole,  such  as  an  oil  hole, 
drilled  normally  into  the  surface  of,  say,  a  twisted  or  bent  circular  shaft.  In 
the  twisted  shaft  the  tensile  stress  concentration  factor  is  4,  while  the  shear 
stress  factor  is  2.  In  the  bent  shaft  the  tensile  and  shear  factors  are  both  3. 
These  results  are  practically  independent  of  the  ratio  of  the  size  of  the  hole  to 
the  size  of  the  shaft,  unless  the  former  is  so  large  as  to  reduce  materially  the 
cross-section  of  the  shaft. 

Similar  concentration  factors  are  found  in  cases  where  holes  are  drilled  in 
members  of  other  shapes. 

As  another  illustration  we  may  take  a  small  groove  of  semi-elliptic  cross- 
section,  cut  in  the  surface  of  a  twisted  shaft.  If  p  is  the  radius  of  curvature 
at  the  bottom  of  the  groove  and  a  is  the  depth,  tne  shear  stress  concentration 
factor  is 

1  + 

This  result  is  true  whatever  the  angle  between  the  direction  of  the  groove  and 
the  axis  of  the  shaft.     The  tensile  stress  factor,  on  the  other  hand,  varies  from 


if  the  groove  is  parallel  or  perpendicular  to  the  axis,  to 


if  the  groove  runs  round  the  shaft  in  the  form  of  a  45°  spiral. 

It  may  be  remarked  here  that  work  on  grooves  of  other  shapes  has  shown 
that  the  stress  concentration  depends  mainly  on  the  depth  and  the  radius  of 
curvature  at  the  bottom,  provided  the  groove  is  not  very  shallow.  For  instance, 
a  60°  Vee  groove  with  a  rounded  corner  gave  factors  only  a  few  per  cent,  less 
than  those  calculated  from  the  foregoing  formulae.  It  appears,  therefore,  that 
the  latter  may  be  used  even  when  the  groove  departs  considerably  from  the 
elliptic  form. 

It  will  be  seen  from  the  above  results  that  stress  concentrations  may 
theoretically  be  quite  large  in  cases  likely  to  occur  in  practice.  Thus,  if  the 
depth  o£  a  groove  is  sixteen  times  the  radius  of  curvature,  the  concentration 
factor  may  vary  from  5  to  9,  while  in  the  extreme  case  of  a  surface  crack  even 
these  values  may  be  greatly  exceeded. 

The  theoretical  stresses  are  lower  in  the  case  of  a  continuously  grooved 
surface  such  as  a  screwed  portion  than  in  an  isolated  groove.  For  example,  in  a 
particular  screwed  rod  the  tensile  stress  concentration  factor  was  estimated  to 
be  3.7,  while  for  an  isolated  groove  of  the  same  shape  the  factor  was  found  to 
be  4.9.  Possibly  this  fact  partly  explains  why  it  is  found  to  be  of  value  to  turn 
down  the  plain  portion  of  a  screwed  bolt  to  the  root  diameter ;  if  this  is  not 
done  the  stress  concentration  in  the  end  groove  is  considerably  greater  than  in 
the  rest  of  the  screwed  part. 

As  in  the  case  of  the  circular  holes,  the  above  results  are  substantially 
independent  of  the  size  of  the  grooves  unless  these  are  so. large  that  a  material 
amount  of  metal  is  removed  in  their  formation. 

A  striking  consequence  of  this  fact  is  that  severe  stress  concentrations  should 
arise  from  the  scratches  and  other  surface  defects  left  on  machine  parts  by  the 
oidinary  processes  of  manufacture.  As  will  be  seen  later,  the  non-fulfilment  of 
this  prediction  is  one  of  the  outstanding  discrepancies  between  the  theoretical 
and  experimental  aspects  of  our  subject. 

4.  Practical  Limitations  of  the  Foregoing  Methods. 

In  practice  it  is  not  infrequently  found  that  the  application  of  the  foregoing 
results,  on  the  basis  of  the  usual  criteria  of  failure,  gives  a  value  for  the  weaken- 
ing due  to  stress  concentration  which  is  misleading  or  entirely  wr<5ng.  In  other 
cases  the  agreement  is  quite  good  enough  for  practical  requirements.  It  is 
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therefore  important  to  find  out  in  what  way  the  assumptions  made  in  applying 
the  theory  are  defective. 

Methods  I.  and  III.,  above,  are  subject  to  the  usual  assumptions  of  the 
mathematical  theory  of  elasticity,  namely,  that  the  material  is  homogeneous  and 
isotropic  (or  possesses  some  particular  kind  of  seolotropy),  that  the  stresses  are 
proportional  to  the  strains,  and  that  the  strains  are  so  small  that  their  squares 
may  be  neglected.  As  regards  II.,  it  has  been  shown  that  this  method  measures 
the  stresses  actually  existing  in  the  xylonite  plates  which  are  usually  employed 
in  the  experimental  work.  Its  use  is  therefore  subject  only  to  the  assumption 
that  the  stress-strain  relations  of  the  material  considered  are  of  the  same  form 
as  those  of  xylonite. 

In  the  first  place,  it  is  to  be  noted  that  method  I.  gives  results  which  are 
in  substantial  agreement  with  those  of  method  II.  in  those  cases  where  a 
comparison  is  possible.  We  may  therefore  conclude  that  the  assumption  of 
linear  stress-strain  relations  involves  no  important  error  provided  the  departure 
from  linearity  is  of  the  same  order  as  that  of  xylonite  under  the  conditions  of 
photo-elastic  tests.  Similarly  we  note  that  the  assumption  of  infinitesimal  strains 
is  not  material  if  the  strains  are  comparable  with  those  in  xylonite  under  the 
above-mentioned  conditions.  Since  xylonite  exhibits,  in  these  two  respects, 
noticeable  departures  from  the  ideal  assumptions,  within  the  ranges  of  stress 
common  in  photo-elastic  tests,  it  may  be  concluded  that  the  effects  of  these 
assumptions  are  of  little  importance  in  practice  except  where  relatively  large 
strains  occur,  as  in  cleavage  slipping  and  viscous  flow. 

Further  consideration  of  the  causes  of  departure  from  theory  may  most  con- 
veniently be  made  with  reference  to  the  type  of  fracture  which  occurs.  For 
this  purpose  fractures  may  be  classified  as  elastic  or  brittle,  if  rupture  occurs 
without  material  inelastic  def ormaFion ;  plastic,  if  cleavage  slipping  occurs; 
viscous,  if  the  type  of  fracture  is  mainly  or  largely  determined  by  viscous 
flow;  and  fatigue  fractures.  The  viscous  type  will  not  be  further  discussed 
here. 

In  the  case  of  a  member  subjected  to  a  steady  load  such  that  the  fracture, 
when  it  occurs,  is  of  the  plastic  type,  it  has  long  been  known  that  stress  concen- 
trations have  little  or  no  weakening  effect.  The  reason  for  this  is  also  well 
known,  namely,  that  an  amount  of  plastic  flow,  at  the  point  of  high  'Stress,  which 
is  small  compared  with  the  amount  of  flow  at  fracture  is  adequate  to  annul  the 
stress  concentration.  In  some  experiments  described  in  a  recent  paper  (4),  the 
.author  showed  that  this  flow  occurs  at  approximately  the  estimated  load  even 
in  the  case  of  extremely  small  surface  scratches  having  a  depth  of  the  order 
of  10-4  inch. 

The  two  practical  cases  in  which  weakening  due  to  stress  concentration  is  of 
real  importance  are  those  of  elastic  fracture  under  a  load  steadily  applied,  and 
fatigue  failure  under  a  periodically  varying  load. 

Much  work  has  been  done  during  the  past  few  years  which  bears,  directly  or 
indirectly,  on  the  question  of  stress  concentrations  in  these  two  cases,  and  it 
may  now  be  said  with  considerable  confidence  that  one  reason  for  the  frequent 
discrepancy  between  the  estimated  concentration  factor  and  the  observed  weaken- 
ing is,  in  both  instances,  the  existence  of  a  scale  effect.  The  general  result  is 
that  stress  concentrations  are  practically  unimportant  if  the  linear  dimensions 
of  the  region  of  high  stress  are  sufficiently  small,  no  matter  what  the  estimated 
value  of  the  concentration  factor  may  be,  provided  that  it  does  not  exceed  a 
certain  high  upper  limit.  As  we  have  seen,  this  scale  effect  is  not  predicted 
by  the  purely  theoretical  methods  of  attack. 

As  regards  elastic  fractures,  the  scale  effect  appears,  so  far  as  existing 
knowledge  goes,  to  be  the  only,  important  cause  of  observed  departures 
from  theory.  In  the  case  of  fatigue  fractures,  on  the  other  hand,  it  may 
be  regarded  as  probable  that  redistribution  of  stress  due  to  cleavage  slipping 
often  occurs  under  ranges  of  stress  within  the  fatigue  limits  of  the  material. 
Under  these  circumstances  an  additional  factor  is  introduced  into  the 
problem. 

The  existing  evidence  on  the  above  points  is  discussed  in  the  three  succeed- 
ing sections. 
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5.  Scale  Effect :  (a)  Elastic  Fractures. 

With  certain  exceptions,  combined  stress  f  tests  which  result  in  -elastic 
fractures,  and  which  do  not  involve  stress  concentrations,  may  be  regarded  as 
satisfying  the  law  that  rupture  occurs  when  the  greatest  (positive)  tensile  stress 
reaches  a  particular  value.  The  exceptional  cases  are  those  in  which  the  greatest 
principal  compression  is  numerically  much  larger  than  the  greatest  principal 
tension. 

'Confining  ourselves  in  the  present  discussion  to  cases  where  the  greatest 
compression  is  not  predominant,  it  appears  that  we  may  reasonably  infer  that 
the  weakening  effect  of,  say,  a  scratch  should  be  obtainable  at  once  from  the 
tensile  stress  concentration  factor.  It  is  found,  however,  that  this  prediction 
is  by  no  means  verified  if  the  scratches  are  small. 

A  familiar  illustration  is  met  with  in  the  operation  of  glass-cutting,  where 
local  weakening  is  produced  by  means  of  a  scratch,  the  object  being  to  determine 
the  direction  of  fracture  when  the  glass  is  subsequently  broken.  It  is  found 
that  the  weakening  effect  is  insufficient  unless  the  depth  of  the  scratch  exceeds 
a  certain  value,  notwithstanding  the  fact  that  the  estimated  concentration  factor 
due  to  the  shallower  scratches  would  appear  to  be  quite  adequate  for  the  purpose 
in  view. 

As  an  example  of  the  same  phenomenon  in  elastic  fractures  of  crystalline 
metals,  some  experiments  on  hardened  cast-steel  may  be  cited. 

The  theory  states  that  if  surface  scratches  be  made  on  a  tensile  test  piece 
there  should  be  no  stress  concentration  in  the  case  of  scratches  whose  direction  is 
parallel  to  that  of  the  tensile  stress,  while  if  they  are  perpendicular  to  the  direc- 
tion of  the  stress  the  concentration  factor  should  be  a  maximum  for  the  particular 
shape  of  scratch.  Experiments  were  therefore  undertaken,1  at  the  author's 
request,  to  find  the  effect  of  the  orientation  of  the  surface  scratches  on  the  tensile 
breaking  load  of  dead  hard  cast-steel  test-pieces.  The  scratches  were  made 
with  No.  0  carborundum  cloth,  and  their  depth  was  of  the  order  of  10-4  inch. 
The  shapes  of  scratches  made  by  this  means  were  examined  micrographically, 
and  it  was  deduced  that  the  concentration  factor  due  to  the  perpendicular 
(circumferential)  scratches  should  be  3  to  4  at  least.  In  the  tensile  tests, 
however,  there  was  no  systematic  difference  between  the  breaking  loads  of  the 
axially  and  circumferentially  scratched  specimens.  Here,  then,  is  a  case  of  a 
severe  stress  concentration  which  appears  to  have  been  completely  annulled  by 
reason  of  the  scale  effect. 

In  another  experiment,  conducted  to  determine  the  effect  of  larger  grooves 
on  the  same  material,  two  strips  0.025  inch  thick  and  0.5  inch  wide  were  broken 
by  bending.  One  of  the  strips  had  a  serrated  edge,  the  depth  of  the  serrations 
being  0.025  inch.  It  was  found  that  the  breaking  couple  of  the  plain  strip  was 
2.48  times  that  of  the  serrated  one.  The  concentration  factor'  due  to  the 
serrations  was  estimated  to  be  2.9,  so  that  in  this  case  the  greater  part  of  the 
theoretical  weakening  was  developed. 

A  theory  of  scale  effect  in  relation  to  elastic  fractures  was  advanced  by  the 
author  in  the  paper  already  mentioned  (4).  It  was  shown  that,  on  the  usual 
assumptions  of  elastic  theory,  the  accepted  strengths  of  materials,  under  tests 
which  result  in  elastic  fracture,  are  incompatible  with  strengths  which  are 
deducible  on  theoretical  grounds  from  other  physical  properties  of  the 
materials,  the  theoretical  strengths  being  by  far  the  greater.  A  detailed  investi- 
gation performed  on  a  certain  kind  of  glass  indicated  that  this  discrepancy  could 
only  be  explained  by  supposing  that  the  material  contained  minute  flaws  which 
gave  rise  to  very  severe  stress  concentrations.  It  was  found  that  in  this  glass 
the  necessary  order  of  magnitude  of  the  flaws*  was  about  10-4  inch.  In  the  case 
of  vitreous  materials  a  process  was  discovered  whereby  these  flaws  could  be 
temporarily  eliminated,  and  the  substances  then  possessed  tensile  strengths 
agreeing  approximately  with  the  theoretical  values.  Thus  for  the  particular 
kind  of  glass  which  was  chiefly  used  the  tensile  strength  was  normally  about 
25,000  Ib.  per  sq.  in.  After  treatment  values  as  high  as  900,000  Ib.  per  sq.  in. 
were  recorded,  so  that  the  concentration  factor  due  to  the  flaws  was  about  36. 

These   results  suggest  immediately  an  explanation  of  the  scale  effect.     If  a 

1  By  Mr.  W.  D.  Douglas,  of  the  Eoyal  Aircraft  Establishment. 
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small  surface  scratch  be  formed  whose  depth  is  of  the  same  order  as  the 
dimensions  of  the  flaws  already  existing,  but  which  gives  rise  to  a  less  severe 
stress  concentration,  the  strength  of  the  piece  must  clearly  be  unaffected,  since 
whether  the  scratch  be  present  or  not,  rupture  must  be  occasioned  by  the  stress 
concentration  due  to  the  Haws.  If,  however,  a  groove  be  made,  of  such  a  size 
that  the  region  of  high  stress  is  large  compared  with  the  size  of  the  flaws,  there 
must  be  a  double  magnification  of  stress,  arising  firstly  from  the  concentration 
of  stress  at  the  bottom  of  the  groove,  and  secondly  from  the  flaws  within  that 
region  of  high  stress. 

Since  the  above-mentioned  paper  was  published,  it  has  been  found  possible 
to  prepare  pure  vitreous  silica  in  a  stable  form  in  which  a  relatively  large  pro- 
portion (at  least  one-third)  of  the  theoretical  tenacity  is  retained.  According  to 
the  theory,  this  implies  that  the  flaws  must  be  reduced  almost  to  molecular 
dimensions,  whence  it  follows 'that  the  scale  effect  should  be  practically  non- 
existent. 

This  conclusion  hns  received  ample  support  from  a  number  of  experiments, 
the  results  of  which  are  sufficiently  illustrated  below. 

Two  thin  circular  rods  of  the  material  were  touched  together  as  lightly 
as  possible,  whereby  a  minute  injury,  invisible  under  a  magnification  of 
250  diameters,  was  inflicted  on  the  surface  of  each.  One  of  the  rods  was  broken 
by  flexure,  with  the  abrasion  on  the  tension  side.  The  degree  of  flexure  at 
rupture  was  only  about  one-eighth  of  that  required  to  break  an  uninjured  rod 
of  the  same  diameter.  The  other  rod  was  broken  with  the  injury  on  the  com- 
pression side,  and  in  this  case  no  weakening  could  be  detected.  This  is  in 
agreement  with  the  hypothesis  of  rupture  under  a  specific  tensile  stress,  since  no 
concentration  of  tensile  stress  should  arise  from  a  surface  defect  on  the  com- 
pression side  of  a  bent  beam.  Further  experiments  were  performed  in  which 
rods  were  lightly  touched  with  other  solid  bodies,  both  hard  and  soft,  and  with 
the  finger  tip.  In  all  cases  a  notable  weakening  resulted.  Even  if  a  rod  was 
left  exposed  to  bombardment  by  the  dust  particles  in  the  atmosphere  there  was  a 
slow  but  perfectly  definite  weakening. 

In  the  case  of  brittle  substances,  then,  it  appears  that  stress  concentrations 
due  to  minute  flaws  and  surface  defects,  far  from  being  negligible,  constitute 
the  controlling  factor  which  determines  the  magnitude  of  the  technically  avail- 
able tenacity  of  these  materials,  a  result  which  may  be  regarded  as  a  remarkable 
vindication  of  the  elastic  theory  of  stress  concentrations. 

6.  Scale  Effect  =  (b)  Fatigue  Fractures. 

The  general  nature  of  available  experimental  results  in  this  branch  of  the 
subject  is  the  same  as  in  the  case  of  elastic  fractures,  that  is  to  say,  fine  scratches 
appear  to  have  little  or  no  effect,  while  sufficiently  large  grooves  often  give  rise 
to  nearly  the  full  theoretical  weakening.  The  evidence  is,  however,  less  com- 
plete, and  much  more  work  is  required  to  place  our  knowledge  of  scale  effect  in 
fatigue  phenomena  on  a  satisfactory  basis. 

On  the  theoretical  side  the  main  difficulty  arises  from  the  lack  of  a  well- 
established  theory  of  fatigue  phenomena  in  general.  Beilby's  original  work  (5) 
on  the  production  of  an  amorphous  phase  in  metals  by  overstrain  showed  that 
large  internal  stresses  might  be  set  up  as  a  result  of  cleavage  slipping,  by  reason 
of  the  difference  in  density  between  the  amorphous  and  crystalline  phases.  This 
result  suggested  as  a  possible  theory  of  fatigue  that  the  internal  stresses  set  up 
by  repeated  cleavage  slipping,  when  combined  with  the  stress  system  externally 
applied,  might  ultimately  be  adequate  to  initiate  a  brittle  fracture  of  the  material. 
Beilby  considered  that  the  amorphous  phase  was  developed  between  the  slipped 
surfaces,  so  that  on  the  theory  derived  from  his  work  fatigue  fracture  should 
always  be  an  ultimate  effect  of  repeated  slipping.  A  similar  conclusion  may  be 
reached  regarding  the  theory  of  Ewing  and  Humfrey  (6),  according  to  which 
fatigue  cracks  are  initiated  by  the  attrition  of  the  slipping  cleavage  surfaces. 

A  difficulty  in  connection  with  both  these  theories  is  that,  as  the  present 
author  has  shown  (see  §4  above)  cleavage  slipping  can  occur  at  the  corners  of 
very  small  scratches  at  loads  very  much  below  those  necessary  to  cause  it  to  take 
place  in  the  remainder  of  the  material.  According  to  the  foregoing  theories, 
therefore,  the  presence  of  such  scratches  on  a  test-piece  should  greatly  reduce 
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its  fatigue  strength,  a  conclusion  which  is  not  borne  out  by  the  results  of 
tests. 

Another  view  was  put  forward  by  the  author  (4),  according  to  which  it  was 
considered  that  the  slipped  surfaces  suffered  no  permanent  change,  and  that  the 
amorphous  phase  was  generated  entirely  at  the  inter-crystalline  boundaries.  On 
this  theory,  no  phase  change,  and  therefore  no  fatigue  cracking,  can  result  from 
cleavage  slipping  within  a  'crystal  unless  the  slipping  extends  at  least  as  far 
as  the  junction  with  neighbouring  crystals.  Hence  this  theory,  unlike  the  other 
two,  does  predict  an  effect  similar  to  the  observed  scale  effect,  and  to  that  extent 
is  more  satisfactory.  It  does  not  seem  likely,  however,  that  the  author's  theory 
will  prove  entirely  adequate  in  its  present  form,  as  many  new  and  so  far  un- 
explained facts  have  come  to  light  since  the  publication  of  the  paper  in  which 
it  was  first  advanced.  For  instance,  some  data  now  available  suggest  on  the  one 
hand  that  fatigue  fractures  can  sometimes  occur,  without  cleavage  slipping,  and 
on  the  other  that  repeated  cleavage  slipping  in  certain  materials  does  not  neces- 
sarily lead  to  fatigue  failure  even  though  it 'extend  over  a  number  of  crystals. 

On  the  experimental  side  the  position  of  this  branch  of  our  subject  appears 
to  be  even  less  satisfactory  than  on  the  theoretical  side.  So  far  as  the  author 
is  aware,  very  few  published  results  exist  which  are  sufficiently  complete  for  a 
close  comparison  with  theory ;  essential  data,  such  as  the  shape  of  the  scratches 
or  the  size  of  the  crystals,  being  omitted. 

7.  Fatigue  Fractures :  Effect  of  Stress  Redistribution. 

It  may  be  regarded  as  proved  that,  in  many  metals,  no  important  plastic  flow 
occurs,  within  the  fatigue  limits  of  stress,  under  an  applied  load  alternating 
between  numerically  equal  positive  and  negative  values.  Under  these  conditions 
the  only  uncertainty  in  applying  the  theory  of  stress  concentrations  is  that  due 
to  the  scale  effect. 

The  case  is  otherwise,  however,  in  the  technically  commoner  examples  of 
fatigue  failure  in  which  the  alternating  load  is  superposed  on  a  relatively  large 
steady  load.  Fatigue  under  such  loads  is,  in  general,  preceded  by  plastic  flow, 
and  in  dealing  with  failure  at  a  point  of  local  high  stress,  such  flow  involves  a 
redistribution  of  stress  in  the  neighbourhood.  The  nature  of  this  redistribution 
cannot  at  present  be  found  by  calculation.  It  follows  that  while  the  range  of 
stress  at  the  fail  point  is  calculable  the  mean  stress  is  not ;  there  is  therefore 
no  basis  on  which  the  ordinary  data  regarding  the  fatigue  strength  of  the  material 
can  be  employed  to  find  the  conditions  of  rupture  at  a  point  of  stress  concentra- 
tion under  such  kinds  of  load. 

In  the  case  of  some  metals,  some  recent  experiments  suggest  that  plastic  flow 
may  be  met  with  even  under  a  pure  alternating  stress. 

It  may  be  remarked  that  the  effect  discussed  above,  though  a  serious  obstacle 
to  theoretical  work,  is  probably  of  great  value  in  practice,  since  the  redistribution 
of  stress  must  decrease  the  mean  stress,  and  therefore  the  weakening  effect  of  fhe 
concentration. 

An  additional  complication  accrues  from  the  fact  that  if  appreciable  plastic 
flow  occurs  at  the  fail  point  the  material  there  must  be  in  a  cold-worked  condition, 
whence  increased  resistance  to  fatigue  is  to  be  anticipated. 

8.  Stress  Concentration  in  Mechanical  Testing. 

It  is  customary,  in  the  usual  mechanical  tests,  including  fatigue  tests,  to  make 
the  working  portion  of  the  test-piece  of  smaller  section  than  the  portions  which 
are  gripped,  in  order  to  localise  the  fracture  and  other  effects  of  the  test.  The 
junction  of  the  two  parts  is  rounded  off  so  as  to  secure  a  gradual  change  of 
section.  Now,  at  such  a  change  of  section  there  must,  in  general,  be  a  stress 
concentration,  so  that  the  actual  maximum  stress  must  be  greater  than  the 
estimated  working  stress  in  the  specimen ;  moreover,  it  is  not  usually  possible 
at  present  to  calculate  the  magnitude  of  the  concentration  factor.  It  appears, 
then,  that  if  the  difference  between  the  estimated  stress  and  the  true  maximum 
stress  is  appreciable  the  ostensible  results  of  the  tests  may  be  materially 
erroneous,  especially  in  the  case  of  fatigue  tests.  It  is  evidently  necessary  in  any 
important  test  work  to  reduce  this  unknown  factor  to  negligible  dimensions. 
Partly  on  theoretical  grounds  and  partly  as  the  result  of  experience  with  test 
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pieces  of  various  forms,  the  author  has  reached  the  conclusion  that  this  object 
may  be  attained  in  most  cases  by  making  the  radius  of  the  fillet  at  the  change 
of  section  at  least  five  times  the  working  diameter  (or  corresponding  dimension) 
of  the  test-piece. 

The  only  other  type  of  stress  concentration  which  is  likely  to  modity  the 
results  of  mechanical  tests  is  that  due  to  surface  scratches.  As  has  been  seen, 
advantage"  may  be  taken  of  the  scale  effect  to  render  this  factor  unimportant,  by 
making  the  scratches  sufficiently  fine.  Scratches  having  a  depth  of  the  order  of 
10- '  inch,  such  as  are  produced  by  No.  0  emery  cloth,  or  by  grinding  with  a 
wheel  of  about  No.  80  grit,  would  appear  to  be  satisfactory  in  all  ordinary  cases. 

9.  Allowance  for  Stress  Concentrations  in  Design. 

The  author  has  often  been  asked  by  designers  what  measures  can  be  taken  to 
guard  against  the  occurrence  of  'dangerous  stress  concentrations  in  practice.  It 
is  clear  from  what  has  already  been  said  that  nothing  like  a  complete  answer  to 
this  question  is  possible  at  present.  All  that  can  be  done  is  to  lay  down 
approximate  rules  based  not  only  on  calculable  cases  and  such  experimental  facts 
as  have  so  far  been  collected,  but  also  very  largely  on  personal  judgment.  The 
subjoined  recommendations,  which  embody  advice  which  the  author  has  given 
from  time  to  time  in  the  past,  may,  perhaps,  be  found  useful  pending  the  develop- 
ment of  more  rational  methods.  It  may  be  remarked  that  they  have  not  infre- 
quently been  successful  in  overcoming  difficulties  encountered  in  practice.  They 
are  intended  to  apply  mainly  to  members  possessing  some  ductility  and  subjected 
to  periodically  varying  loads. 

In  the  first  place,  the  radius  of  the  fillet  in  a  re-entrant  corner  should  be  at 
least  a  quarter  of  a  certain  dimension  which  may  be  called  the  ruling  dimension. 
For  giooves,  such  as  screw-threads  and  keyways,  the  ruling  dimension  is  the  t 
depth  of  the  groove.  Alternatively,  in  a  continuously  grooved  surface,  such  as 
occurs  in  a  gear  wheel,  the  root  thickness  of  the  teeth  or  ridges  may  be  used  if 
this  is  less  than  the  depth.  For  sudden  bends,  as  in  crankshafts,  the  ruling 
dimension  may  be  taken  to  be  the  radius  of  the  shaft  (or  the  corresponding 
dimension).  At  a  change  of  section  on  a  shaft,  the  smaller  radius  or  the  differ- 
ence between  the  radii,  whichever  is  the  less,  may  be  taken. 

The  above  values  of  the  radius  of  fillet  should  be  regarded  as  minima  and 
should  be  exceeded  wherever  possible. 

In  addition,  an  allowance  should  be  made  for  the  weakening  effect  of  the 
stress  concentration.  It  is  impossible  at  present  to  give  any  very  definite  idea 
of  the  magnitude  of  this  factor,  as  it  appears  to  vary  so  greatly  with  the  nature 
of  the  material  and  the  kind  of  load.  With  the  worst  materials  and  pure 
alternating  loads  it  may  be  necessary  to  allow  the  full  theoretical  concentration 
factor,  which  with  the  above  values  of  the  radius  of  fillet  will  usually  be  about 
2  to  6.  In  more  favourable  cases,  e.g.,  where  the  load  fluctuates  but  does  not 
change  sign,  much  lower  values  may  be  used.  Probably  a  factor  of  2  will  be 
found  satisfactory  in  the  majority  of  such  cases. 

As  has  been  seen  above,  plastic  flow  can  be  of  great  value  in  reducing  the 
effect  of  local  stress  intensifications.  Hence  two  qualities  are  desirable  in  the 
material  of  a  member  liable  to  break  at  a  point  of  stress  concentration.  In  the 
first  place,  the  fatigue  limits  should  be  as  wide  as  possible,  and  in  the  second 
the  range  of  stress  within  which  plastic  flow  is  absent  should  be  as  small  as 
possible. 
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V. 

The  Strain  Energy-Function  and  the  Elastic  Limit. 

By  Professor  B.  P.  HAIGH,  D.Sc.,  M.B.E.,  Royal  Naval  College,  Greenwich. 

In  a  contribution  to  the  report  of  this  Committee,  for  1919,  the  author 
analysed  published  data,  with  the  object  of  comparing,  for  ductile  metals, 
the  different  elastic  limits  under  simple  and  complex  stresses.  Experimental 
values  of  the  ratios  between  the  elastic  limits  for  different  kinds  of  stress 
were  compared  with  the  corresponding  values  predicted  by  applying  three 
alternative  hypotheses,  due  to  Lame  and  Rankine,  to  Darwin,  Tresca  and 
Guest,  and  to  Saint  Venant ;  and  with  values  deduced  from  a  noval  hypothesis, 
viz.  that  the  relation  between  the  elastic  limits  is  governed  by  the  strain  energy 
per  unit  volume  which,  at  the  elastic  limit,  reaches  a  definite  limiting  value 
independent  of  the  simple  or  complex  nature  of  the  applied  stress. 

Diagrams  were  plotted  in  which  the  marked  points  represented  experimental 
determinations  of  the  ratios  between  the  elastic  limits.  On  comparing  the  loci 
of  these  experimental  points  with  the  graphs  representing  the  alternative 
hypotheses,  it  was  observed  that  the  strain-energy  graph  was  in  fair  agreement 
with  experiment  throughout  the  field  of  investigation,  and  that  the  graphs 
representing  other  hypotheses  were  in  agreement  in  narrower  fields,  e.g.  in  the 
case  of  Saint.  Venant' s  hypothesis,  when  the  ratio  between  the  principal  stresses 
was  less  tharT+  0.30. 

The  diagram  shown  in  fig.  13  illustrates  the  method  of  comparison  adopted 
for  two-dimensional  stresses;  the  axes  OX  and  OY  representing  the  principal 
stresses  on  the  material  and  the  lengths  OA  and  OB  the  elastic  limits  in  simple 
tension.  The  co-ordinates  of  points  in  the  quadrants  represent  the  principal 
stresses  for  complex  combinations.  The  strain-energy  hypothesis  is  represented 
by  an  ellipse  whose  eccentricity  varies  slightly  for  different  values  of  Poisson's 
Ratio  <r.  Other  hypotheses  are  represented  by  figures  composed  of  intersecting 
straight  lines.  The  diagram  includes  also  a  fifth  hypothesis,  published  in  1916 
by  Dr.  Albert  Becker,1  viz.  that  the  elastic  limit  is  determined  by  a  dual 
condition — limiting  maximum  shear  stress  and  limiting  maximum  strain.  This 
is  represented  by  a  ten-sided  figure  which,  for  combinations  approximating  to 
shear  stress,  nearly  coincides  with  the  ellipse  for  the  strain-energy  hypothesis. 
Where  two  nearly  equal  like  stresses  are  combined,  as  in  turbine  discs,  Dr. 
Becker's  hypothesis  appears  to  overestimate  the  elastic  limit,  although  not  so 
greatly  as  does  Saint  Venant's. 

In  the  earlier  report  it  was  explained  that  the  strain-energy  hypothesis 
represents  not  merely  an  arbitrary  assumption,  roughly  endorsed  by  published 
data,  but  an  attempt  to  apply  established  thermodynamic  principles  to  the 
now  generally  accepted  theory  that  the  production  of  non-elastic  strain  is 
associated  with  a  change  in  the  state  of  the  metal,  from  the  crystalline  to  the 
*  vitreous  '  or  '  amorphous '  phase.  The  earlier  report,  however,  was  chiefly 
an  analysis  of  published  data,  without  any  attempt  at  explanation.  In  what 
follows,  the  object  is  to  trace  the  theoretical  bearing  of  thermodynamic  prin- 
ciples on  the  relation  between  the  elastic  limits  under  different  kinds  of  stress. 

Physical  Theory  of  Non-elastic  Strain. 

The  current  theory  for  permanent  strain  rests  on  two  main  experimental 
observations  :  (1)  that  permanent  strain  is  the  cumulative  result  of  numerous 
gliding  movements  in 'individual  grains  of  metal,  the  reality  of  these  movements 
being  demonstrated  by  the  slip-bands  observed  when  polished  faces  of  a  block 
strained  beyond  the  elastic  limit  are  examined  under  the  microscope ;  and 

1  A.  J.  Becker,  Bulletin  85,  University  of  Illinois. 
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(2)  the  development  of  permanent  strain  is  associated  with  a  change  of  physical 
state,  from  the  crystalline  to  the  undercooled  liquid  or  '  vitreous  '  state.  In 
the  developed  theory  the  change  from  the  elastic  crystalline  state  to  the  vitreous 
is  held  to  be  essential  for  the  explanation  of  the  gliding  movement. 

DIAGRAM  FOB  ELASTIC-LIMIT  :  Two  PRINCIPAL-STRESSES. 


FIG.  13. 
Strain-energy  hypothesis  is  represented  by  ellipse  (Equation  X2-f  Y2— 2cr.X.Y=OA72 

Poisson's  ratio  <r=-). 

3 

Rankine  hypothesis  for  max.  principal-stress,  figure  A1B2A'3B'4. 
St.  Venant  hypothesis  for  max.  principal -strain,  figure  A5B6A'7B'8. 
Guest's  hypothesis  for  max.  shear-stress,  figure  A1BA'3B'. 
Becker  hypothesis  for  max.  shear-stress  and  strain,  figure  A9,  10,'  11,  B12,  13,  A'. 

Before  this  theory  was  advanced  by  Beilby,  and  developed  by  Ewing, 
Kpsenhain.  and  others,  there  was  no  satisfactory  explanation  of  permanent 
strain.  The  difficulty  is  clearly  defined  in  Professor  James  Thomson's  state- 
ment, written  in^  1861  2  :  '  I  have  hot  any  conception  of  continuous  crystalline 
structure  admitting  of  what  may  be  called  ductile  bending — that  is,  bending 

2  James  Thomson,  Proc.  Roy.  Soc.,  1861. 
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beyond  the  limits  of  elasticity  .  .  .  and  still  remaining  of  the  nature  of  one 
Continuous  crystal.  What  .  .  .  may  be  the  nature  of  the  change  of  molecular 
arrangement  induced  by  bending  them  I  cannot  say;  but  I  suppose  that,  in 
their  yielding,  their  crystalline  structure  is  materially  altered,  and  rendered 
discontinuous  where,  before,  it  was  continuous.' 

Where  Thomson  used  the  word  '  discontinuous,'  current  theory  substitutes 
'vitreous'  or  '  amorphous,'  .with  many  associated  ideas  and  more  or  less 
satisfactory  definition.  The  physical  and  chemical  properties  of  vitreous  metals 
are  still  in  some  measure  uncertain,  because  the  vitreous  phase  tends  to 
recrystallise ;  but  it"  is  known  that  the  physical  characteristics  include  great 
hardness,  the  limited  mobility  of  a  highly  viscous  fluid,  and,  in  the  case  of 
iron,  low  magnetic  permeability.  The  change  is  not  merely  a  matter  of 
pulverising.  Without  entering  on  controversial  matter  that  has  collected  round 
the  definition  of  an  '  allotropic '  change,  it  may  be  accepted  that  the  change 
•that  occurs  when  a  ductile  metal  is  strained  is  purely  physical  in  the  sense 
that  it  involves  only  a  rearrangement  of  the  molecules  without  change  of 
internal  molecular  structure. 

If  this  theory  be  accepted  as  an  explanation  of  permanent  strain,  it  is 
evident  that  the  change  from  the  crystalline  to  the  vitreous  state  must  precede 
the  gliding  action ;  the  essential  condition  for  gliding  movement  cannot  be 
merely  a  consequence.  Without  denying  the  probability  that  further  quantities 
of  metal  may  suffer  the  same  change  during  the  gliding  movement,  it  is  clear 
that  the  initial  change  of  state  must  be  a  direct  result  of  the  preceding  elastic 
stage  of  straining.  It  follows  that,  in  the  course  of  any  alternative  process 
by  which  the  metal  can  be  strained  to  its  different  elastic  limits  by  the  applica- 
tion of  different  stresses,  the  metal  that  suffers  the  change  must  absorb — as  heat 
or  as  mechanical  work — such  quantities  of  work  as  enable  it  to  change  from  the 
stable  crystalline  state  to  the  vitreous  phase,  which,  at  ordinary  temperatures, 
is  known  to  be  '  metastable,'  i.e.  unstable  but  restrained  from  change  by 
'  internal  viscosity.' 

Thermodynamic  Principles. 

When  any  physical  or  chemical  change  is  produced  mechanically,  by  forces 
exerted  by  external  bodies  acting  on  a  mass  of  the  substance  in  question 
(e.g.  by  a  piston  acting  on  vapour  enclosed  in  a  compressor  cylinder),  the  work 
done  by  'the  forces  during  the  change  is  ordinarily  greater  than  the  quantity 
that  can  be  regained  by  allowing  the  change  to  occur  in  the  reverse  direction 
at  the  same  temperature.  The  difference  between  the  two  quantities  is  con- 
verted to  heat  by  the  action  of  friction  or  other  '  irreversible  '  effects.  Under 
ideal  frictionless  conditions  the  two  quantities  of  work  may  be  equal,  and  the 
process  of  change  is  then  said  to  be  '  reversible.' 

A  well-known  application  of  the  second  law  of  thermodynamics  states  that 
when  a  given  change  can  be  produced  by  different  reversible  processes,  carried 
out  at  the  same  temperature,  the  quantity  of  work  that  must  be  done  by  the 
forces  acting  on  the  substance  must  be  definitely  constant,  i.e.  independent 
of  how  the  forces  act.  If  one  and  the  same  change  of  physical  state  can  be 
caused  by  stresses  of  different  kinds,  pull  or  shear  or  combinations  of  these, 
the  quantities  of  work  done  by  these  forces  on  unit  mass  of  substance  suffering 
the  change  must  always  be  the  same,  provided  that  the  action  is  reversible 
and  is  carried  out  at  the  same  temperature.  If  the  action  be  not  perfectly 
reversible,  somewhat  more  work  may  have  to  be  done ;  but  the  quantity  will 
approximate  to  the  ideal  constant  value  if  the  action  is  nearly  reversible. 

It  is  clear,  that  this  thermodynamic  law  affords  a  theoretical  relation  between 
the  elastic  limits  of  a  ductile  metal.  To  develop  such  a  relation,  however,  it 
is  necessary  to  equate  expressions  for  the  work  done  by  different  forces  in 
changing  unit  mass  of  metal  from  the  crystalline  to  the  vitreous  state  in 
reversible  and  isoMiermal  processes. 

Prof.  James  Thomson's  Work  on  'Regelation.' 

In  1849.3  in  anolying  thermodynamic  principles  to  problems  involving  change 
from  the  crystalline  to  the  fluid  state,  Professor  James  Thomson  showed  that 
the  melting-point  of  ice  should  be  lowered  by  the  application  of  fluid  pressure  : 

8  James  Thomson,  Trans.  Roy.  Soc.,  Edinburgh,  1849. 
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and  in  the  following  year  the  theorem  was  verified  in  experiments  carried  out 
by  Lord  Kelvin,  then  Professor  William  Thomson. 

Considering  the  Carnot  cycle  of  operations  in  an  ice-water  heat  engine 
operating  in  a  small  range  of  temperature  near  the  freezing-point,  Thomson 
showed  that  the  work  done  in  a  cycle  using  a  difference  of  pressure  dp,  acting 
through  the  small  change  of  volume  (V  —  v)  that  occurs  when  ice  changes  to 
water,  could  be  equated  to  the  work  done  by  an  ideal  heat  engine  receiving 
the  latent  heat,  L,  of  the  water  produced  and  working  through  a  range  of 
temperature  equal  to  the  corresponding  depression  of-  the  freezing-point.  Thus  : 


dp  (V-t>)  = 
whence  dT/dp  =  (T/L)  (V  -  v) 

where  dT/dp  is  the  rate  of  fall  of  the  freezing-point  with  rise  of  pressure. 

In  establishing  this  expression,  Thomson  considered  only  the  influence  of 
fluid-pressure,  acting  equally  in  all  directions  ;  but  in  a  later  paper,  of  1861,* 
the  application  of  thermodynamic  principles  was  greatly  widened.  He  then 
states  :'...!  soon  positively  formed  the  opinion  that  any  stresses  whatever, 
tending  to  change  the  form  of  a  piece  of  ice  in  ice-cold  water  (whether  these 
stresses  be  of  the  nature  of  pressures  or  tensions  —  that  is,  -pulls  or  pushes  —  and 
whether  they  be  in  one  direction  alone,  or  in  more  directions  than  one),  must 
impart  to  the  ice  a  tendency  to  melt  away  and  to  give  out  its  cold,  which 
will  tend  to  generate,  from  the  surrounding  water,  an  equivalent  quantity  of 
ice  free  from  the  applied  stresses.'  The  italics  are  Thomson's,  and  appear 
noteworthy  in  relation  to  the  straining  of  metals.  It  would  seem  that  the 
action  of  gliding  may  be  not  unlike  the  action  pictured  by  Thomson. 

It  does  not  appear  that  Thomson  proceeded  to  apply  to  metals  the  wider 
principle  enunciated  in  1861  ;  and  it  has  often  been  pointed  out  that  the  regelation 
theorem  of  1849  does  not  in  itself  suffice  as  an  explanation  of  the  change  of 
state  that  occurs  in  metals.  Unlike  water,  most  metals  show  a  slight  expansion 
when  they  change  from  the  crystalline  to  the  fluid  form. 

The  Thermodynamic  Reversibility  of  the  Stages  of  Strain. 

In  applying  the  second  law  of  thermodynamics  to  establish  any  relation 
between  different  forces  that  can  effect  a  specified  change,  it  is  essential  first 
to  define  ideal  conditions  such  that  the  process  of  change  may  be  reversible 
and  isothermal.  No  real  process  being  strictly  reversible  under  working 
conditions,  it  is  inevitable  that  some  tax  is  made  on  the  imagination. 

The  process  of  permanent  strain  may  be  regarded  as  comprising  two  distinct 
stages  :  (a)  the  initial  stage,  in  which  the  metal  is  strained  elastically  until 
the  first  molecule  leaves  the  '  continuous  lattice  '  of  the  crystal  to  enter  the 
'  discontinuous  assembly  '  of  the  vitreous  phase;  and  (b)  the  subsequent  stage, 
in  which,  in  consequence  of  the  changed  physical  properties  of  the  metal, 
gliding  movements  are  produced  by  any  shear-stress  that  is  applied. 

It  is  not  suggested  that  the  second  stage  approximates  to  a  reversible 
process.  Usually  pictured  as  viscous  gliding,  but  probably  more  complex  in 
action  and  possibly  akin  to  the  process  of  melting  and  recrystallisation  pictured 
by  Thomson's  description  of  the  action  in  the  case  of  ice,  the  action  certainly 
results  in  large  quantities  of  work  being  converted  to  heat,  raising  the  tempera- 
ture of  the  metal.  Such  an  irreversible  action  affords  no  basis  for  an  application 
of  the  second  law  of  thermodynamics. 

It  is  important  to  note,  however,  that  this  irreversible  gliding  movement 
is  only  a  consequence  of  the  preceding  change  of  state,  and  comparable  with 
other  consequences,  such  as  the  recognised  changes  in  thermal  and  electrical 
conductivity,  magnetic  permeability,  and  chemical  activity.  The  motion  that 
produces  permanent  strain  is  possible  only  when  the  change  of  state  has  been 
produced  by  an  earlier  cause. 

The  initial  stage,  on  the  other  hand,  exhibits  the  characteristics  of  a  reversible 
process,  although,  admittedly,  experimental  reversal  may  be  difficult  of 
demonstration.  Considerations  bearing  on  this  point  may  be  summarised  as 
follows  :  (1)  Although  the  work  done  in  straining  the  metal  up  to  and  beyond 

4  James  Thomson,  Proc.  Roy.  Soc.,  London,  1861. 
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the  elastic  limit  cannot  be  wholly  regained  by  releasing  the  applied  stresses, 
there  is  evidence  to  show  that  part  of  the  difference  may  be  regained  by 
other  methods.  'For  example,  since  the  electrolytic  potentials  of  the  strained 
and  unstrained  metals  are  different,  a  quantity  of  work  may  be  regained  by 
using  electrodes  of  the  two  in  a  suitable  cell.  When  a  given  mass  of  vitreous 
metal  is  replaced  by  an  equivalent  mass  of  crystalline,  a  quantity  of  electrical 
work  is  obtained,  and  any  deficiency  —  such  as  would  doubtless  be  found  in 
experiment  —  may  fairly  be  regarded  as  the  equivalent  of  heat  generated  during 
the  gliding  motion  of  the  second  stage  of  straining.  (2)  While  the  great  internal 
viscosity  of  the  vitreous  metal  prevents  recrystallisation  at  ordinary  tempera- 
tures, the  process  may  be  accelerated  by  raising  the  temperature  (above  400°  C. 
in  the  case  of  steel).  (3)  General  experience  in  the  widest  fields  of  research 
indicates  that  every  change  of  physical  or  chemical  state  is  inherently  reversible, 
actually  or  by  means  of  appropriate  imaginary  methods  of  operation.  (4)  In 
straining  the  metal  elastically,  up  to  the  elastic  limits  at  which  the  change 
occurs  under  different  complex  stresses,  quantities  of  work  are  done  on  the 
metal  in  strictly  reversible  manners  ;  and  these  quantities  may  be  expressed 
in  terms  of  the  applied  stresses  and  the  elastic  constants  of  the  metal. 

Expression  for  the  work  done  in  the  Reversible  Stage  of  Strain. 

The  reversible  process  now  considered  is  the  initial  stage  of  strain,  com- 
mencing when  the  crystalline  metal  —  initially  free  from  strain  —  is  subjected  to 
a  gradually  increasing  stress,  and  terminating  when  the  first  molecules  are 
projected  out  of  the  'crystalline  lattice'  into  the  'vitreous  assembly.'  It  is 
evident  that  the  process  can  be  carried  out  isothermally,  thus  completing  the 
conditions  required  for  the  application  of  thermodynamic  principles.  We  have 
to  attempt  the  problem  of  obtaining,  in  terms  of  the  applied  stresses,  an 
expression  for  the  work  done  on  unit  mass  of  metal  suffering  the  change  ;  or, 
alternatively,  an  expression  proportional  to  this  quantity. 

If  the  whole  mass  suffered  the  change,  and  were  isotropic,  and  if  the  change 
involved  no  alteration  of  volume  other  than  the  elastic  compression  or 
expansion,  the  work  done  by  the  applied  forces  would  be  identical  with  the 
elastic  strain-energy  at  the  elastic  limit  ;  and  would  be  given  by  the  expression 


,  —       _  (YZ+ZX+XY) 

2E  E 

Part  of  this  quantity  of  work,  '  however,  would  be  stored  in  the  vitreous  metal 
produced,  in  virtue  of  its  elastic  compression  or  expansion  under  the  applied 
stresses.  Thus  the  nett  work  corresponding  to  the  change  of  state  would  be 


where  K    denotes  the   (unknown)  modulus  of  compressibility  of  the  vitreous 
metal. 

On  the  assumption  that  the  change  results  in  a  change  of  volume,  in  ratio  a 
as  measured  free  from  stress,  a  further  term  may  be  introduced  —  analogous  to 
the  .  expression  used  in  Thomson's  '  regelation  '  theorem  :  Thus 


The  relation  between  the  elastic  limits  would  then  be  expressed  by  an  equation 
stating  that  the  quantity  W"  is  constant,  i.e.  independent  of  the  ratios  between 
the  stresses  X,  Y,  and  Z. 

The  problem  is  appreciably  complicated  by  the  non-isotropic  character  of 
actual  crystalline  metal  ;  and  by  the  fact  that  only  a  small  proportion  of  the 
total  mass  actually  suffers  the  change  of  state.  On  the  probable  assumption  that 
the  increase  of  volume  a  is  only  slight,  the  relative  importance  of  the  terms 
introduced  for  W'  and  W"  largely  disappears  ;  or  on  the  not  improbable 
hypothesis  that  the  change  involves  a  reduction  of  volume,  the,  two  terms  may 
wholly  disappear,  leaving  the  strain-energy  W  as  the  constant  governing  factor. 
On  this  latter  hypothesis  we  may  picture  the  change  of  state  as  accompanied 
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by  the  formation  of  numerous  small  cavities  along  the  gliding  surface,  so  that 
the  action  may  be  not  inconsistent  with  the  reduction  of  density  so  commonly 
recorded  in  experiment.  That  such  cavities  are  formed  on  the  gliding  surfaces 
may  be  inferred  also  from  many  phenomena  associated  with  mechanical  fatigue. 

It  is  highly  probable  that  the  work  done  on  the  small  masses  that  suffer  the 
change  is  drawn  in  part  from  the  strain-energy  of  surrounding  crystalline 
masses  that  suffer  no  change  of  state.  Were  this  not  so,  the  gliding  surfaces 
would  doubtless  spread  more  widely  than  they  appear  to  do.  The  smaller  the 
grain  size,  the  smaller  the  volume  from  which  a  given  change  can  draw  strain- 
energy,  and  the  greater  the  stresses  required  to  attain  the  elastic  limit. 

Few  metals  being  free  from  internal  strain,  quantities  of  work  must  be 
stored  in  individual  grains  before  external  stresses  are  applied.  These  quan- 
tities doubtless  contribute  to  the  total  constant  energy  required  to  effect  the 
change ;  but  at  present  they  lie  beyond  the  possibility  of  estimation. 

The  microstructure  of  most  ordinary  metals  being  highly  complex,  the  distri- 
butions of  strain-energy  between  individual  grains,  or  parts  of  grains,  must  be 
exceedingly  irregular,  particularly  when — as  in  the  case  of  Ferrite  and  Cementite 
— the  elastic  constants  differ  widely.  Even  in  pure  metals  and  solid  solutions, 
the  non-isotropic  properties  of  the  grains  must  result  in  irregular  distributions 
of  energy,  so  that  some  grains  suffer  permanent  strain  before  others.  As  a 
general  rule,  however,  we  may  assume  that  the  ratio  between  the  maximum  and 
mean  intensities  of  strain-energy  will  be  nearly  constant  unless  the  grain  size  is 
abnormal,  so  that  the  mean  value  may  be  expected  to  follow  the  law  governing 
the  maximum.  The  ratio  between  the  elastic  limit  and  the  yield-point  may  be 
closely  related  to  this  irregular  distribution  of  energy. 

In  view  of  these  considerations,  and  in  spite  of  the  numerous  uncertainties 
of  detail,  it  seemed  to  the  author  probable  that  the  strain-energy  per  unit 
volume  might  attain,  at  the  elastic  limit,  a  nearly  constant 'limiting  value  inde- 
pendent of  the  simple  or  complex  nature  of  the  applied  stresses ;  or  if  the 
strain-energy  were  found  to  vary  appreciably  and  definitely,  that  its  changes 
might  afford  a  serviceable  clue  to  the  investigation  of  the  more  complex 
phenomena  of  strain.  On  plotting  the  graphs  to  represent  the  hypothetical 
values  of  the  elastic  limit  given  by  the  equation 

(X2+Y2+ Z2)— 2ff  (YZ  +  ZX+XY)=F2 

where  F  signifies  the  elastic  limit  in  simple  tension,  it  was  found  that  the 
experimental  loci  were  very  close  to  the  hypothetical,  not  only  for  compara- 
tively purs  metals,  but  also  for  those  of  more  complex  structure,  such  as 
Pearlitic  steels.  Judging  from  the  results  obtained  in  a  variety  of  applications, 
the  above  relation  gives,  in  the  opinion  of  the  author,  a  very  reliable  measure 
of  the  relation  between  the  elastic  limits. 

It  is  not  suggested  that  the  non-ductile  fracture  of  a  brittle  material,  or 
even  the  ultimate  strength  of  a  ductile  metal,  is  governed  by  the  thermo- 
dynamic  considerations  that  have  been  set  forth.  Even  the  yield-points  are  only 
approximately  concerned,  and  the  stricter  application  is  limited  to  the  relation 
between  the  elastic  limits. 


VI. 

Alternating  Combined  Stress  Experiments. 
By  W.  MASON,  D.Sc.,  and  W.  J.  DELANEY,  B.Eng. 

In  1914-15  some  tests  on  the  dead  mild  steel  furnished  by  the  British  Associa- 
tion Committee  on  Complex  Stress  Distribution  were  published.1  These  tests 
included  alternating  torsion  and  alternating  bending,  but  these  straining  actions 
were  each  applied  separately.  The  experiments  now  described  have  been  made 
during  the  past  ten  months  upon  the  same  batch  of  steel  j  but  in  these  later 
tests  the  alternating  bending  and  torsion  were  simultaneously  applied.  The 

1  '*  Alternating  Stress  Experiments/  Proc.  Inst.  Mech.  Eng.,  January-May, 
1917. 
1921  A  A 
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main  feature  of  all  these  experiments  was  measurement  of  the  ranges  of  cyclic 
bending  and  twisting  strain. 

The  machine  used  was  the  one  described  in  the  paper  above  referred  to,  and 
the  frequency  of  the  cycles  of  stress  was,  as  previously,  200  per  minute.  The 
apparatus  for  measurement  of  the  strains  is  also  substantially  the  same,  and  the 
ranges  of  strain  simultaneously  measured  are  the  angle  of  twist,  and  the  angular 
deflection  due  to  bending  between  two  sections  about  4  in.  apart.  The 
bending  moment  was  uniform  over  this  length,  which  included  2^  in.  turned 
parallel  between  shoulders.  The  mirrors  for  measuring  these  strains  were  fixed 
to  the  squared  shoulders  of  the  specimen. 

All  the  conditions  were  the  same  as  for  the  experiments  previously  made  in 
bending  and  twisting  separately;  the  specimens  were  cut  from  the  same  two 
bars  A  and  B,  the  number  of  specimen  denoting  the. position  of  the  specimen 
in  the  uncut  bar.  These  conditions  allow  a  systematic  comparison  to  be  made 
between  these  tests  and  those  done  previously — a  procedure  of  obvious  advantage 
in  saving  of  time  and  labour,  especially  since  simultaneous  measurement  of  the 
two  sets  of  cyclic  strains  entails  considerable  patience. 

Phase  Angles  of  Bending  and  Twisting  in  Combined  Cycles. 

A  set  of  six  specimens  (Tables  I.  and  IV.)  was  tested,  four  of  which  were 
solid  and  two  hollow.  These  were  of  similar  shape  and  dimensions,  and  the 
surfaces  were  finished  in  the  same  manner  as  previously.  It  will  be  seen  that 
in  four  of  these  the  phase  angle  (i.e.  the  angle  between  the  crank  giving  the 
simple  harmonic  alternating  bending  moment  and  the  crank  which  applied  the 
S.H.  alternating  torque)  was  90°;  and  in  two  of  them  the  phase  angle  of  them 
was  zero.  In  the  case  of  Specimens  A6,  A13,  and  B31,  the  maximum  direct 
stress  calculated  from  the  bending  moment  was  approximately  twice  the 
maximum  shear  stress  calculated  from  the  torque.  In  calculating  the  stresses, 
the  usual  simple  formulae  for  bending  and  twisting  were  used. 

TABLE  IV. 

PRELIMINARY  DEAD-WEIGHT  TESTS  AND  DIAMETERS  OF  SPECIMENS. 


Specimen 

Preliminary  Dead-weight  Test 

k 

1 

0  q0^^ 

I&H 

=  11 

Si 

02   o     +] 

.5  ® 

1 

i 

**-*  CO      II 
O  ""^      1* 

j'gS 

SDJJ+I 

8>1^ 

Remarks 

1 

B 

|ij  n 

§>-2  S 

lli-2 

||  " 

fw    +1 

Ii« 

<B    * 

S3     || 

Sg 

T3 

^^ 

inch 

lb,in. 

tons  per 
sq.  in. 

cm.  on 
scale 

lb.-in. 

tons  per 
sq.  in. 

cm.  on 
scale 

/  '  Alternating 
See    pre-    1  stress  ^xperi- 

B31 

0-3125  solid 

±   54 

±  8-00 

5-70 

±   54 

±4-00 

7-90 

vioustest.'  ments/ 

P-  U1           Meek.  Eng. 

\  Jan.-May,l917 

A13 

0-4375  solid 

±183 

±10-00 

6-60 

±183 

±5-00 

7-20 

Accidentally  bent  and 
afterwards  straightened 

A6 

hollow,  see 

±172 

±10-00 

4-90 

±172 

±5-00 

5-60 



below 

B24 

hollow,  see 

±149 

±  8-54 

4-20 

±149 

±4-27 

4-90 



below 

A9 

0-4376  solid 

±183 

+10-00 

5-90 

±183 

+  5-00 

7-20 



B17 

0-4371  solid 

±179 

±  9-75 

5-95 

±179 

±4-87 

6-80 

*  1=  Moment  of  Inertia  of  Section  about  a  diameter,  (inch)4, 
•j-  J=  Polar  Moment  of  Inertia  of  Section,  (inch)4. 

Number 
of 
Specimen 

External 
Diameter 
(inch) 

Internal 
Diameter 
(inch) 

Thickness  of  Wall 

Max. 
(inch) 

Min. 
(inch) 

Mean 
(inch) 

Remarks 

A6 
B24 

0-5760 
0-5770 

0-5050 
0-5050 

0-0360 
0-0380 

0-0355 
0-0340 

0-0357 
0-0360 

Bore  eccentric 

COMPLEX  STRESS  DISTRIBUTIONS  IN  ENGINEERING  MATERIALS.  333 

With  phase  angle  90°,  the  stress  p  due  to  bending  and  the  shear  stress  q  due 
to  twisting  may  at  any  instant  t  be  represented  thus — 


p=pi  cos  Zir-         q=qi  cos  I  ^+*' 

where  pl  is  the  maximum  direct  stress  due  to  the  cycle  of  bending, 
q]      „  „  shear  „  „  „          twisting, 

and  T  is  the  period  of  either  cycle. 

The  principal  stresses  at  the  time  t  will  be — 

2  cos82irl+4g!2  sin3 

Now  if  pl  =  2q1=2a  (say),  as  in  the  three  tests  just  mentioned,  the  principal 
stresses  are 

acos2ir^  +  a 

and  the  maximum  shear  stress  due  to  the  continued  action  is  thus  always  equal 
to  '  a.' 

The  regions  at  which  the  maximum  shear  stress  '  a '  always  exists  are,  of . 
course,  at  the  diametrically  opposite  skins  of  the  specimen  in  the  plane  of  the 
bending  moment ;  the  planes  of  the  shear  '  a '  are  at  right  angles  to  the  skin 
and  rotate  uniformly,  half  a  revolution  made  in  the  time  T  of  one  cycle.  On 
other  planes  at  45°  to  a  principal  stress  (the  principal  planes  also  make  one  half 
revolution  in  time  T),  the  maximum  shear  stresses  of  the  cycle  vary  between 

'  a '  and  —  according  to  the  position  of  the  45°  plane  with  respect  to  the  plane 

of  the  principal  stresses.  Comparing  this  state  of  stress  with  the  stress-system 
induced  by  an  alternating  bending  or  an  alternating  twisting  separately  applied, 
the  number  of  planes  exposed  to  the  maximum  shear  stress  '  a '  is  indefinitely 
increased,  and  correspondingly  the  gliding  planes  of  a  much  larger  number  of 
crystals  are  subjected  to  alternating  shear  stress  of  maximum  value  'a.' 

In  the  tests  of  Nos.  A9  and  B17  the  phase  angle  was  zero.  The  range  of 
twisting  moment  on  A9  was  constant  throughout,  while  the  bending  moment 
was  increased  by  stages.  B17  was  subjected  to  constant  ranges  of  both  bending 
and  twisting  until  fracture  occurred.  Both  tests  were  continuous,  the  machine 
being  run  throughout  without  stopping. 

Discussion  of  Tables  II.  and  III. 

Tables  II.  and  III.  have  been  prepared  in  order  that  a  comparison  may  be 
made  conveniently  ^between  the  stresses  and  strains  (i)  in  the  90°  phase  com- 
bined stress  tests,  and  in  (ii)  tests  wherein  the  bending  and  twisting  were 
applied  separately.  Hollow  and  solid  specimens  are  put  in  separate  Tables  in 
view  of  a  previous  result  of  one  of  the  authors,  viz.  that  a  solid  specimen  under 
alternating  bending  or  twisting  is  apparently  considerably  stronger  than  a 
hollow  one.2  The  result  of  the  comparison  of  stresses  may  be  summarised 
thus — 

If  P  be  the  maximum  direct  stress  of  bending  cycles  separately  applied, 
if  Q  be  the  maximum  shear  stress  of  torsion  cycles  separately  applied, 
and  if  p1t  q1t  a,  have  the  meaning  previously  assigned, 

then  for  the  ranges  of  stress  which  produced  fracture,  it  was  found  that 
P>pi  and  Q>q\. 

The  maximum  shear  stress  due  to  bending  separately  applied  being  £P,  the 
relations  between  the  maximum  shear  stresses  are  therefore  |P  >  a  and  Q>a. 

It  should  be  observed  that  these  results  are  got  by  comparing  only  two 
specimens  tested  with  combined  alternating  stress  at  phase  angle  90°.  A  con- 
sideration of  B31  also  supports  these  conclusions  if  the  smaller  number  of  cycles 
endured  by  it  is  considered.  It  may  be  remarked,  also,  in  view  of  the  small 

2  'Alternating  Stress  Experiments,'  Proc.  Inst.  Mech.  Eng.,  January-May, 
1917,  p.  151.  The  stresses  are  calculated  (see  Table  I.)  from  formulas  which 
assume  that  stresses  are  proportional  to  strains,  which  gives  a  calculated 
maximum  stress  greater  than  the  actual,  more  especially  for  solid  specimens. 
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number  of  tests  made,  that  measurement  of  the  cyclic  strains  furnishes  a  check 
on  the  accuracy  of   the  applied  stresses   as  well    as  on  the   normality  of  the 
specimen.     The  columns  headed  '  Preliminary  Deadweight  Test '  in  Table  IV. 
are  inserted  so  that  this  and  other  comparisons  may  be  made. 
It  is  probable  that  the  differences  of  shear  stress 

|P— a          and          Q— a 

would  have  been  rather  larger  if  the  number  of  cycles  at  the  stresses  giving 
fracture  in  the  Combined  tests  had  been  as  large  as  the  number  of  cycles  in 
these  separate  tests.  It  would  appear,  then,  that  there  is  strong  prima-facie 
evidence  for  believing  that  for  mild  steel  in  combined  tests  with  phase  angle  90Q 
and  with  /p^=2q^  the  limiting  range  of  stress  is  less  than  the  limiting  ranges 
for  separate  tests  by  about  10  per  cent. 

With  regard  to  the  cyclic  strains  (phase  difference  90°)  it  will  be  observed 
that  the  percentage  increases  of  the  component  bending  and  twisting  strains 
measured  in  the  combined  tests  are  practically  the  same  as  those  measured  for 
the  same  stresses  in  the  separate  tests;  in  other  words,  the  total  (i.e.  elastic 
+  non-elastic)  strains  pertaining  to  the  bending  and  twisting  applied,  simul- 
taneously are  practically  the  same  as  if  cycles  of  the  same  ranges  of  stress  had 
been  applied  separately.  But  the  ranges  of  stress  causing  fracture  in  combined 
tests  are  about  10  per  cent,  less  than  the  ranges  that  would  have  produced 
fracture  had^the  bending  and  twisting  been  applied  separately;  thus  the  com- 
ponent ranges  of  strain  at  stress  ranges  approaching  fracture  are,  as  seen  in 
Tables  II.  and  III.,  much  less  in  combined  than  in  separate  tests.  Comparing 
solid  specimens,  the  percentage  decrease  of  total  strain  comes  out  to  be  36  per 
cent,  for  the  torsional  strain,  and  20  per  cent,  for  the  bending  strain ;  and  com- 
paring hollow  specimens,  the  respective  figures  are  approximately  the  same,  viz. 
34  and  24.  To  obtain  these  figures,  A6  and  A13  were  compared  respectively 
with  A14,  A16,  B22,  B27,3  and  A7,  All,  A15.3 

Remarks  on  Component  Strains  in  Tests  with  Phase  Angle  90°. 

An  interesting  point  observed  in  the  tests  with  phase  difference  90°  was  the 
following.  After  a  run  under  a  pair  of  simultaneous  stresses,  during  which 
there  was  cyclic  non-elastic  strain,  an  increase  of  one  of  the  pair  of  stresses, 
say  the  bending  stress,  caused  a  slight  contraction  in  the  range  of  torsional 
strain,  and  vice  versa.  The  test  of  Specimen  B24  was  specially  directed  towards 
the  study  of  this  effect.  To  obtain  an  enlarged  effect,  one  of  the  simultaneous 
stresses  was  made  zero  and  a  run  of  several  thousand  cycles  given.  Then, 
without  stopping  or  alteration  of  speed,  the  other  stress  was  reimposed. 

For  example,  in  course  of  test  of  B24, 

(a)  With   range  of  direct  stress  due  to  bending  =  ±  13'4  tons  per  sq.    in. 

and  range  of  torque  zero,  the  range  of  bending  strain  increased 
from  7'35  to  7'58  cm.4  during  40,000  cycles.  On  adding  torque  to 
give  +  5'30  tons  per  sq.  in.  of  shear  stress  and  retaining  the  same 
range  of  bending  moment,  the  bending  strain  decreased  to,  and 
remained  at,  7*33  cm.  during  12,000  cycles. 

(b)  Again,  with   range    of  shear    stress   due   to   torsion  = +6'70   tons/p" 

and  range  of  bending  moment  zero,  the  range  of  torsional  strain 
increased  from  10'9  to  13*5  cm.  during  a  run  of  108,000  cycles. 
On  adding  bending  moment  to  give  +11'95  tons  per  sq.  in.  and 
retaining  the  same  torque,  the  range  of  torsional  strain  decreased  to, 
and  remained  at,  11-80  cm.  during  25,000  cycles. 

This  effect  is  not  due  to  friction  occasioned  by  the  extra  loading  on  the 
parts  of  the  machine  when  the  additional  stress  is  put  on.  In  case  (b)  the  effect 
of  friction  would  be  the  reverse  of  that  observed ;  and  while  in  (a)  the  effect  of 
friction  would  be  similar  to  that  observed,  the  magnitude  of  the  friction  (of 
ball  bearings)  would  not  be  sufficient  to  account  for  the  decrease.  Moreover, 

3  See  'Alternating  Stress  Experiments'  (Proc.  Inst.  Mech.  Eng.,  January- 
May,  1917)  for  record  of  these  tests. 

4  On  image  of  scale  reflected  from  mirrors. 
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with  phase  difference  90°,  the  maximum  of  bending  moment  coincides  with  zero 
torque,  and  vice  versa-,  so  that  friction  must  have  diminished  to  zero  at  the 
time  when  the  effect  is  measured. 

But  another  circumstance  must  be  considered.  Suppose  the  specimen  to  be 
simultaneously  stressed,  and  consider  the  instant  when  the  torque  is  a  maximum 
and  the  bending  moment  zero.  There  still  remains  the  bending  due  to 
hysteresis ;  and  application  of  the  torque  to  the  bent  specimen  would  give  a 
reduqed  torsional  stressing  and  straining  at  the  parts  deflected  away  from  the 
axis  of  the  specimen  by  this  hysteresis  bending.  But  a  rough  calculation  shows 
that  the  observed  effect  is  too  large  to  admit  of  this  explanation,  the  hysteresis 
bending  (though  this  was  not  measured)  being  evidently  too  small.  The  authors 
cannot  venture  to  attempt  an  explanation  of  this  phenomenon.  Since  the  ranges 
of  strains  in  combined  tests  agree  substantially  in  magnitude  with  those  under 
the  same  ranges  of  stress  applied  separately,  it  might  appear  that  the  work 
done  in  a  combined  cycle  would  be  equal  to  the  sum  of  that  done  'in  the  same 
cycles  applied  separately.  But  this  cannot  be  assumed,  since  the  width  and 
shape  of  the  corresponding  hysteresis  curves  may  not  be  similar.  The  authors 
have  not  as  yet  attempted  to  measure  the  width  of  the  hysteresis  diagram ;  and 
until  this  is  done  it  seems  too  early  to  attempt  any  explanation. 

Rernar/cs  on  Component  Strains  with  Phase  Angle  Zero. 

Turning  now  to  the  two  tests  with  the  bending  and  twisting  in  phase  with 
each  other — in  B17  the  range  of  direct  stress  due  to  bending  was  twice  the  range 
of  shear  stress  due  to  twisting ;  while  in  A9  the  ratio  between  these  respectively 
was  increased  during  the  test,  and  during  the  last  stage,  at  the  end  of  which 
fracture  occurred,  the  ratio  was  approximately  2'5.  Comparison  of  the  ranges 
of  shear  stress  with  those  of  Al,  A7,  and  All,  A15,  shows  that  the  average 
values  of  the  induced  maximum  shear  stresses  producing  fracture  are  almost 
exactly  the  same  for  both  ^multaneous  and  for  separate  applications  of  bending 
and  twisting.  The  material  in  this  as  in  other  respects  already  found,  con- 
forms approximately  to  an  extension  of  Guest's  Law  to  alternating  stresses.5  6 
These  shear  stresses  are  thus  about  10  per  cent,  greater  than  that  of  the  average 
of  the  shear  producing  fracture  when  the  phase  difference  is  90°  and  p^=2qr 

The  component  cyclic  strains  (elastic  +  non-elastic)  induced  towards  the  end 
of  the  last  stage  of  the  combined  tests  (phase  angle  zero)  were  much  below  those 
in  the  final  stages  of  comparable  tests  in  separate  bending  and  twisting.  Thus, 
comparing  A9  with  A7,  the  reduction  in  range  of  total  torsional  (including 
both  elastic  and  non-elastic)  cyclic  strain  is  50  per  cent.  ;  and  comparing  A9 
with  All  and  A15,  the  reduction  in  total  cyclic  bending  strain  is  23  per  cent. 
Comparing  A17  with  A7,  the  reduction  in  range  of  total  torsional  strain  is 
47  per  cent. ;  and,  again,  comparing  A17  with  All  and  A15,  the  reduction  in  range 
of  total  cyclic  bending  strain  is  36  per  cent.  It  should  be  pointed  out,  however, 
that  the  components  of  the  alternating  combined  stresses  at  fracture  range  were 
individually  much  less  than  the  stresses  at  fracture  range  due  to  alternating 
bending  and  twisting  applied  separately. 

Principal  Strains. 

Although  the  component  strains  at  fracture  ranges  in  the  combined  tests 
are  less  than  the  respective  strains  in  bending  or  torsion  separately 'applied  to 
different  specimens,  it  does  not  follow  that  the  maximum  principal  strains  will 
be  correspondingly  smaller.  It  is  therefore  of  interest  to  make  comparison  of 
the  maximum  principal  strains.  Having  calculated  the  latter,  it  is  easy  to 
obtain 7  the  greatest  strain  difference.  These  quantities  have  an  additional 
interest  inasmuch  as  strains  may  be  assumed  to  follow  a  linear  law  of  distribution 
from  axis  to  skin  of  specimen  even  when  the  strains  are  no  longer  elastic, 
whereas  in  the  latter  condition  the  stresses  have  not  a  linear  distribution  and 
are  not  accurately  known. 

5  See  'Alternating  Stress  Experiments,'  Proc.   Inst.  Mech.  Eng.,  January- 
May,  1917,  p.  149. 

6  B.  P.  Haigh.     Proc.  Inst.  Mech.  Eng.,  January-May,  1917,  p.   100. 

7  As  suggested  by  Prof.  L.  N.  G.  Filon. 
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The  distances  of  scale  from  indicating  mirrors  are  the  same  in  the  present 
tests  as  in  those  made  four  years  ago,  and  are — 

For  bending  strain  apparatus        ....     165  cm. 
,,     torsional      „  ,,  ....     155    ,, 

For  bending,  the  mirrors  registered  a  scale  displacement  proportional  to  the 
angle  between  the  tangents  at  the  ends  of  the  specimen.  The  actual  value  of 
this  angle  is 

*  Range  of  Strain'  in  cm._Rb    (say) 
2x2x165  ~660 

As  mentioned  on  page  331,  the  mirrors  were  fixed  to  the  square  shoulders  of 
the  specimen,  and  the  length  I  over  which  this  angle  is  measured  is  not  known 
accurately  If  I  be  regarded  as  an  equivalent  length 8  of  specimen,  and  if 
r  stands  for  the  (outside)  radius  of  specimen,  then  the  maximum  fibre  strain 
will  (provided  the  strains  have  a  linear  distribution  and  are  proportional  to 
distance  from  axis  of  specimen)  be 


-=  (angle  between  tangents)  =-    — - 


The  angle  of  twist  over  length  I  will  be 

*  Range  of  Strain'  in  cm._RT  (say) 
2  x  2  x  155  ~620 

and  the  shear  strain  will  be  accordingly    -  .  —~ 

In   the    Preliminary   Dead-weight    Tests    (see   Table    IV.),    in    which   the 
material  is  elastic,  the  bending  strain 

=  ^  «?*=^  (1) 

I    660      E 

where  +pe  is  the  range  of  stress  which  gives  the  range  of  elastic  strain  eRb  (em.), 
and  E  =  Young's  Modulus. 

The  shear  stress  in  torsion  similarly 


_ 
I   620      C 


where   +  qe    is   the   range   of    shear    stress   giving   the    range   of    elastic   strain 
eRT    (cm.),  and  C  is  the  modulus  of  rigidity. 

If  we  consider  a  specimen  on  which  preliminary  dead-weight  tests  were  made 
both  in  bending  and  torsion,  and  if  we  assume  that  the  fillet  at  the  shoulders 
affected  the  bending  and  twisting  stresses  and  strains  in  the  same  proportion, 
then  — 


Pe  *&?  660 
qe  eRb  620 


8  Equivalent  in  the  sense  that  the  angle  of  bending  over  a  parallel-turned 
length  I  is  to  be  the  same  as  that  measured  on  the  length  between  the  square 
shoulders. 
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In  the  '  preliminary  dead-  weight  tests  '  of  A6,  B24,  A9,  BIT,    2±  =  2,    and 

E  » 

we  obtain  the  following  values  of  ~-:  — 


For  B24=2xx       =  2-48 
For  A9       =  2x        X  §    =  2-59 


giving  a  mean  value  for    —   of  2-48>  which,  incidentally,  makes  Poisson's  Ratio 

C 
=  0.24. 

Assuming   E  to  be  30xl06   Ib.   per  sq.    in.,  this  ratio    of   E  to   C    makes 
C  =  12.1xl06. 

E 

For  specimen  B24,  which  happens  to  have  the  mean  of  the  above  values  of  ^, 

r  .bending       660     8'54x  2240     n.inn 
I  (t^tj  =  4=2  X     30X10*    = 

Taking  any  test   in  bending,  let  R   (cm.)   be  any    range  of   strain,   elastic  or 
otherwise,  then  strain  corresponding  to  R 

=r  J* 

I  660' 

and   substituting  the  value  of—   for  the  specimen  from  equation  (1),  we  have 
this  strain 


For  any  test  in  Torsion,  the  greater  principal  strain  corresponding  to  a  range 
of  strain  R  (cm.) 

=r     R    =9e  620        R 

~l  1240  ~  C   A  '  1240 

=  Rif  X9-27X10-5  .    (4) 

eRT 

The  strains  given  by  formulae  (3)  and  (4)  will  be  exact  if  the  following  assump- 
tions are  valid  :  — 

(1)  that  E  and  C  are  the  same  for  the  material  of  all  the   specimens,  and 
that  E  =  30xl06lb./n"; 

(2)  that  the  fillets  connecting  the  parallel  part  with  shoulders  of  a  specimen 
have,   in    that   specimen,   the   same  proportional    effect   on   both   bending   and 
twisting  strains  and  stresses  ; 

(3)  that  the  strains  have  a  linear  distribution  from  axis  outwards. 

Using  formulae  (3)  and  \4>),  the  greater  principal  strains  and  the  greatest 
strain  difference  have  been  worked  out,  and  set  forth  in  Table  V.,  for  the 
separately  applied  bending  and  twisting  tests  previously  done. 
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TABLE   V. 
BENDING  OR  TORSION  SEPARATELY  APPLIED. 


Specimen 

Kind  of  Test 

Eange  of 
'  Strain.' 
Cms.  on 
Scale 
=  R 

Data  from  • 
Dead-weight  Test 

Maximum 
Greater 
Principal 
Strain 

Max!  mum 
of 
Greatest 
Strain 
Difference 

Maxi- 
mum 
Shear 
Stress 

Pe°I<Je 

.VA 

tons/D" 

cms.  on 
scale 

A15  solid 

Bending  only 

10-33 

10-00 

6-12 

•00126 

•00157 

6-87 

All       „ 

55 

11-97 

10-50 

6-60 

•00141 

•00176 

7-15 

A14  hollow 

9-87 

10-00 

5-36 

•00136 

•00171 

6-37 

A16       „ 

» 

8-18 

10-00 

5-13 

•00118 

•00147 

6-37 

A7  solid 

Torsion  only 

14-67 

5-00 

7-01 

•000975 

•00195 

7-65 

A10  hollow 

*  i           * 

9-75 

4-75 

5-27 

•000813 

•00163 

6-50 

B22       „ 

10-72 

5-00 

5-86 

•00085(5 

•00170 

6-25 

B27       „ 

11-36 

5-00 

6-03 

•000877 

•00175 

6-10 

B28       „ 

»>            " 

10-30 

'5-00 

5-93 

•000805 

•00161 

6-25 

B30       „ 

»»            » 

10-65 

5-00 

5-87 

•000842 

•00168 

6-25 

NOTE. — The  data  in  the  first  five  columns  are  taken  from  Tables  2  and  3,  pp.  136-146,  Proc.  Inst. 
Mech.  Eng.,  January-May,  1917. 

It  will  be  observed  that  the  values  of  the  greater  principal  strains  are  widely 
divergent.  The  values  of  the  greatest  strain  difference  are  not  so  uniform 
as  the  maximum  shear  stresses,  especially  if  account  be  taken  of  the  enhanced 
values  of  the  latter  for  solid  specimens,  which  greater  values  are  presumably  due 
to  the  formula  of  calculation. 

It  should  be  observed  that  the  greatest  strain  differences  for^  bending,  in 
Table  V.,  are  calculated  on  the  assumption  that  the  minor  principal  strains  ez 
are  equal  to  <re\  where  <r  is  Poisson's  Ratio  and  ei  is  the  (greater  principal) 
strain.  This  would  be  exact  if  the  stresses  were  elastic,  but  since  there  is 
a  non-elastic"  element,  the  results  for  e2  are  approximate.  The  non-elastic  ele- 
ment is,  however,  only  about  one-sixth  (at  most)  of  the  total  strain,  so  that  the 
error  in  assuming  that  the  greatest  strain  difference  (bending  tests  only), 

viz.  61 — e2,  is  equal  to  e\  ~f~T    makes,  at  the  most,  an  error  of  4  per  cent,  in  the 

greatest  strain  difference. 

To  apply  formulae  (3)  and  (4)  to  calculate  the  strain  components  for  the 
combined  tests,  the  assumption  is  made  that  these  components  are  linearly 
distributed  as  in  separately  applied  testing.  The  effect  described  on  pp.  334 
and  335  with  respect  to  decrease  of  one  component  strain  when  the  other  is 
increased  may  possibly  be  due  to  non-linear  distribution  of  these  strains.  It  is 
probable,  however,  that  there  will  be  no  marked  departure  from  linearity,  and 
formulae  (3)  and  (4)  have  been  used  to  evaluate  the  component  strains  for 
fracture  ranges. 

Given  two  direct  strains  e\  and  e2  at  right  angles  to  each  other,  and  a  slide 
<J>  on  planes  parallel  to  BI  and  ez — which  are  the  circumstances  corresponding 
to  the  localities  of  greatest  stress  in  the  specimens — then,  if  8j  and  52  are  the 
principal  strains, 

\ 


and  the  greatest  strain  difference  is 

Si-52=\/(ei 
where 


and 


-X  7-43  X  KM 

etib 

=— <rR-^-  X7-43X10"5 


(5) 
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TABLE   VI. 

COMBINED*  TESTS. 

NOTE. — For  tests  at  90°  phase  angle  the  greatest  values  of  the  strain  components  and  of  the  Greater 
Principal  Strains  and  Greatest  Strain  Difference  are  given  both  for  (a)  Maximum  Bending  and  Zero 
Torque,  and  (6)  Maximum  Torque  and  Zero  Bending 


Width  of 

' 

a 

3b 

If 

Jl  o 

Hysteresis 

Strain  Components 

Greater 

l§ 

alt" 

M 

Specime  i 

'o'a 
<D  <u 

.21 

Cm.  on  Scale 

Principal 
Strain 

ill 

1° 

a 

§  § 

JH 

Bend- 

Tor- 

XlO» 

9  <2 
"P  £ 

« 

rt  ° 

ing 

sion 

«1X103 

4.X  103 

o 

|S 

(a)  0-92 

(a)  0-06 

(a)  0-92 

(a)  1-15 

/A6 

6-05 

6-75 

0-42 

0-36 

•"  

5-75 

hollow 

(6)  0-064 

(6)  1-12 

(6)  0-58 

(6)  1-12 

fe 

o=> 

(a)  1-115 

(a)  negligible 

(a)  1-115 

(a)  1-39 

£  S3 

^  8 

&  *«  J 

A13 

8-85 

9-35 

1-22 

0-17 



^™" 

6-48 

is 

si 

solid 

(6)  0-154 
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The  values  of  8j  and  5,— S2  taken  in  Table  VI.  have  been  calculated  from  formulas  (5)r 
Since  <r  has  been  taken  as  — ,  we  have 


5 

i— 62= - 


and  formulae  5  reduce  to 


and 


'-8'=A/ff 


In  working  out  the  values  of  5j  and  Sl— 82  for  the  specimens  tested  with 
phase  angle  90°,  it  has  been  kept  in  mind  that  the  maximum  strain  of  bending 
will  coincide  with  only  the  residual  (or  hysteresis)  strain  at  zero  torque,  and 
vice  versa.  These  hysteresis  strains  could  not  be  measured,  but  the  width  of 
hysteresis  loop  may  be  taken  (since  the  non-elastic  element  is,  at  most,  only 
about  one-sixth  of  the  total  strain)  as  the  difference  between  the  range  of  strain 
measured  and  the  range  that  would  have  existed  if  the  material  had  remained 
elastic.  Considering  two  loops,  one  for  bending  and  the  other,  differing  in 
phase  by  90°  from  the  first,  for  the  torsion,  it  is  possible  to  correlate  exactly 
only  four  points  of  the  one  loop  with  four  points  of  the  other.  These  four 
points  are,  of  course,  the  maximum  of  strain  of  one  kind  with  the  strain 
at  zero  stress  of  the  other.  If  the  strains  were  elastic,  and  therefore  simple- 
harmonic  like  the  stresses,  any  one  point  of  the  one  stress-strain  curve  could 
of  course  be  correlated  with  its  corresponding  point  on  the  other  curve.  To 
find  the  position  of  crank,  which  gives  a  maximum  for  the  greatest  strain  differ- 
ence, a  simple  calculation  was  made  on  the  assumption  that  the  strains  were 
simple-harmonic.  On  this  assumption  the  maximum  for  this  quantity  comes 
out  to  be  at  either  maximum  bending  or  maximum  torque.  Since  the  residual 
(hysteresis)  strains  are  small  (see  Table  VI.),  the  one  or  the  other  of  these  two 
epochs  has  been  assumed  to  be  that  of  the  maximum  of  greatest  strain-difference. 
The  values  of  this  quantity  have  been  calculated  for  both  these  epochs,  and 
are  given  in  Table  VI. 
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In  Specimens  A6  and  A13  the  maximum  value  (underlined  in  Table  VI.)  of 
greatest  stress-difference  was  found  to  correspond  to  the  maximum  of  bending; 
and  in  B31,  to  correspond  to  the  maximum  of  torque.  It  is  clear  that  the 
maximum  value  of  .the  greater  principal  stress  in  A6  and  A13  must  also  corre- 
spond to  maximum  bending,  because  to  obtain  this  principal  strain  we  have  to 

take  the   numerical  sum  of     76!  and  A/  — — +^>2,   and  the  latter  of  these   two 

"   16ei2 

3 

expressions,  as  well  as  the    -  el    happens  to  be  a  maximum  at  the  maximum  of 

bending  moment. 

For  B31,  having  regard  to  the  relative  values  of  et  and  <p,  it  appears  that 
the  maximum  value  of  the  greater  principal  strain  is  also  at  the  maximum  of 
bending. 


delation  between  Maximum  Shear  Stresses  and  Greatest  Strain-difference. 

An  examination  of  Tables  V.  and  VI.  shows  that  the  values  of  the  maximum 
principal  strains  are  not  nearly  so  uniform  as  the  maximum  shear  stresses; 
but  they  are  not  so  widely  divergent  as  the  principal  stresses  (not  tabulated). 
As  might  be  expected,  the  values  of  the  greatest  strain- difference  are  rather 
more  uniform,  those  for  the  combined  tests  appearing  rather  smaller  than  those 
for  the  separate  tests. 

Fig.  14  has  been  plotted  from  Tables  V.  and  VI.  It  will  be  noted  that  the 
solid  specimens  group  themselves  to  the  right  for  the  reason  that  their  calculated 
stresses  are  too  large — i.e.  more  than  the  real  stresses.  Making  use  of  a  former 
result  of  one  of  the  authors,  tbat  the  stresses  calculated  for  solid  specimens 
(using  formulas  true  only  for  perfectly  elastic  conditions)  are  about  15  per 
cent,  too  high,  the  stresses  of  the  solid  specimens  have  been  reduced  in  the 
ratio  ±$%,  and  the  reduced  values  have  been  plotted  against  the  greatest  strain 
difference  in  fig.  15.  In  fig.  15  the  points  corresponding  to  the  hollow  specimens 
appear  in  the  same  positions  with  respect  to  the  axes  as  in  fig.  14.  Fig.  15 
illustrates  also  the  relation  between  the  magnitudes  of  the  maximum  shear 
stresses  of  the  differently-tested  specimens. 


VII. 
On  Some  Problems  Relating  to  the  Design  of  High-Speed  Discs. 

By  R.  V.  SOUTHWELL,  of  the  National  Physical  Laboratory. 

Recent  investigations  by  Prof.  H.  Lamb  and  by  the  present  author  l  appeared 
to  have  some  bearing  on  the  practical  problem  of  vibrations  in  turbine  discs,2 
and  the  work  has  accordingly  been  continued  at  the  National  Physical  Laboratory, 
with  the  assistance  of  Miss  B.  S.  Gough.  The  following  summary  of  the  progress 
made  has  been  written  in  the  hope  that  discussion  will  reveal  the  directions  in 
which  further  extension  could  be  most  usefully  directed.  The  accuracy  of  the 
calculations  depends  in  every  instance  upon  the  validity  of  the  theory  of  thin 
plates,  as  applied  to  the  problems  treated ;  and  since  the  limitations  of  this 
theory  are  not  as  yet  fully  understood,  comparative  experiments  will  be  of  the 
greatest  value. 

1  'The  Vibrations  of  a  Spinning  Disc,'  Proc.  Roy.  Soc.  (A),  vol.  99  (1921), 
pp.  272-280. 

2  Cf.  K.  Baumann,   '  Some  Recent   Developments  in  Large  Steam  Turbine 
Practice,'  Jour.  Inst.  Elect.   Eng.  (1921). 
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The  following  problems  have  been  considered,  in  the  belief  that  they  are  the 
most  important  which  confront  the  designer  : — 

I.  The  calculation  of  Centrifugal  Stresses. 

II.  The  calculation  of  Transverse  Deflections  in  Diaphragms  and  Discs. 
UI.  The  calculation  of  Critical  Speeds  (or  of  the  Normal  Frequencies  of 
Free  Transverse  Vibration)  in  Diaphragms  and  Discs. 


FIG.  16. 

Two  types  of  disc  have  been  considered,  as  shown  in  the  above  figure  :  type  A 
is  a  disc  of  uniform  thickness,  and  type  B  a  disc  of  which  the  thickness  varies 
as  some  (negative)  power  of  the  radial  distance.  The  fact  that  theory  demands 
an  infinite  thickness  at  the  axis  is  not  a  practical  objection,  since  an  actual  disc 
is  bored  to  receive  a  shaft  of  finite  diameter,  and  allowance  for  this  circumstance 
can  be  made  by  means  of  an  auxiliary  stress -system. 


Problem  I.  Calculation  of  Centrifugal  Stresses  in  a  Disc  of  Specified  Profile. 

The  analysis  of  these  stresses,  in  a  disc  of  uniform  thickness,  has  been  worked 
out  by  C.  Chree,3  and  a  discussion  of  the  problem,  with  an  extension  to  discs 
of  type  B,  is  contained  in  A.  Morley's  Strength  of  Materials.*  Recent  investi- 
gations indicate  that  the  centrifugal  stresses  in  a  disc  of  any  specified  profile 
can  be  determined  without  much  difficulty  by  means  of  graphical  methods,  but 
it  is  doubtful  whether  these  methods  are  urgently  needed  by  the  designer,  who 
will  probably  be  content  with  analysis  which  is  applicable  to  discs  of  type  B. 

Strictly  speaking,  the  loading  due  to  the  blades  will  riot  be  distributed 
uniformly  along  the  rim,  but  will  vary  from  point  to  point.  If  the  stress-distri- 
bution near  the  roots  of  the  blades  is  a  matter  of  practical  importance,  it  will 
be  desirable  to  investigate  this  auxiliary  stress  system  :  the  requisite  analysis 
should  not  prove  difficult,  for  since  the  blades  are  uniformly  spaced,  we  may 
conveniently  analyse  the  load  system,  as  regards  its  variation  along  the  circum- 
ference, into  Fourier  components,  and  of  these  it  is  probable  that  only  two  or 
three  will  be  important. 

3  Proc.  Comb.  Phil.  Soc.,  vol.  vii.   (1891),  part  iv. 

4  Chap.  xi. 
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Problem  II.  Calculation  of  the  Transverse   Deflections  Induced  in  a  Disc  by 
a  Specified  iboad  System. 

Two  cases  of  this  problem  arise:  — 

(a)  when  the  disc  is  stationary; 
(6)  when  the  disc  is  rotating. 

It  is  evident  that  rotation  must  increase  the  effective  flexural  rigidity  of  the 
disc,  since  work  has  to  be  done  against  the  centrifugal  stresses  by  any  agency 
which  produces  transverse. deflections. 

Under  the  heading  (a),  we  have  exact  solutions  for  the  deflection  of  a  disc  of 
uniform  thickness  (type  A),  under  various  systems  of  transverse  load  which  are 
symmetrical  about  its  axis.5  Recent  investigation  has  shown  that  the  analysis 
can  be  extended  without  difficulty  to  any  system  of  transverse  loading  which 
varies  in  intensity  as  some  power  of  the  radial  distance,  and  that  a  further 
extension  may  be  made  to  discs  of  type  B.  Further,  in  this  problem  also  it 
appears  that  graphical  methods  can  be  employed  to  find  the  deflection  due  to  any 
load  system  (axial  symmetry  only  being  postulated),  in  a  disc  of  any  specified 
profile  :  but  here  again  it  seems  doubtful  whether  any  real  demand  exists  for 
methods  of  such  generality. 

Under  the  heading  (b),  we  have  to  consider  the  effect  of  centrifugal  force, 
and  the  difficulty  confronts  us  that  the  form  of  the  deflected  disc,  as  well  as  the 
absolute  magnitude  of  the  deflections,  will  in  general  be  altered  by  the  rotation  : 
had  a  change  in  magnitude  been  the  only  effect,  we  could  very  easily  have  deter- 
mined its  amount.  But  since  the  centrifugal  forces  will  tend  to  reduce  the 
deflections,  the  exact  magnitude  of  their  effect  is  of  less  practical  interest  than 
would  have  been  the  case  if  it  had  increased  the  danger  of  rubbing  between 
fixed  and  moving  parts,  and  in  view  of  the  analytical  difficulties  associated  with 
an  exact  solution,  it  will  probably  be  sufficiently  accurate  to  start  with  the  exact 
solution  for  a  non-rotating  disc;  to  calculate  the  additional  transverse  loading 
which  would  be  required  to  maintain  this  deflection  against  the  centrifugal 
stresses  corresponding  to  the  specified  speed  of  rotation;  and  to  estimate  a 
mean  coefficient  whereby  this  may  be  represented  as  a  fractional  addition  to  the 
original  loading  :  we  shall  then  have  an  idea  of  the  centrifugal  effect,  expressed 
roughly  in  the  form  of  a  fractional  addition  to  the  nexural  rigidity. 

Problem  III.     Calculation   of   Critical   Speeds   (or  of  the  Normal  Frequencies 
of  Free  Transverse  Vibration)  in  a  Disc  of  Specified  Profile,  Stationary  or 

Rotating. 

There  is  evidence  that  resonance  is  liable  to  occur  between  the  frequencies  of 
free  vibration  which  are  natural  to  the  turbine  disc  and  the  frequencies  of  the 
small  periodic  disturbing  forces  which  come  into  play  when  the  turbine  is 
running.  The  designer's  problem  is  to  arrange  that  such  resonance  shall  not 
occur,  at  all  events  within  the  practical  speed  range,  and  evidently  there  are, 
theoretically  speaking,  two  ways  of  effecting  this  result  :  either  he  must  eliminate 
the  periodic  disturbing  forces,  or  he  must  arrange  that  their  frequencies  shall  be 
less  than  any  of  the  frequencies  natural  to  the  disc.  In  practice,  the  disturbing 
forces  are  not  always  avoidable,  and  their  origin  is  sometimes  obscure  :  but 
their  nature  can  best  be  guessed  from  the  nature  and  frequency  of  the  vibration 
which  they  excite,  and  thus  our  most  important  problem  is  the  determination  of 
the  natural  modes  and  frequencies  of  free  transverse  vibration,  for  a  disc  of  any 
specified  form. 

As  in  problem  II,  a  distinction  must  be  drawn  between  the  effects  of  the 
flexural  rigidity  and  of  the  effective  rigidity  induced  by  rotation,  and  the  neces- 
sity for  taking  rotation  into  account  is  much  greater  in  the  present  instance  :  but 
it  fortunately  happens  that  the  requisite  allowance  for  it  can  be  made  much 
more  easily.  From  experience  gained  in  similar  problems  we  may  expect  that  the 
nature  of  the  vibrations  will  be  different  according  as  the  disc  is  entirely  free 
or  to  some  extent  constrained,  and  that  the  former  conditions  will  yield  the 
simpler  solutions. 

5  Of.  A.  E.  H.  Love,  Theory  of  Elasticity,  §314  (a);  and  A.  Morley, 
Strength  of  Materials,  chap  xiii. 
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The  free  transverse  vibrations  of  a  disc  of  uniform  thickness,  entirely  uncon- 
strained, were  investigated  by  Kirchhoff  in  1859. 6  The  normal  modes  are 
characterised  by  concentric  nodal  circles  and  by  equally  spaced  nodal  diameters  : 
for  any  given  number  (s)  of  nodal  diameters,  the  frequency  rises  with  the 
number  (n)  of  the  nodal  circles,  and  in  regard  to  the  higher  frequencies  it  is 
found  that  the  effect  of  increasing  n  by  unity  is  approximately  the  same  as  that 
of  increasing  s  by  two.  The  frequency  of  the  free  vibration  depends  to  a  slight 
extent  upon  the  value  of  Poisson's  Ratio,  of  which  a  representative  value  for 
steel  is  0.3. 

Kirchhoff's  solution  fails  in  respect  of  those  types  of  vibration  in  which  the 
number  of  nodal  diameters  is  unity  or  zero,  when  the  disc  is  attached  at  its 
centre  to  a  shaft  and  rotor  of  inertia  sufficient  to  prevent  any  change  in  the 
position  or  slope  of  the  disc  at  that  point.  The.  necessary  modifications  have 
now  been  investigated,  and  the  graver  frequencies  of  vibration  calculated  for  all 
types  of  vibration  possible  to  a  non-rotating  disc.7  It  is  found  that  substantially 
exact  values  could  have  been  obtained  with  much  simpler  analysis,  by  using  a 
method  due  to  the  late  Lord  Rayleigh,  in  which  the  type  of  the  deflection 
occurring  in  the  vibration  is  assumed,  and  an  estimate  of  the  frequency  derived 
from  the  corresponding  expressions  for  the  potential  and  kinetic  energies.  The 
method  is  particularly  valuable,  in  that  it  may  be  applied  with  equal  accuracy 
(assuming  the  validity  of  the  theory  of  thin  plates)  to  discs  of  curved  profile  :  we 
are  thus  in  a  position  to  calculate  the  graver  frequencies  of  vibration  for  any 
given  disc. 

The  effects  of  rotation,  in  a  disc  of  uniform  thickness,  were  discussed  by 
Prof.  Lamb  and  the  present  author  in  the  paper  to  which  reference  has  been 
made  above.  When  the  disc  is  stressed  by  centrifugal  forces,  its  effective 
flexural  rigidity  is,  as  we  have  seen,  increased  :  even  though  it  were  practically 
a  membrane,  without  any  flexural  rigidity  of  its  own,  it  would  vibrate  with 
finite  periods  under  the  restoring  effect  of  the  centrifugal  system.  The  problem 
first  considered,  therefore,  was  that  of  finding  the  frequencies  of  free  transverse 
vibration  when  the  centrifugal  system  acts  alone,  and  this  was  solved  without 
difficulty  in  the  case  of  a  uniform  disc.  The  modes  are  generally  similar  in  their 
nature  to  those  of  the  non-rotating  disc,  and  the  frequency,  in  any  given  mode, 
varies  directly  as  the  speed  of  rotation.  No  alteration  in  frequency  is  entailed 
by  constraints  acting  at  the  centre. 

In  the  general  case,  where  vibration  occurs  under  a  restoring  system  to 
which  both  the  flexural  rigidity  and  the  centrifugal  stresses  contribute,  an  exact 
solution  would  be  very  difficult  to  obtain,  principally  for  the  reason  that  the 
type  of  the  vibration  will  itself  change  with  «.  But  the  need  for  an  exact 
solution  is  obviated  by  Lord  Rayleigh's  theorem,  that  a  small  error  in  the 
assumed  type  of  vibration,  when  the  frequency  is  estimated  by  his  method,  leads 
to  an  error  in  the  frequency  which  will  be  of  the  second  order ;  in  fact,  there  is 
no  need  to  do  any  further  calculation  whatever,  since  it  can  be  shown  that  the 
gravest  frequency  corresponding  to  any  definite  number  of  nodal  diameters  will 
be  given,  with  close  approximation,  by  an  expression  of  the  form 

P*=Piz+P22, 

where  /pl    is  the  frequency  found  by  neglecting   the  rotation,    and  p2   is  the 
frequency  found  by  neglecting  the  flexural  rigidity. 

It  has  been  shown,  further,  that  the  gravest  frequency  will  be  underestimated 
by  the  formula  just  given.  This  result  suggests  a  very  simple  method  for  investi- 
gating the  gravest  frequencies  natural  to  a  disc  of  curved  profile  :  for  the  method 
of  Lord  Rayleigh  may  be  applied  with  equal  accuracy  to  the  calculation  both  of 
p,z  and  of  p^2,  and  in  each  case,  if  the  gravest  frequency  is  under  investigation, 
it  will  give  figures  for  these  quantities  which  are  either  exact  or  over-estimated. 
Thus  the  two  errors  involved  in  these  approximate  methods  tend  to  cancel  one 
another  as  regards  their  effect  on  the  estimated  value  of  pz,  and  a  closely 
approximate  result  may  be  expected. 

6  '  Ueber    die    Schwingungen    einer    elastischen    Scheibe,'    Crelle's  Journal, 
vol.  40  (1850),  and  Pogg.  Ann.,  vol.  81  (1850).     The  essentials  of  the  analysis  are 
reproduced  in  Lord  Rayleigh's  Theory  of  Sound,  §§218,  219. 

7  It  is  these   graver  frequencies  which   are  of    primary    importance,    since 
vibrations  of  high  frequency  will  in  practice  be  eliminated  by  air  damping,  etc. 
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VIII. 

On  the  Stability  of  a  Rotating  Shaft,  Subjected  Simultaneously  to 

End  Thrust  and  Twist. 

By  R.  V.  SOUTHWELL  and  BARBARA  S.  GOUGH,  of  the 
National  Physical  Laboratory. 

§1.  Nearly  forty  years  ago,  Sir  George  Greenhill  worked  out  the  criterion  of 
stability  for  a  long  shaft,  '  simply  supported  '  (i.e., by  journals  which  impose  no  con- 
straint upon  its  direction)  at  each  end,  and  subjected  simultaneously  to  end  thrust 
and  twist.1  His  results  showed  that  the  influence  of  twist  is  relatively  unimportant, 
and  may,  for  most  practical  purposes,  be  neglected  :  the  criterion  then  reduces  to 
Euler's  familiar  formula  for  the  critical  load  of  a  free-ended  shaft. 

When  the  shaft  is  not  subjected  to  twist,  but  revolves,  end  thrust  may  combine 
with  the  rotation  in  bringing  about  instability  of  another  type,  which  reveals  itself 
as  a  tendency  to  '  whirl.'  'The  two  factors  in  this  result  are  far  more  comparable  in 
importance,  and  the  criterion  of  instability  (fortunately  simple)  must  be  used  in  its 
complete  form.2  The  question  then  arises — when  end  thrust,  twist  and  rotation  act 
simultaneously,  as  will  often  occur  in  practice,  are  we  justified  by  GreenhilPs  result  in 
neglecting  the  influence  of  twist,  or  must  we  allow  for  a  tendency  on  the  part  of  the 
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FIG.  17 

rotation  to  combine  with  the  twist  on  more  or  less  equal  terms,  as  it  combined  with 
the  thrust  in  the  second  of  the  problems  just  described  ? 

§2.  An  answer  to  this  question  is  afforded  by  the  accompanying  diagrams. 
Instability  will  occur  when  a  certain  relation  obtains  between  the  thrust  P,  the  rota- 
tional speed  a>  and  the  torque  T,  and  the  stability  criterion  is  exhibited  graphically 
by  means  of  curves  which  connect  critical  values  of  the  first  two  factors  for  a  series 
of  values  of  the  third.  The  contours  given  in  the  figures  are  the  results  of  calcula- 
tion :  additional  contours  may  be  obtained  from  these,  if  required,  by  cross-plotting 
and  interpolation. 

Fig.  17  relates  to  a  shaft  of  which  the  ends  are  *  simply  supported,'  and  fig.  18  to 
an  '  encastr6 '  shaft,  in  which  change  of  direction,  as  well  as  displacement,  is  prevented 

1  Proc.  Tnst.  Meek.  Eng.,  1883,  pp.  182-225.     The  case  of  clamped  ends  was  also 
treated,  with  results  which  do  not  agree  with  the  calculations  of  this  paper :  the  ex- 
planation appears  to  be  that  the  condition  of  zero  displacement  at  either  end  was 
not  realised.    Some  discussion  was  also  given  of  the  general  equations  employed  in 
this  paper,  but  solutions  were  not  obtained. 

2  An  investigation  of  this  problem  is  given  by  A.  Morley, '  Strength  of  Materials,' 
§166.      It  was  briefly  discussed  by   Greenhill  in  p.  208  of  the  paper  referred  to 
Vide  equation  (18)  of  the  present  paper. 
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at  the  ends.  The  co-ordinates  actually  employed  are  non-dimensional  quantities, 
A,  B  and  C,  related  to  the  dimensions  and  material  of  the  shaft,  and  to  P,  «  and  T, 
thus : — 


A-  Tl 
~2ET 


4EI 


and  C  = 


where 


(1) 


I  denotes  the  length  of  the  shaft  between  bearings, 
El  denotes  the  flexural  rigidity  of  the  shaft, 

•\X7 

and    —  denotes  the  mass  of  the  shaft  per  unit  length. 

Any  self -consistent  system  of  units  may  be  employed. 

The  diagrams  have  been  extended  into  the  region  in  which  B  is  negative,  for  the 
purpose  of  illustrating  the  effect  of  end  tension  in  maintaining  stability.  Negative 
values  of  C  are  of  course  inadmissible,  since  they  correspond  to  imaginary  values  of 
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FIG.  18. 

«,  and  the  sign  of  A  (i.e.,  of  the  torque)  is  immaterial.  In  fig.  17  dotted  lines  have 
been  drawn  to  represent,  for  values  of  A  corresponding  to  the  calculated  contours, 
the  contours  which  would  be  given  by  the  formula 

4j5       loC-/       i  xo\ 

-^  +  —^  =  1 (2) 

This  formula  can  be  shown  to  be  exact  when  either  A  or  C  is  zero,  and  the  close 
agreement  between  the  two  sets  of  contours  shows  it  to  be  a  very  close  approximation 
in  all  cases  where  B  is  positive.  A  similar  formula,  with  the  terms  on  the  left  hand 
side  altered  by  the  addition  of  coefficients  (H9, 0-25  and  0-195  respectively,  represents, 
with  rather  less  close  approximation,  the  results  of  fig.  18. 

§3.  Perhaps  the  most  striking  feature  of  the  problem  is  the  relative  complexity  of 
the  analysis  which  leads  to  such  simple  results.  Although  shown  by  actual  com- 
parison to  be  closely  equivalent  to  a  formula  of  the  type  (2),  the  exact  relation 
obtaining  in  either  case  between  A,  B  and  C  is  of  a  highly  transcendental  nature, 
and  its  interpretation  entails  much  laborious  calculation.3  For  practical  purposes, 
the  diagrams  will  be  quite  sufficiently  accurate,  and  more  convenient  than 


3  The  authors'  thanks  are  due  to  Miss  S.  W.  Skan  for  her  valuable  assistance  in 
this  work. 
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algebraic  formulas  of  the  type  (2),  ^besides  applying  over  a  wider  range ;  but 
the  establishment  of  these  formulas  enables  us  to  assert  that  the  relative  importance 
of  end  thrust,  rotation  and  twist  is  of  the  same  order,  when  the  three  actions  exist 
simultaneously,  as  it  has  been  shown  to  be  when  they  are  taken  as  acting  in  pairs : 
that  is  to  say,  any  torque  which  a  shaft  of  normal  dimensions  can  sustain  without 
damage  to  the  material  will  be  negligible  in  its  effect  on  the  stability  of  that  shaft. 

§4.  Equations  of  Neutral  Equilibrium. — In  considering  the  stability  of  the 
system,  we  begin  by  assuming  that  the  shaft  is  initially  straight  and  centrally 
loaded,  but  that  conditions  of  neutral  stability  are  eventually  reached  under  which 
the  shaft  can  be  held  in  equilibrium  in  a  configuration  of  slight  distortion,  the 
centre-line  taking  a  curved  form,  either  plane  or  tortuous.  The  influence  of  gravity 
is  neglected.  We  take  the  origin  of  co-ordinates  at  the  central  section  of  the  shaft, 
the  axis  Oz  coinciding  with  the  axis  of  the  shaft  in  the  unstrained  configuration  ; 
and  we  choose  axes  Ox,  Oy,  which  rotate  with  the  shaft,  perpendicular  to  each  other 
and  to  Oz.  We  define  the  distorted  centre-line  by  the  co-ordinates  x,  y  and  z, 
relative  to  these  axes,  of  a  point  P  on  the  centre-line  which  was  originally  distant  s 
from  the  origin  0. 

The  bending  of  the  shaft  at  any  section  P  is  in  the  osculating  plane,  at  P,  to  the 
curved  centre-line  ;  and  it  is  produced  by  a  couple,  M,  of  which  the  axis  is  parallel 
to  the  bmormal  at  P.  The  direction-cosines  of  the  binormal  are 

dy      dzz  _  dz      d*y  \      x 
dg2      ds  '  ds'1  j 


(dz      dyx  _  dx 
P\ds  '   ~d*      Sl 


(3) 


where  p,  the  radius  of  curvature,  =  El/ M  ; 

and  the  couple  M  can  therefore   be  resolved  into  component   couples,  with   axes 
parallel  to  Ox,  Oy,  Oz  respectively,  of  magnitudes 

dzz  _  dz 
ds~  '  d?      ds 

d*x      dx      d*z 


Now,  in  considering  the  stability  of  the  straight  shaft,  we  are  concerned  only 
with  a  configuration  of  infinitesimal  distortion,  and  we  may  therefore  regard  x  and 
y  as  small  quantities  of  the  first  order.  Neglecting  small  quantities  of  the  second 
order,  we  may  take  dz,ds  as  unity,  and  reduce  the  expressions  for  the  component 
couples  to  the  simpler  forms 


Again,  the  direction-cosines  of  the  tangent  to  the  strained  centre-line  at  the 
point  P  are  (if  we  again  neglect  small  quantities  of  the  second  order) 

dx 
3T 

Thus,  if  Tf  is  the  torsional  couple  acting  on  the  shaft  at  the  section  P,  it  may 
be  resolved  into  components,  with  axes  parallel  to  Oa;,  Oy,  Oz  respectively,  of 
magnitudes 
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The  total  magnitudes  of  the  component  couples  which  act  at  any  section  are, 
therefore, 


,  . 

ds          ds*        ds         ds2 

These  component  couples  represent  actions  produced  by  the  stresses  across  the 
section  P  upon  that  part  of  the  shaft  which  lies  between  0  and  P.  In  addition, 
the  stresses  will  produce  component  forces,  in  the  direction  of  the  axes  Ox,  Oy, 
Oz,  which  we^may  represent  by  X,  Y  and  Z.  Now  the  component  accelerations, 
due  to  rotation,  of  an  element  of  the  shaft  situated  at  P  are 

-Wfti'a?,  _Ww«y  andO, 

9  9 

in  the  directions  O<s,  Oy,  Oz   respectively;    hence,  the  equations   of   equilibrium 
for  the  element  are 

f  +^  =  °>  ......   (*) 


~  =o,  .     .     .     .     .     .  (6) 

Y  +  Zf  =0,     .        .  ,     (7) 

ds 

X-Z^  =  0,  (8) 

ds\     ds         Aft  ds 

and^-X^  +  Y^  =  0,    ,  (9) 

ds          ds         ds 

From  (4)  and  (5)  it  is  evident  that  X  and  Y  are  small  quantities  of  the  first 
order.    Hence,  neglecting  small  quantities  of  the  second  order  in  (9),  we  may  write 

5-ft 

ds 

whence 

T'  =  const.  =  T  (by  end  conditions).         .        .        .        .     (10) 

We  have  similarly,  from  (6), 

Z  ss  const.  =  —  P  (by  end  conditions)  ;  .  .  .  .  (11) 
and  on  substituting  from  (10)  and  (11)  in  (7)  and  (8),  differentiating,  and 
substituting  for  ^  and  -j-  from  (4)  and  (5),  we  obtain,  finally,  as  the  equations  of 
neutral  stability, 

^(T^-EI^U^-P^O  ......  (12) 

ds*\    ds         d&)     g  ds2 

and 


s        g  s 

These  two  equations  replace  the  six  equations  (4)-(9). 

S6.  Special  Case  of  Zero  Torque.—  When  T  is  zero,  a  solution  of  (13)  may  be 
obtained  by  assuming  that  a>  =  0  identically.     Equation  (12)  then  reduces  to 


and  a  solution  is  obtainable  in  the  form 

^sK,  sinew  +K2  cos  <w+K,sinh/3*  +  K4  cosh  /3«,    .        .        .    (15) 
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where  K,,  K2,  K3,  K4  are  arbitrary  constants,  and  where 

.W««.  _P_      \ 
*KI  +  2EI' 

P 


=  v/  / 
'  V 


4E312      ^EI     2EI'     * 

The  values  of  K,,  K2,  K3,  K4  —  or  rather,  their  relative,  magnitudes—  are  determinabie 
from  the  end  conditions  of  the  problem.  But  since  they  cannot  all  vanish 
simultaneously  (for  the  shaft  would  then  remain  straight),  we  have  only  three 
ratios  at  our  choice  for  satisfying  four  specified  conditions  ;  hence,  by  elimination  of 
the  four  constants,  we  obtain  the  criterion  of  neutral  stability. 

It  will  be  convenient  in  this  section  of  the  paper  to  take  our  origin  at  one  end. 
Then,  when  the  ends  are  simply  supported  we  have  the  relations 


and  substituting  from  (15)  we  find  that 

K2  =  K3  =  K4  =  0, 
K,  sin  al  =  O. 

It  follows  that  we  must  have 

sinaZ=--0  ......    .        .     (17) 

and,  by  (16),  the  criterion  for  neutral  stability  in  this  case  is 


Z2    2EI      V  4E2I2 
PZ2 


When  the  ends  are  clamped,  so  that  we  have  the  relations 

^  =  ^  =  0, 
as 

we  find,  on  substitution  from  (15),  that 


=      =  0,  when  s  =  O  or  I, 
as 


so  that 

y  =  K,  (sin  as  —  5  sinh  &s)  +  K^  (cos  as — cosh  £*). 


The  criterion  then  takes  the  form 

Z,      cos  al  —  cosh  £Z, 

=  0, 


sin  aZ  — ?  sinh  £Z,      cos  aZ  — cosh  £ 
P 

cos  aZ  —  cosh  £Z,  —  (  sin  aZ  +  _  sinh 
\  a 

and  on  expanding  the  determinant  we  obtain,  finally,  as  the  criterion  for  neutral 
stability, 

2(1  -  cos  al  cosh  /3Z)  =  °^H^2  sin  al  sinh 
a  and  p  being  given  by  (16),  so  that  (19) 


where  B  and  C  are  the  quantities  defined  in  equation  (1). 
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§6.  General  Case.—  Reverting  to  equations  (12)  and  (13),  and  eliminating  x  or  y, 
we  obtain  the  equation 


The  solution  may  be  written  in  the  form 

(or,  y)  =  2[A,  sin  A,s  +  B,  cos  Ar<?] 
where  A,,  A2,  .      .  .  etc.,  are  the  roots  of  the  equation 

(EIA4  -  PA2  -  ^0,2  V  1  T2A6  _  O 
V  9      ' 

Writing  this  equation  in  the  form 

PA2-W«2  =  ±TA3,  • 


we  see  that  the  roots  obtained  by  taking  the  negative  sign  on  the  right-hand  side 
are  equal  and  opposite  in  sign  to  those  which  are  obtained  by  taking  the  positive 
sign.  The  most  general  expression  for  x  is  therefore  given  by 

x  =  A,  sin  A,*  +  A2  sin  \^s  +  A,,  sin  \3s  +  A4  sin  A4s 

+  B,  cos  A!«  +  B2  cos  AjS  +  B3  cos  \3s  +  B4  cos  A4$,        .        .        .     (21) 

where  A,,  A.J,  A*,,  A4  are  roots  of  the  equation 

W 
%.;  EIA4-TA3-PA2-   v-»2=0;  .         ...     (22) 

and  the  corresponding  expression  for  y  can  be  obtained  from  (13)  in  the  form 

y  —  B!  sin  \s  +  B2  sin  A-jS  +  B3  sin  \3s  +  B4  sin  \ts 

—  A,  cos  AjS  —  A.J  cos  A^  —  A3  cos  A^  —  A4  cos  A4«.         ,         .         .     (23) 

These  expressions  are  legitimate  whether  A,,  A.2,  A3,  A4  be  real,  imaginary,  or 
complex. 

Tli  e  Criterion  for  'Simply  Supported'  Ends.—Qu  the  assumptions  (a)  that  the 
journals  prevent  movement  of  the  ends  of  the  shaft,  but  without  constraining  the 
dirqctwn  of  the  bent  centre-line,  and  (5)  that  the  twisting  couple  applied  at  the 
ends  of  the  shaft  has  its  axis  parallel  to  Oz,  the  terminal  conditions  may  be  expressed 
as  follows  :  — 


(24) 


The  last  two  equations  represent  the  condition  that  at  the  ends  the  resultant  couple 
on  the  shaft  has  Oz  for  its  axis. 

If  we  substitute  for  x  and  y  from  (21)  and  (23),  these  conditions  yield  the 
relations 

A;  sin  A,     +  A2  sin  A2  L  +  A3  sin  A3  _  +  A4  sin  A4  -  =  O,  .         .         .     (25) 
A,  cos  A,  -  +  A2  cos  A2     +  A3  cos  A3  -9  +  A  ,  cos  A4  -  =  O,  .         .         .     (26) 


.     (27) 


(EIA,2  —  TAi)  A,  sin  A,  i.  +  (ElA,,2 —  TA2)  A2  sin  A2  - 
2  2 
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+  (El  V  -  TAO  A3  cos  A3     +  (El  V -  TA4)  A,  cos  \4    =  O,  . 


.     (28) 


and  a  similar  set  of  four  relations,  which  may  be  obtained  by  writing  B,,  B2,  .... 
for  A,,  A2,  ....  in  the  above.  Either  set  of  four  relations  obviously  yields  the 
same  equation  in  I,  if  we  eliminate  the  coefficients  Ap  A^,  ...  or  Bp  B2,  .... 

Subtracting  (25),  multiplied  by  P,  from  (27),  and  remembering  that,  by  (22), 

when  A  has  any  of  the  values  A,,  A2,  A3,  A4,  we  see  that  for  (27)  we  may  substitute 
the  relation 

l-  sin  A,  -  +  ^  sin  A2  -  +  ^  sin  A3  -  +  ^  sin  A4  -  =  O.    .        .     (30) 


Similarly,  for  (28)  we  may  substitute  the  relation 


.(31) 


Then  from  (25),  (26),  (30)  and  (31),  eliminating  A,,  A2,  Ag,  A4,  we  have,  as  the 
required  criterion  for  a  condition  of  neutral  stability  in  the  free-ended  shaft, 


sin  A,  -, 

1 

sin  A3  --, 

sin  A4  -g 

1 

1 

1 

1 

COS  A,     » 

cos  A2-» 

cos  A3  -, 

COSA4_2, 

,                  1 

1                j 

1              1 

1       .            1 

v  sin  AI  2' 

V  Sln  *2  2 

-   sm  A,    , 
V             2 

V  si:  L  A<  2' 

^COSA,^ 

1               lf 

1          A  l^ 

1  COSA  2, 

,  V           '2 

A/           22 

V           J  2 

A42            42 

=  o. 


The  expansion  of  the  determinant  on  the  left-hand  side  of  (32)  is 


cosecx, 


_x4t)  sin  Aj  *  8in  (Ai-Af)  I  sin  (Ai_X4)  « 


-(V-V)  (V-  V)  sin  A,     sin  (A.-A,)     sin  (A,~A4)  l. 


^-A,*)  sin  A4     sin  (\- 


and  since 

(V-V) 

and  the  quantity 


-O  (V-V), 


(32) 


—which,  by  (22),  is  equal  to   /  — cannot,  by  hypothesis,  be  zero,  it   is  easily 
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shown  that  the  condition  (32)  may  be  written  in  the  form 


3-4 

a 

-(V-V)  (V-V)   COS  (A2  +  A^A.-A,)! 

-  (  V  ~  V)  (  V-  V)   COS  (A,  +  A,  -  A2  -  X4)?  =  0, 


(33) 


which   may   therefore  be  taken  as   the  criterion  of  neutral  stability  for  'simply 
supported '  ends. 

The  Criterion  for  '•Clamped  Ends.1 — On  the  assumption  that  the  bearings  prevent 
both  displacement  and  change  of  direction,  the  terminal  conditions  take  the  form 


.     (34) 


Substituting  for  x  and  y  in  these  equations  from  (21)  and  (23),  we  obtain  the 
relations 

A,  sin  A,  -  +  A2  sin  AS  -  +  A3  sin  A3  -  +  A4  sin  A4  -  =  0,       .         .     (35) 


2 


A,  cos  A,  -  +  A2cosA2  -  +  A3cosAs  2 


.     (36) 


A,  A,  sin  A,  -  +  A.jA2  sin  A2  -  +  A3A8  sin  A3  -  +  A4A4  sin  A4  -  =  O,   .     (37) 


4-  =  0,   .     (38) 


and  a  similar  set  of  four  relations,  which  may  be  obtained  by  writing  B,,  B2,  . 
for  A,,  A2,  .  .  .  in  the  above.     Either  set  of  four  relations  yields  the  same  equation 
in  I,  if  we   eliminate  the  coefficients  A,,   A2,  .  .  .  or   B,,   B2,  .  .  .  etc.,   and  the 
condition  for  neutral  stability  may  therefore  be  written  in  the  form 


7 

I 

^ 

I 

in  *,  4- 

sin  A,  -, 

sin  As  -, 

sin  A4-^, 

7 

7 

^ 

^ 

COS  A,  - 

COS  A2  ~9 

COS  A,-, 

COS  A4  -, 

I 

I 

I 

I 

A,  sin  A,  -, 

A2  sin  A2  ^, 

A,  sin  A3-, 

A4sinA4-, 

I 

7 

/ 

I 

A,  COS  A,     , 

A.J  COS  A2-, 

A8COS  A3-, 

A4  COS  A4  -, 

O. 


(39) 


Expanding  the  determinant,  we  obtain,  as  the  criterion  in  this   instance,  the 
equation 

(A,  -  A2)  (A,-  A4)  COS  (A,  +  A2-  A3-A4)  - 


-(^1-^3)  (^2-^4)  COS 

in  which  A,,  A.,,  A3,  A4  are,  as  before,  the  roots  of  equation  (22). 


.      (40) 
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§7.  Equal  Roots  of  Equation  (22)  not  Admissible.  —  At  first  sight  it  would  appear 
that  the  conditions  (33)  and  (40)  of  neutral  stability  would  both  be  satisfied  if 
equation  (22)  had  equal  roots  :  and  real  values  can  be  found  for  P,  T  and  «2  which 
will  satisfy  this  requirement.  But  if  we  merely  substitute  A,  for  A2  (say)  in  the 
expressions  (21)  and  (23)  for  x  and  y,  we  shall  in  effect  be  reducing  the  number  of 
arbitrary  constants  involved  in  these  expressions  from  eight  to  six,  and  it  is  not 
difficult  to  show  that  the  boundary  conditions  can  only  be  satisfied  by  making  a?  and 
y  vanish  for  all  values  of  s  :  the  shaft  then  remains  straight,  and,  although  this  is  a 
configuration  of  equilibrium,  it  is  of  no  interest  for  present  purposes.  If,  on  the 
other  hand,  we  write  down  the  complete  solution  (with  eight  constants)  under  the 
condition  of  equal  roots,  we  are  in  effect  substituting  A,  for  \z  in  expressions  for  x 
and  y  which  differ  from  (21)  and  (23)  in  having  sin  \^  and  cos  A^  replaced  by  s  cos 
A^-  and  s  sin  A^  respectively  ;  but  when  this  latter  modification  is  made  in  (21)  and 
(23),  and  the  criterion  obtained  by  the  same  methods  as  before,  it  is  found  to  be  no 
longer  satisfied  by  making  A.J  equal  to  Ar 

§8,  Solution  for  '  Simply  -supported  '  Ends.  —  We  proceed  to  interpret  the  solutions 
(33)  and  (40)  in  terms  of  the  physical  constants  of  our  problem.  Writing 


for  —  ,  we  may  throw  equation  (22)  into  the  form 


(41) 


where  A,  B  and  C  are  the  quantities  defined  in  (1)  ;  and  we  may  notice  that  two  of 
the  roots  are  necessarily  real,  since  the  constant  term  in  this  equation  is  negative. 
If  we  express  the  roots  in  the  form 


.     (42) 

3,     4  , 

equation  (41)  must  factorise  as  follows  :  — 

{/t2 
and  we  have 


Again,  the  criterion  (33)  may  be  written  as 

Oi2  ~  M2a)  (X.2  -  M42)  sin  fa,  -  /*,)  sin  (^  -  /*4)    1  f.. 

'(Mi1-*1)  Os'-M*')  sin  (M,-M2)  sin  Os  -M«),  J 

and  if  we  substitute  from  (42)  this  becomes 

IQabcd  sin  (a  +  b—  c—  d)  sin  (a—  b  -c  +  d) 

=  (a2  +  &2  -  c2  -  d*  +  2ab  -  2cd)  (a?  +  b*-c*-d?-  2ab  +  2cd)  sin  2*  sin  2d, 
=  {  Sabcd  +  (a?  -  &2)8  +  (c2-  d2)2  -  2(a2  +  ft*)  (6*  +  rf2)}  sin  2b  sin  2d, 


Sabcd  {cos  2b  cos  2d  -cos  2(a-c)\ 
=  {  (a?  -  ft2)8  +  (c*  -  d*)*  -  2(02  +  ft»)  (c2  +  d*)  }  sin  2b  sin  2d.      .        .     (46) 

In  this  form,  we  can  obviously  deal  with  either  real  or  imaginary  values  of  b  or  d, 
observing  that,  if  6  =  0i, 

sin  26     sinh  2)8 

26  2/3 

and 

cos  26  =  cosh  2)3. 


354  REPORTS   ON   THE  STATE   OF  SCIENCE,   ETC. 

The  solution  of  (45)  has  been  obtained  by  trial,  and  is  exhibited  in  tig.  17.  The 
procedure  adopted  is  to  trace  a  contour  line  (T  =  const.)  on  a  diagram  of  which  the 
co-ordinates  are  P  and  o>2.  We  are  then,  by  (43),  given  the  value  of  (a  +  c),  but  can, 
with  this  limitation,  vary  the  separate  values  :  choosing  any  pair  of  values,  we  have 
from  the  third  of  (43)  a  range  of  corresponding  values  of  &2  and  d2.  The  sign  is 
immaterial,  and  for  each  pair  we  calculate  the  left-  and  right-hand  expressions  in 
(45).  Proceeding  in  this  way  by  trial,  when  the  expressions  are  equal  we  know 
«,  J,  0,  d,  and  hence  B  and  C,  from  (43) 

It  may  be  noticed  that  a  possible  solution  of  the  criterion  (45)  is  given  by  a  =  c  -  0, 
2  b  =  rnr.  This  gives 

A  =  0, 

B  =  #  +  **, 


and  hence 

4B          J6C 


PI2         Wa>2Z< 


which  reduces  to  (18)  when  »=l.     When  T  is  small,  therefore,  we  may  expect  a 
solution  round  about  the  value  2b  =  nir,  and  for  practical  purposes  n  will  be  1. 

§9.  Solution  for    Clamped  Ends.  —  Employing   as   before    the   substitutions   of 

72 

(41)-(43),  -  and  multiplying  (40)  throughout  by      ,  we  may  write  this  criterion  in 

the  form 

(Mi  -  M2)  (M3  ~  M4)  cos  (/*,  +  M2  ~  M3  -  M4) 

~  (Mi  -  M3)  (M2  ~  M4)  COS  (M3  +  M!  -  M2  -  M4) 
=  (M2-M3)  (/t4-Mi)cos(/x2  +  )u3-/ti4-/Li1), 

=  1  (Mi  -  M2)  (M3  -  M4)  -  (Mi  -  Ms)  (M2  -  M4)  }  cos  (^  +  ^  -  jt4  -  /*,), 
or 

(M,  -  M2)  0*s  ~  M4)  sin  (M,  -  M3)  sin  (M2  -  M») 

/x4),          .         .         .     (46) 


which  corresponds  with  equation  (44)  of  the  last  section.     Hence,  if  we  substitute 
from  (42),  we  obtain  the  relation 


kbd  sin  (a+l>  —  c  —  d)  sin  (a  —  b  — 

=  (a  +  b-c-d)  (a  —  b  —  c  +  d)  sin  2bsin  2d, 
or 

2bd  {  cos  2(b  -  d)  -  cos  2(a  -  c)  \ 

=  {  (a  -c)2  -(b-d)*}  sin  2&sin  2d, 
which  may  also  be  written  in  the  form 

2bd  {cos  2b  cos  2d  -cos  2(a-c)} 

....     (47) 


for  convenience  in  dealing  with  imaginary  values  of  b  or  d. 

Equation  (47)  has  been  solved,  like  (45),  by  trial,  and  the  results  are  exhibited 
in  fig.  18. 
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IX. 

The  Stresses  in  Cylinders  and  Pipes  with  Eccentric  Bore. 
By  G.  B.  JEFPERY,  M.A.,  D.Sc. 

In  a  recent  paper1  the  author  developed  a  general  method  for  the  solution  of 
two-dimensional  elastic  problems  in  which  the  stresses  are  given  over  two  circular 
non-intersecting,  non-concentric  boundaries.  One  of  the  problems  which  most 
readily  yields  to  this  method  of  solution  is  that  of  a  cylinder  or  pipe,  whose  internal 
and  external  sections  are  both  circular  but  not  concentric,  in  stress  under  a  uniform 
hydrostatic  pressure  either  internal  or  external.  The  problem  is  soluble  in  finite 
terms  and  reference  may  be  made  to  the  original  paper  for  a  discussion  of  the  stress 
distribution.  These  results  may  be  applied  to  the  determination  of  the  diminution 
in  the  strength  of  a  pipe  or  cylinder  due  to  a  small  error  of  centring  as  between  the 
external  surface  and  the  bore.  If  the  radius  of  the  bore  is  rv  and  that  of  the 
external  surface  rz  and  if  the  distance  between  their  centres  is  d  the  maximum 
stress  when  the  pipe  is  under  an  internal  pressure  P  is  on  the  internal  surface  at 
the  thinnest  part  if  d  <  %r{  and  is  of  magnitude. 


f 


This  depends  only  upon  the  ratios  of  rv  r2,  d,  and  we  may  therefore  take  a  =  ratio  of 
mean  thickness  of  cylinder  wall  to  internal  diameter,  &  =  ratio  of  centre  distance 
to  mean  thickness  of  cylinder  wall,  or 


The  maximum  stress  them  becomes 

P   "2(1  +  2a)2  {  (1  +  2a)2  +  1  -  4a)8  - 


In  the  appended  tables  the  maximum  stress  in  the  material  is  given  in  Ibs.  per 
sq.  in.  per  1000  Ibs.  per  sq.  in.  internal  pressure  for  values  of  o,  0  lying  between  0 
and  -2  in  each  case. 

1  Plane  Stress  and.  Plane  Strain  in  Bi-polar  Co-ordinates.     Phil.  Trans.  Royal 
Society,  Vol.  221  A,  p.  265. 
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